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The  goal  of  this  dissertation  was  to  experimentally  investigate  the  high  rate 
and  high  pressure  mechanical  response  of  a mortar  and  concrete  mix  and  use  or 
develop  a constitutive  model  able  to  describe  the  observed  behavior.  Triaxial 
compression  tests  at  a strain  rate  of  10~6/s,  and  for  confining  pressures  ranging 
from  0 to  0.5  GPa  were  conducted.  Dynamic  tests  in  the  range  60/s  to  about 
160/s  under  both  unconfined  and  confined  conditions  were  conducted  using  the 
University  of  Florida’s  7.62  cm  diameter  split  Hopkinson  pressure  bar  (SHPB). 

The  data  obtained  in  such  tests  allowed  the  quantification  of  the  combined  effects 
of  confining  pressure  and  strain  rate  on  the  deformation  and  strength  of  the 
materials.  For  mortar,  dilatancy  has  been  observed  at  high  levels  of  the  principal 
stress  difference  for  both  dynamic  and  quasi-static  conditions.  The  unconfined 
dynamic  compressive  strengths  are  approximately  double  those  of  the  quasi-static 
compressive  strengths.  Most  of  the  confined  SHPB  mortar  specimens  showed  very 
little  damage  post-test  other  than  some  chipping  around  the  top  edges,  most  likely 
due  to  localized  tensile  effects. 
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For  the  concrete  selected  for  this  research,  WES5000,  quasi-static  hydrostatic 
tests  conducted  up  to  a pressure  of  0.5  GPa  allowed  for  the  accurate  determination 
of  the  dependence  of  the  bulk  modulus  on  pressure  and  the  correct  estimation  of 
the  material’s  compaction  properties  when  subjected  to  pressures  in  the  range 
encountered  in  dynamic  events.  For  confined  quasi-static  conditions,  the  material 
exhibited  hardening  behavior  up  to  failure.  Both  compressibility  and  dilatancy 
regimes  of  the  volumetric  behavior  were  observed,  the  dilatancy  threshold  being 
highly  dependent  on  the  level  of  confinement.  The  unconfined  dynamic  strength 
is  as  high  as  1.5  times  the  quasi-static  strength,  the  material  generally  exhibiting 
far  more  cracking  under  similar  loading  conditions  than  was  observed  in  mortar. 
The  confined  dynamic  tests  showed  similar  stress-strain  response  as  the  quasi- 
static tests  conducted  at  the  same  level  of  confinement.  A decrease  in  strain  rate 
sensitivity  with  increasing  confining  pressure  was  observed. 

A new  elastic/ viscoplastic  model  that  captures  compressibility  and  dila- 
tancy, as  well  as  strain  rate  effects  has  been  developed  for  concrete.  As  a general 
framework,  the  elastic-viscoplastic  formulation  of  Cristescu  was  chosen  due  to  its 
apparent  capability  to  capture  the  behavior  of  interest  and  the  fact  that  there 
are  no  a priori  limitations  or  restrictions  regarding  the  specific  expressions  ot  the 
yield  function  and  viscoplastic  potential.  However,  it  was  found  that  in  order  to 
capture  the  peculiarities  of  the  volumetric  response,  namely  the  reduced  dilatancy 
at  very  high  pressure,  the  hypothesis  of  existence  of  a viscoplastic  potential  has 
to  be  abandoned.  A new  flow  rule  was  proposed  for  WES5000.  Procedures  for 
determination  of  the  constitutive  functions  were  developed.  Comparison  between 
the  model  predictions  and  data  showed  that  the  proposed  model  describes  with 
very  good  accuracy  the  high-pressure  behavior  of  concrete. 


xvi 


CHAPTER  1 
INTRODUCTION 

The  past  century  has  seen  tremendous  advances  in  the  use  of  concrete  as  well 
as  numerical  techniques  for  the  design  and  analysis  of  concrete  structures.  It  is 
reasonable  to  say  that  it  has  become  the  most  ubiquitous  of  modern  construction 
materials,  used  in  everything  from  roads  and  bridges  to  modern  office  structures 
and  massive  hardened  military  command  and  control  facilities  designed  to  survive 
blast  and  impact  loadings.  Yet,  it  is  perhaps  the  least  understood  of  construction 
materials.  Concrete  is  a composite  material  made  up  of  fine  aggregate  (sand), 
coarse  aggregate  (rock),  Portland  cement  and  water.  On  a micro-scale,  it  is  a 
heterogeneous  and  anisotropic  material.  Yet  on  a macro  scale,  it  can  be  treated 
as  a continuum.  It  has  been  reported  that  concrete  exhibits  both  compressible 
and  dilatant  properties,  pressure  dependencies,  and  highly  rate  sensitive  strength 
and  failure  properties  in  both  tension  and  compression.  The  rate  effects  are 
truly  pronounced  in  the  regimes  typical  of  impact  and  blast  loading  problems. 

It  is  in  this  highly  dynamic  realm  that  the  need  exists  for  improved  concrete 
models.  A material  constitutive  model  needs  to  be  capable  of  replicating  the 
behavior  of  interest  and  be  based  on  data  gathered  in  stress  and  strain  rate  regimes 
representative  of  in  situ  material.  The  strain  rates  of  interest  for  various  classes 
of  impact  problems  were  tabulated  by  Zukas  et  al.  [1]  and  are  presented  in  Table 
1-1.  Wilbeck  [2]  reported  that  the  impact  pressures  commonly  seen  for  explosively 
driven  problems  or  kinetic  energy  penetrators  range  from  on  the  order  of  the 
strength  of  the  target  material  to  the  strength  of  the  projectile  respectively.  This 
implies  that  the  pressure  range  of  interest  varies  from  tens  of  Mega-Pascals  (MPa) 
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to  the  order  of  a Giga-Pascal  (GPa)  for  a case  involving  a high-strength  steel 
penetrator. 

Table  1-1:  Strain  rate  response  as  a function  of  impact  velocity 


Striking  Velocity,  km/s 

Strain  Rate,  mm/mm/s 

> 12 

108 

3tol2 

107 

lto3 

105 

0.5tol 

104 

0.05fo0.5 

102 

< 0.05 

10° 

This  program  of  research  is  composed  of  concomitant  efforts  of  theoretical 
analysis,  modeling,  and  experimental  test,  the  main  goals  of  which  are 

1)  To  characterize  the  behavior  of  WES5000  concrete  over  a broad  range  of 
stress  and  strain  rate  conditions. 

2)  To  develop  a model  that  could  capture  the  main  features  of  the 
observed  behavior. 

The  remainder  of  this  document  is  outlined  as  follows.  Chapter  2 presents 
a review  of  testing  devices  and  experimental  techniques  used  to  characterize  the 
quasi-static  (strain  rate  ~ 10_6/s)  behavior  of  geologic  materials  and  geologically 
derived  materials  such  as  concrete.  The  concrete  material  used  in  the  program 
is  then  presented  along  with  experimental  results  under  unconfined  and  confined 
conditions  in  the  range  of  50  to  500  MPa. 

Chapter  3 discusses  dynamic  material  characterization.  A general  background 
on  dynamic  material  characterization  is  provided  followed  by  a description  of 
the  confined  dynamic  Split  Hopkinson  Pressure  Bar  technique  developed  for  this 
research  effort.  Results  are  presented  for  strain  rates  on  the  order  of  lO1^  to  102/s 
under  both  confined  and  unconfined  conditions. 

In  Chapter  4,  after  reviewing  different  approaches  to  concrete/geologic 
material  constitutive  modeling,  an  elastic/ viscoplastic  formulation  is  proposed. 
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The  structure  of  the  model  is  presented,  followed  by  derivation  of  the  constitutive 
functions,  and  determination  of  the  model  parameters  based  on  experimental 
results  generated  in  this  effort.  Finally,  comparisons  are  made  between  theoretical 
predictions  and  experimental  data.  The  summary  and  conclusions  are  presented  in 
Chapter  5. 


CHAPTER  2 

EXPERIMENTAL  INVESTIGATION  OF  THE  QUASI-STATIC  BEHAVIOR 

2.1  Background  and  Overview  of  Testing  Devices  and  Quasi-Static  Experimental 

Techniques 

To  fully  understand  the  mechanical  behavior  of  a material,  tests  must  be 
conducted  that  correspond  to  the  range  of  pressures,  strains,  and  strain  rates 
encountered  in  the  applications  of  interest.  For  dynamic  events  such  as  projectile 
penetration,  pressures  may  vary  from  the  gigapascal  range  in  close  proximity  to  the 
nose,  to  near  zero  at  nearby  free  surfaces  [3],  Likewise,  strain  rates  may  vary  from 
106  s-1  near  the  event  to  virtually  quasi-static  response  many  penetrator  diameters 
away.  There  is  no  one  mechanical  testing  device  with  which  it  is  possible  to 
achieve  this  broad  range  of  pressures  and  strain  rates.  Figure  2-1  shows  the  typical 
mechanical  testing  devices  and  the  strain  rates  at  which  they  are  applicable.  It  is 
reproduced  from  a work  by  Lindholm  (U.  S.  Lindholm,  High  Strain  Rate  Testing, 

In  R.  F.  Bunshah,  editor,  Techniques  in  Metals  Research,  volume  5.  Interscience, 
New  York,  1971  [4]).  Creep  to  quasi-static  response  characterization  is  generally 
accomplished  through  use  of  a load  frame  device,  frequently  hydraulically  driven. 
Dynamic  response  is  investigated  using  Hopkinson  bars,  split  Hopkinson  pressure 
bars,  or  flyer  plate  experiments.  These  testing  devices  and  techniques  are  presented 
in  Chapter  3.  To  fully  characterize  the  quasi-static  mechanical  behavior,  tests  that 
induce  a general  3-D  stress  or  strain  in  the  material  specimen  are  required.  Such 
tests  require  the  use  of  some  type  of  load  cell.  This  may  take  the  form  of  a cube 
cell.  In  a cube  cell  (sometimes  referred  to  as  a “true’’  triaxial  cell),  the  specimen 
is  held  between  6 platens.  The  advantage  to  this  type  of  test  is  that  the  load  in 
all  three  principal  directions  may  be  controlled  independently.  The  experimental 
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apparatus  is  by  necessity  more  complex  and  costly.  Specimen  fabrication  may 
also  be  more  complicated  and  the  pressure  levels  obtainable  are  generally  below 
that  of  interest  for  penetration  problems.  The  most  commonly  used  device  is  the 
Karman  cell,  named  after  Theodore  von  Karman,  who  in  1911  conducted  triaxial 
compression  experiments  on  marble  [5].  In  a Karman  experiment,  a cylindrical 
specimen  is  held  between  two  platens  of  a load  frame  and  placed  in  a hydrostatic 
confinement  cell. 

A typical  experiment  consists  of  two  phases,  hydrostatic  and  deviatoric.  In 
the  hydrostatic  phase,  the  specimen  is  subjected  to  an  all-around  pressure  which  is 
gradually  increased  up  to  a predetermined  level.  The  confining  pressure  is  applied 
via  a hydraulic  fluid.  In  the  deviatoric  phase,  the  radial  pressure  is  held  constant 
while  an  axial  load  is  applied  at  a constant  rate  to  the  specimen.  The  axial  load, 
applied  via  the  platens,  is  increased  until  the  specimen  fails.  Measurements 
typically  recorded  include  axial  and  radial  strains,  confining  pressure,  and  axial 
load.  With  modern  electronic  servo-controls,  it  is  possible  to  closely  monitor  and 
control  the  load,  loading  rate  and  the  confining  pressure.  The  experiment  may  be 
run  under  load  or  strain  control  and  may  be  run  monotonically  or  with  creep  or 
relaxation  cycles.  Figure  2 2,  taken  from  Chinn  and  Zimmerman  [6],  displays  a 
cut-away  view  of  a typical  Karman  device. 

By  performing  a series  of  experiments  at  several  confining  pressures,  it  is 
possible  to  study  the  influence  of  the  maximum  and  minimum  principal  stresses 
on  material  behavior.  Based  on  data  from  such  tests,  continuum  models  can  be 
developed  and  further  used  to  predict  material  behavior  for  loading  conditions 
which  are  not  accessible  through  direct  testing.  Detailed  discussions  on  how  to 
interpret  the  data  and  evaluate  the  constitutive  parameters  are  well  documented  in 
the  literature  [7,  8,  9].  Hydrostatic  pressure  in  such  experiments  is  limited  by  both 
material  considerations  and  the  design  restrictions  of  the  particular  Karman  cell 
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Figure  2 -1:  Dynamic  aspects  of  mechanical  testing 


utilized.  Karman  cells  have  been  designed  to  accommodate  hydrostatic  pressures 
on  the  order  of  one  gigapascal  (GPa)  [10]. 

Higher  confining  pressures,  beyond  1 GPa,  can  be  achieved  through  the  use  of 
a uniaxial  strain  compression  device.  Such  a device  is  essentially  a rigid  cylindrical 
collar  in  which  the  specimen  is  compressed.  Uniaxial  strain  is  assumed  and  only 
the  axial  strain  and  applied  load  are  measured.  Descriptions  of  such  devices  are 
given  by  Bridgman  11]  and  Cristescu  [12].  Drawbacks  to  this  method  include 
the  exacting  tolerances  required  of  test  specimens,  often  difficult  to  achieve  with 
geologic  material’s,  and  the  difficulty  in  quantifying  the  effects  of  friction  in  the 
experiment. 

In  this  research  effort,  Karman  experiments  were  used  to  characterize  the 
quasi-static  behavior  since  this  technique  offered  the  possibility  of  testing  over  a 
broad  range  of  confining  pressures  while  allowing  the  measurement  of  both  axial 
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Figure  2-2:  Section  through  triaxial  Karman  device 

and  radial  strains,  and  hence  characterization  of  the  volumetric  behavior.  The 
remainder  of  this  chapter  is  organized  as  follows. 

Section  2.1  describes  the  concrete  material  (WES5000)  under  investigation 
and  the  experimental  plan.  Sections  2.2  to  2.11  detail  the  testing  protocol  and  the 
test  results.  Section  2.12  deals  with  evaluation  of  the  elastic  properties  from  the 
unload-reload  cycles.  Finally,  Section  2.13  summarizes  the  quasi-static  material 
characterization  effort. 

2.2  Material  Description 

2.2.1  Historical  Note 

The  first  use  of  what  might  be  called  concrete  was  by  the  the  early  Greeks  who 
learned  to  combine  lime  and  water  with  sand,  crushed  stone  or  brick  and  broken 
tiles.  Archeological  digs  have  shown  that  the  Greeks  utilized  highly  durable  mixes 
as  early  as  1700  B.C.  [13].  The  Romans  later  learned  to  add  a volcanic  ash  to  the 
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mix.  This  ash,  called  pozzolan  after  the  town  of  Pozzuoli  Italy,  gave  the  concrete 
added  durability  and  strength.  This  volcanic  ash  is  for  all  intents  and  purposes  the 
same  as  the  fly  ash  often  added  to  modern  concrete  mixes.  The  first  major  use  of 
concrete  in  large  scale  construction  projects  occurred  in  Rome  approximately  2000 
years  ago.  Perhaps  the  greatest  example  of  the  versatility  and  durability  of  Roman 
concrete  is  the  Pantheon,  constructed  by  the  Roman  Emperor  Hadrian  during  the 
period  118  to  128  A.D..  It  still  stands  in  Rome  and  is  used  as  a church  to  this  day. 
The  Pantheon  is  truly  a marvel  of  ancient  engineering  and  architecture.  In  fact, 
the  142.5  foot  span  of  the  Pantheon  dome  was  not  surpassed  until  the  19th  century 

[14]- 

The  secret  of  Roman  concrete  was  lost  in  410  A.D.  when  Alaric  the  Goth 
sacked  Rome.  It  remained  essentially  lost  until  1824,  when  Joseph  Aspdin  received 
his  patent  for  Portland  cement.  Aspdin  heated  a mixture  of  finely  ground  clay  and 
limestone  until  C02  had  been  driven  off,  but  at  a temperature  much  lower  than 
that  required  for  clinker ing.  The  prototype  of  modern  cement  did  not  appear  until 
1845  when  Isaac  Johnson  burnt  a mixture  of  clay  and  chalk  until  clinkering,  so  that 
the  reactions  necessary  for  the  formation  of  strongly  cementitious  compounds  took 
place  [15].  Joseph  Monier  deserves  the  credit  for  making  the  first  practical  use  of 
reinforced  concrete  in  1867,  while  the  first  cast-in-place  concrete  structures  in  the 
United  States  began  appearing  in  the  1870s  [16]. 

The  French  development  of  the  explosive  shell  from  1885-90  led  to  the  intro- 
duction of  concrete  into  military  construction  as  traditional  masonry  fortifications 
proved  no  match  for  the  new  offensive  technology  [17].  World  War  II  saw  the  tran- 
sition away  from  forwardly  positioned  defensive  fortifications  such  as  the  French 
Maginot  line,  designed  to  counter  a land  based  adversary,  towards  the  use  of  hard- 
ened command  control  and  communications  facilities  designed  for  defense  against 
air  delivered  munitions.  Since  that  time,  the  proliferation  of  concrete  bunkers 
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has  spread  worldwide  as  nations  have  sought  to  protect  vital  national  assets  from 
conventional  and  nuclear  threats.  Today,  the  targeting  of  civilian  structures  by 
terrorist  groups  as  well  as  an  increased  emphasis  worldwide  on  building  more  earth 
quake  resistant  structures  has  made  the  understanding  of  construction  materials, 
such  as  concrete,  to  dynamic  loads  far  more  than  a matter  of  purely  military 
interest. 

2.2.2  Material  Studied,  WES5000 

The  concrete  selected  for  study  is  referred  to  as  WES5000.  The  aggregate 
size  is  such  that  it  must  pass  through  a 3/8  inch  (9.5  mm)  sieve.  The  5000  stands 
for  5000  psi  (34.5  MPa)  unconfined  compressive  strength,  the  original  design 
specification.  In  reality,  the  90  day  unconfined  compressive  strength  is  nominally 
7000  psi  (48.3  MPa).  The  ingredients  and  proportions  for  the  mix  are  given  in 
Table  2-1 


Table  2-1:  Ingredients  and  mixture  proportions  for  WES5000  concrete. 


Item 

Mixture  proportions, 
saturated  surface-dry 

Type  I Portland  cement 

264.0  kg/m a 

Flyash 

55.8  kg/m 3 

9.5  mm  local  unprocessed  chert  coarse  aggregate 

1037.6  kg/m 3 

Local  unprocessed  chert  fine  aggregate 

840.7  kg/m 3 

Water 

145.9  kg/m6 

Water  reducing  admixture  “300N” 

0.65  l/m6 

High-range  water  reducing  admixture  “Rheobuild  716” 

1.6  //m3 

WES5000  was  first  formulated  at  the  U.S.  Army  Corps  of  Engineer  Research 
and  Development  Center  Waterways  Experiment  Station  (WES)  as  a standard 
test  concrete  for  phenomenology  studies.  As  such,  a great  deal  of  data  exist  for 
validation  of  the  model.  The  WES5000  used  in  this  study  was  poured  by  WES 
on  6 May  1998.  The  material  was  cast  into  two  barrels  and  allowed  to  cure  for 
approximately  30  days  prior  to  coring.  Two  different  diameter  cores  were  taken. 
The  cores  for  triaxial  quasi-static  testing  were  nominally  50mm  in  diameter  and 
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were  cut  into  lengths  of  approximately  110mm.  The  cores  for  Split  Hopkinson 
Pressure  Bar  (SHPB)  testing  were  nominally  75mm  in  diameter  and  were  cut  into 
lengths  of  approximately  76mm.  The  end  faces  on  each  sample  were  then  ground 
to  assure  parallelism.  Both  unconfined  compression  and  splitting  tension  tests  were 
conducted  at  28  days.  Results  of  these  experiments  are  presented  in  Table  2 2. 


Table  2-2:  Properties  at  28  days  for  WES5000  concrete. 


Test  # 

Unconfined  compression 

Splitting  tension 

1 

45.85  MPa  (6650  psi) 

3.31  MPa  (480  psi) 

2 

45.85  MPa  (6650  psi) 

3.07  MPa  (445  psi) 

3 

45.16  MPa  (6550  psi) 

3.00  MPa  (435  psi) 

Avg. 

45.62  MPa  (6617  psi) 

3.13  MPa  (453  psi) 

2.3  Experimental  Design 

One  of  the  objectives  of  the  experimental  study  was  to  determine  the  quasi- 
static behavior  of  WES5000  under  a variety  of  loading  conditions.  Of  particular 
interest  was  the  influence  of  confining  pressure  on  strength,  and  the  material’s 
deformation  properties.  The  simplest  characterization  experiment  which  can  be 
conducted  is  a uniaxial,  or  unconfined  compression  test.  For  most  civil  engineering 
applications,  this  is  the  only  experiment  which  is  run.  Figure  2-3  shows  a repre- 
sentation of  a uniaxial  compression  test  in  which  both  axial  and  radial  strains  have 
been  recorded.  The  pre-peak  portion  of  the  axial  stress-strain  curve  is  essentially 
linear  over  the  first  one  third.  The  Young’s  modulus  (or  tangent  modulus)  may  be 
taken  as  the  slope  of  this  line.  An  estimate  of  the  Poisson  ratio  may  also  be  deter- 
mined from  the  ratio  of  the  radial  to  axial  strains  over  this  linear  portion.  These 
moduli  would  only  be  applicable  in  static  civil  engineering  applications  where 
pressure  or  compression  are  not  important.  The  most  frequently  reported  measure- 
ment from  such  a test  is  the  peak  stress,  or  unconfined  compressive  strength,  often 
denoted  as  This  term  is  more  useful  as  a quality  control  measure  than  an  actual 
material  property. 
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Figure  2-3:  Stress-strain  curve  for  concrete  uniaxial  compression  test 

ASTM  C 801-91  [18]  defines  the  American  standard  test  method  for  deter- 
mining the  mechanical  properties  of  concrete  under  triaxial  loads.  It  calls  for  a 
series  of  monotonic  Karman  type  experiments  conducted  under  one  of  three  loading 
histories. 

• Type  I . Hydrostatic  pressure  is  increased  monotonically  to  a predetermined  level, 

then  held  constant  while  the  axial  stress  is  increased  monotonically  until 
failure  occurs. 

• Type  II . Hydrostatic  pressure  is  increased  monotonically  to  a predetermined  level, 

then  the  axial  stress  is  held  constant  while  the  lateral  stress  is  decreased 
monotonically  until  failure  occurs  (lateral  expansion  test). 

• Type  III.  The  ratio  of  axial  to  lateral  stresses  is  held  constant  and  the  stresses  are 

increased  monotonically  until  failure  occurs  (proportional  loading  test). 

Under  the  conditions  of  a Karman  experiment,  the  stress  tensor  is  given  by 
Equation  2.1,  where  er3  is  the  lateral,  or  radial  load  applied  to  the  specimen  via  the 
hydraulic  confining  fluid.  This  radial  load  is  also  simply  referred  to  as  confinement 
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or  cr3  throughout  the  text. 
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Hydrostatic  pressure,  or  mean  normal  stress,  is  defined  as  one-third  the  trace  of  the 
stress  tensor,  or 


P = 


-(crj  + 2cr3) 


(2.2) 


where  p is  used  to  represent  the  hydrostatic  pressure,  or  mean  normal  stress. 

In  all  of  these  tests,  axial  and  radial  strains  are  measured,  as  well  as  applied 
axial  and  radial  loads.  Thus,  information  regarding  the  volumetric  response  of 
the  material  and  the  failure  surface  as  a function  of  confining  pressure  may  be 
gathered.  Tests  of  Type  I are  usually  performed  for  characterization  purposes  and 
provide  all  the  information  necessary  for  the  formulation  of  continuum  level  models. 
Tests  of  Type  II  and  III  are  generally  run  for  validation  purposes.  Information 
regarding  the  elastic  moduli  may  be  gathered  in  the  same  means  as  discussed  for 
unconfmed  compression.  In  practice,  unload-reload  cycles  are  generally  included 
in  the  loading  profiles  for  determination  of  the  elastic  properties.  The  moduli  are 
taken  as  the  slope  of  the  unload-reload  curves.  Figure  2-4  displays  the  results 
of  a hydrostatic  compression  test  on  WES5000  with  one  unload-reload  cycle. 
Significant  hysteresis  is  evident  in  the  unload-reload  cycle.  This  hysteresis  can 
make  determination  of  the  unload-reload  slopes  highly  subjective.  Cristescu  [12,  7] 
has  shown  that  the  introduction  of  a creep  stage  (i.e.  stress  state  held  constant) 
prior  to  the  unload-reload  cycle  eliminates  the  viscous  effects  such  that  unloading 
is  purely  elastic.  In  a creep  experiment,  the  specimen  is  allowed  to  deform  under 
a state  of  constant  stress.  The  evolution  of  the  strain  with  time  is  monitored  until 
the  strain  rate  of  the  material  approaches  zero.  This  is  indicative  that  the  material 


has  reached  a stabilized  state.  Unloading  and  reloading  from  this  stabilized  state 
should  not  be  influenced  by  viscous  properties  and  the  material  response  should 
be  purely  elastic.  In  order  to  evaluate  the  viscosity  of  the  material,  as  well  as 
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Figure  2-4:  Monotonic  hydrostatic  compression  test  exhibiting  hysteretic  behavior 

minimizing  its  effect  on  any  unloading-reloading  cycles,  the  decision  was  made  to 
include  a creep  step  prior  to  any  unloading-reloading.  Furthermore,  by  including 
several  such  cycles  (creep-unload-reload)  in  each  phase  of  the  experiment  the 
quantity  of  information  gathered  could  be  maximized.  Figures  2 5 and  2-6  show 
the  proposed  hydrostatic  and  deviator ic  experiments  respectively.  Figure  2-5  is 
presented  in  terms  of  volume  strain  as  a function  of  mean  normal  stress,  where 
volume  strain  is  defined  as  the  trace  of  the  strain  tensor,  or  under  the  conditions  of 
a Karman  experiment  as: 

sv  — £\  + 2^3  (2-3) 

where  e„,  eq,  and  £3  are  the  volume  strain,  axial  strain  and  radial  strain  respec- 
tively. The  terms  iv  and  E\  in  Figures  2-5  and  2-6  are  the  time  rate  of  change, 
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or  strain  rate  for  the  volume  strain  and  axial  strain  respectively.  The  over-dot 
notation  will  be  utilized  throughout  this  text  to  signify  time  rate  of  change. 


Figure  2-5:  Schematic  of  hydrostatic  creep  experiment 

A series  of  such  experiments  was  conducted  at  the  U.  S.  Army  Corps  of  En- 
gineer Engineering  Research  and  Development  Center,  Waterways  Experiment 
Station  in  Vicksburg  Mississippi  to  fully  characterize  the  material.  The  experi- 
mental matrix  is  given  in  Table  2-3.  In  this  table  under  test  type,  TXC  refers  to  a 
triaxial  compression  test,  while  UC  refers  to  uniaxial  compression.  Creep  implies 
creep  steps  were  included.  One  experiment,  9902B016,  was  run  with  relaxation 
steps.  The  presence  of  asterisks  in  the  data  columns  implies  that  no  data  of  that 
type  was  taken  in  the  experiment.  This  was  usually  due  to  failure  of  the  membrane 
separating  the  material  specimen  from  the  hydraulic  fluid. 
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Figure  2-6:  Schematic  of  deviatoric  shear  creep  experiment 

Prior  to  each  test,  the  WES5000  specimen  was  carefully  weighed  and  mea- 
sured. Post  test,  the  specimens  were  dried  in  an  oven  to  remove  all  moisture.  This 
allowed  calculation  of  the  specimen  as  tested  density  as  well  as  the  dry  density, 
and  hence  the  as-tested  moisture  content  of  the  material  could  be  evaluated.  Axial 
and  radial  P-wave  and  S-wave  velocity  measurements  were  also  taken.  The  results 
are  given  in  Table  2 3.  WES5000  is  inherently  rather  porous  and  prone  to  having 
voids,  making  it  a difficult  material  with  which  to  work.  This  porous  nature  makes 
it  impractical  to  directly  instrument  the  specimens  with  strain  gauges  except  at 
very  low  confining  pressures.  At  higher  pressures,  the  thin  gauge  substrate  would 
be  prone  to  “popping”  into  the  tiny  voids  or  pores  over  which  it  had  been  applied, 
thus  rendering  it  non-functional.  Rather,  Linear  Variable  Differential  Transformers 
(LVDTs)  were  utilized  to  measure  the  axial  and  radial  deformation  of  the  material. 
The  LVDT  is  an  electromagnetic  device  that  produces  an  electrical  voltage  propor- 
tional to  the  displacement  of  a movable  magnetic  core.  These  LVDTs  were  housed 
in  a “bird  cage”  apparatus  surrounding  the  specimen.  The  radial  LVDTs  were 
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connected  via  two  metal  studs  glued  to  the  specimen  midpoint  at  diametrically 
opposite  points.  Since  pressure  is  transmitted  to  the  specimen  via  a confinement 
fluid,  precautions  must  be  taken  to  prevent  fluid  infiltration  into  the  specimen. 

The  specimen  is  generally  encased  in  a latex  membrane.  Care  must  be  taken  in 
specimen  preparation  to  open  and  patch  all  surface  voids  in  order  to  prevent  the 
membrane  from  being  forced  down  into  these  voids  and  punctured.  A gypsum 
paste  was  used  to  patch  all  voids.  Subsequent  to  patching  of  the  voids,  the  speci- 
men was  encased  in  the  latex  membrane,  and  this  membrane  was  then  surrounded 
by  a layer  of  plumbers  putty.  The  putty  served  the  function  of  self  sealing  any 
potential  punctures  of  the  inner  membrane.  A second  latex  membrane  was  then 
applied  over  the  putty  and  the  membrane  coated  with  a sealant  to  prevent  its 
degradation  by  the  petroleum  based  hydraulic  fluid.  Figure  2 7 displays  a sample 
prepared  for  testing  with  radial  LVDTs  attached.  The  entire  “bird  cage”  is  then 
lowered  into  the  confinement  cell  for  testing.  WES’s  Secodyne  500K  pound  load 
frame  (see  Figure  2-8)  coupled  with  a 600  MPa  capacity  Karman  cell  was  utilized 
for  all  the  experiments.  The  loading  rate/profile  was  computer  programmable, 
allowing  for  precise  control  of  the  experiments  under  load,  displacement  or  strain 
control. 

2.4  Unconfined  Compression 

The  first  experiment  conducted  (9902A02)  was  a monotonic  unconfined  com- 
pression test.  While  not  generally  revealing  from  a model  development  standpoint, 
the  unconfined  compressive  strength,  or  f'c,  is  an  often  quoted  figure  of  merit  when 
comparing  various  concrete  mixes  or  individual  batches.  Figure  2-9  shows  the 
results  of  the  experiment.  In  this  case,  the  specimen  was  instrumented  with  both 
axial  and  radial  strain  gauges  as  well  as  axial  and  radial  LVDTs.  Note  the  very 
close  correlation  apparent  between  the  LVDT  and  radial  strain  gauge  signals.  For 
the  axial  strain,  there  appears  to  be  an  initial  difference  in  the  signals  with  the 
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Figure  2-7:  Sample  prepared  for  testing 
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Figure  2-8:  Waterways  Experiment  Station  Secodyne  500K  pound  load  frame 
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LVDTs  indicating  a greater  amount  of  axial  deformation.  This  initial  discrep- 
ancy is  most  likely  due  to  localized  effects  at  the  platen-specimen  interface  at  the 
onset  of  the  experiment.  The  strain  gauges  would  not  record  these  effects  since 
they  could  measure  deformations  only  in  the  vicinity  of  the  point  where  they  are 
attached.  The  LVDTs  on  the  other  hand  record  the  deformation  of  the  entire  speci- 
men. However,  after  this  initial  stage,  which  shall  be  referred  to  as  “seating” , strain 
gauge  and  LVDT  signals  are  in  excellent  agreement.  The  compressive-dilatant 
nature  of  the  material  can  be  determined  from  examination  of  the  volume  strain 
curve.  A positive  volume  strain  slope  is  indicative  of  material  compressibility, 
while  a negative  slope  indicates  material  dilatancy.  Therefore,  from  the  volume 
stress-strain  curve,  it  may  be  concluded  that  the  material  is  initially  compressible, 
but  becomes  dilatant  as  the  stress  increases.  The  transition  point  from  compress- 
ibility to  dilatancy  is  referred  to  as  the  compressibility-dilatancy  boundary.  The 
unconfined  compressive  strength  at  121  days  curing  time  was  found  to  be  60  MPa 
(8685  psi). 


Figure  2-9:  Unconfined  compression  test 
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2.5  300  MPa  Experiment 

The  initial  test  matrix  developed,  called  for  a series  of  hydrostatic  compres- 
sion phase  followed  by  deviatoric  phase  experiments  with  several  loading-creep- 
unloading-reloading cycles.  The  loading  rates  were  to  be  held  such  that  the  axial 
strain  rate  would  be  on  the  order  of  1CT6.  Creep  cycle  duration  was  to  be  on  the 
order  of  30  minutes  to  allow  the  material  to  achieve  stabilization.  Tests  were  run 
under  these  conditions  at  150,  300  and  450  MPa  confinement.  Results  of  the  300 
MPa  experiment  (9902B01)  are  presented  in  Figures  2-10  through  2-15.  The 
convention  observed  in  this  dissertation  is  that  compressive  stresses  and  strains 
are  positive.  Measurements  were  recorded  at  one  second  intervals.  Figure  2 10 
shows  the  hydrostatic  compression  phase  of  the  experiment.  At  the  onset  of  load- 
ing, there  appears  to  be  a sudden  negative  jump  in  the  radial  strain.  This  is  not 
believed  to  be  indicative  of  any  material  behavior,  rather  it  is  likely  associated 
with  specimen  “seating” . The  axial  and  radial  strains  correlate  well,  hence  the 
material  can  be  considered  isotropic.  If  a correction  were  made  for  the  initial  offset 
on  the  radial  LVDT,  the  correlation  between  the  two  signals  would  be  even  greater. 
Upon  examination  however,  the  response  of  the  unload-reload  cycles  for  the  radial 
signal  appears  far  too  stiff.  The  degree  of  hysteresis  in  the  signal  also  appears  to  be 
significant.  This  is  especially  peculiar  since  there  is  virtually  no  hysteresis  apparent 
in  the  axial  signal.  Figure  2-11  displays  only  the  first  two  cycles  of  the  hydrostatic 
loading.  Here  the  initial  radial  offset  is  more  visible.  It  is  also  more  clearly  visible 
that  there  is  virtually  no  hysteresis  present  in  the  axial  unload-reload  cycles.  This 
is  indicative  of  the  fact  that  the  material  has  had  time  to  reach  stabilization  via 
creep  (i.e.  it  is  on  the  stabilization  boundary)  and  that  it  is  responding  elastically 
upon  unloading.  This  is  not  observed  in  the  radial  signal,  further  bringing  into 
question  the  validity  of  the  radial  unload-reload  data.  Figures  2-12  through  2-15 
present  the  deviatoric  response  of  the  material.  The  quality  of  the  axial  signal 
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appears  quite  good,  with  little  apparent  hysteresis  in  the  unload-reload  cycles.  A 
peak  principal  stress  difference  (PSD)  of  504  MPa  was  achieved  in  the  experiment. 
Figures  2-13  and  2-14  show  the  radial  deformeter  data.  Inspection  of  the  first 
creep  cycle  indicates  that  the  material  is  creeping  in  the  positive  direction  (radial 
contraction).  According  to  Cristescu  [8]  this  behavior  has  been  observed  in  creep 
testing  of  highly  porous  rocks.  At  the  second  creep  stage,  the  material  is  exhibiting 
oscillatory  behavior  between  expansion  and  contraction.  On  subsequent  cycles, 
the  material  exhibits  only  radial  expansion  in  the  creep  cycles.  There  is  still  some 
degree  of  hysteresis  apparent  in  the  radial  unload-reload  cycles,  though  far  less 
than  observed  in  the  hydro  phase  of  the  experiment.  The  first  three  cycles  alone 
are  displayed  in  Figure  2-14.  Here  the  creep  behavior  can  more  clearly  be  seen. 
There  is  a good  deal  of  noise  in  the  data  surrounding  the  second  creep  cycle.  This 
may  have  been  due  to  localized  aggregate  cracking  or  void  collapse  influencing  the 
radial  signals.  There  was  also  a certain  amount  of  oscillation  recorded  at  the  onset 
of  the  deviatoric  test.  This  is  not  believed  to  be  indicative  of  any  true  material  be- 
havior. Overall,  the  confidence  in  the  radial  LVDT  data  is  considered  low.  Figure 
2-15  depicts  the  axial  radial  and  volumetric  curves  together.  It  can  be  seen  from 
the  volume  curve  that  the  material  exhibits  purely  compressible  behavior  up  to  the 
onset  of  failure.  The  complete  experiment  was  conducted  over  a period  of  6 hours 


and  38  minutes. 


Mean  Normal  Stress  (MPa)  tg  Mean  Normal  Stress  (MPa) 
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ire  2-10:  Hydrostatic  compression  phase,  300  MPa  Karman  experiment 


Figure  2-11:  Hydrostatic  compression  phase  (low  pressure  range),  300  MPa 
Karman  experiment 
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Axial  Strain  (mm/mm) 


Figure  2—12:  Deviatoric  phase  of  300  MPa  Karman  experiment,  axial  strain 


Figure  2-13:  Deviatoric  phase  of  300  MPa  Karman  experiment,  radial  strain 
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Radial  Strain  (mm/mm) 


Figure  2-14:  Deviatoric  phase  of  300  MPa  Karman  experiment,  radial  strain  (low 
PSD  range) 


Figure  2-15:  Deviatoric  phase  of  300  MPa  Karman  experiment,  all  strains 
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2.6  450  MPa  Experiment 

Results  of  the  tests  run  at  450  MPa  radial  confinement  (9902B02)  are  pre- 
sented in  Figures  2-16  through  2-21.  Figure  2-16  shows  the  full  hydrostatic 
compression  phase  of  the  experiment.  At  first  glance,  the  behavior  of  the  material 
does  not  appear  to  be  isotropic  in  that  the  axial  and  radial  signals  do  not  appear 
to  correlate  well.  However,  the  difference  is  actually  rather  small  (on  the  order  of 
a few  tenths  of  a percent  strain)  and  may  well  be  within  the  accuracy  limits  of  the 
LVDTs.  The  overall  behavior  of  the  radial  signals  is  much  improved  over  those  in 
the  300  MPa  creep  experiment,  with  little  apparent  hysteresis  in  the  unload- reload 
cycles.  Figure  2-17  shows  a detailed  view  of  the  first  2 cycles.  Some  evidence  of 
the  axial  seating  phenomena  is  visible,  but  it  is  not  pronounced.  Figures  2-18 
through  2-21  display  the  deviatoric  phase  of  the  experiment.  There  was  again  an 
apparent  anomaly  in  the  data  at  the  onset  of  the  deviatoric  phase,  much  like  that 
discussed  for  9903B01.  However,  as  presented  here,  the  data  has  been  adjusted 
to  remove  the  initial  offsets.  The  axial  LVDT  results,  as  depicted  in  Figure  2-18, 
appear  quite  good  with  a peak  principal  stress  difference  of  671  MPa  achieved. 
Figure  2-19  presents  the  complete  radial  signal.  Once  again,  there  is  some  degree 
of  hysteresis  apparent  in  the  unload-reload  cycles  which  was  not  seen  in  the  axial 
signals.  Figure  2-20  displays  in  detail  the  first  two  cycles.  The  initial  oscillatory 
behavior  is  again  observed  at  the  first  creep  cycle,  however  the  creep  direction 
is  negative  for  the  second  and  all  subsequent  cycles.  Finally,  Figure  2-21  shows 
the  axial,  radial,  and  computed  volumetric  curves  together.  Again  it  is  seen  that 
the  material  is  compressible  until  very  near  the  onset  of  failure.  The  complete 
experiment  was  conducted  over  8 hours  and  12  minutes. 
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Strain  (mm/mm) 


Figure  2-16:  Hydrostatic  compression  phase,  450  MPa  Karman  experiment 


Strain  (mm/mm) 


Figure  2-17:  Hydrostatic  compression  phase  (low  pressure  range),  450  MPa 
Karman  experiment 
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Axial  Strain  (mm/mm) 


Figure  2-18:  Deviatoric  phase  of  450  MPa  Karman  experiment,  axial  strain 


Radial  Strain  (mm/mm) 


Figure  2-19:  Deviatoric  phase  of  450  MPa  Karman  experiment,  radial  strain 
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Figure  2-20:  Deviatoric  phase  of  450  MPa  Karman  experiment,  radial  strain  (low 
PSD  range) 


Strain  (mm/mm) 


Figure  2-21:  Deviatoric  phase  of  450  MPa  Karman  experiment,  all  strains 
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2.7  150  MPa  Experiment 

Results  of  the  150  MPa  experiment  (9902B08)  are  presented  in  Figures 
2-22  through  2-27.  Figure  2-22  displays  the  hydrostatic  compression  phase  of 
the  experiment.  Results  seem  to  indicate  that  the  material  is  not  behaving  very 
isotropically.  Even  taking  into  account  the  initial  platen-specimen  interface  seating 
phenomena,  the  radial  response  appears  much  stiffer.  This  can  be  seen  clearly  in 
the  slopes  of  the  corresponding  unload-reload  cycles.  Figure  2-23  presents  the  first 
two  cycles  of  the  hydrostatic  compression  experiment.  Here  the  platen-specimen 
interface  phenomena  discussed  earlier  for  the  unconfined  experiment  can  be  clearly 
seen.  Figures  2-24  through  2-27  present  the  deviatoric  response  of  the  material. 
The  axial  strain  results  appear  quite  good,  with  a peak  principal  stress  difference 
of  approximately  390  MPa  achieved.  Once  again,  the  radial  signals  indicate  a 
certain  amount  of  hysteresis.  Radial  contraction  is  again  seen  in  the  first  creep 
cycle  with  the  onset  of  the  oscillatory  phenomena  visible  in  the  second  cycle.  This 
can  be  more  clearly  seen  in  Figure  2-  26.  Figure  2 27  displays  all  the  strain  signals 
together.  Again  the  transition  to  dilatancy  is  not  apparent  in  the  volume  strain 
signal  until  the  onset  of  failure. 
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Figure  2-22:  Hydrostatic  compression  phase,  150  MPa  Karman  experiment 


Figure  2-23:  Hydrostatic  compression  phase,  150  MPa  Karman  experiment,  first  2 
cycles 
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Axial  Strain  (mm/mm) 


2-24:  Deviatoric  phase  of  150  MPa  Karman  experiment,  axial  strain 


Figure  2-25:  Deviatoric  phase  of  150  MPa  Karman  experiment,  radial  strain 
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Figure  2-26:  Deviatoric  phase  of  150  MPa  Karman  experiment,  radial  strain  (low 
PSD  range) 


Figure  2-27:  Deviatoric  phase  of  150  MPa  Karman  experiment,  all  strains 
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2.8  500  MPa  Hydrostatic  Compression  Experiment 

Following  completion  of  the  150  MPa  experiment,  a decision  was  made  to 
re-examine  the  experimental  procedure.  Consistent  anomalies  had  been  apparent 
in  the  radial  strain  measurements  for  the  experiments  conducted  so  far.  These 
were  believed  to  be  at  least  partly  due  to  calibration  of  the  LVDTs.  For  a full 
experiment,  the  expected  level  of  radial  deformation  in  the  deviatoric  phase  can  be 
expected  to  be  on  the  order  of  10%  (see  Figure  2-25),  while  the  total  deformation 
in  the  compression  phase  is  generally  on  the  order  of  1%.  Therefore,  LVDTs  ranged 
for  the  total  expected  level  of  deformation  in  the  compression-deviatoric  shear 
experiment  may  not  be  sensitive  enough  to  record  the  very  small  strains  in  the 
hydrostatic  phase.  The  duration  of  the  experiments  was  also  proving  to  be  taxing 
on  the  equipment  and  technicians,  this  in  turn  was  driving  up  the  cost  of  the 
experimentation.  In  order  to  resolve  the  issue  of  instrumentation  sensitivity,  it  was 
decided  to  conduct  one  additional  purely  hydrostatic  compression  experiment  with 
multiple  creep-unload-reload  cycles.  Since  only  the  hydrostatic  portion  would  be 
conducted,  the  LVDTs  could  be  more  appropriately  calibrated.  All  subsequent 
experiments  would  be  run  with  a monotonic  hydrostatic  compression  phase, 
followed  by  a full  deviatoric  test  (with  creep-unload-reload  cycles).  Since  at  the  end 
of  the  creep  cycle  the  material  response  is  theoretically  purely  elastic,  the  unloading 
could  be  conducted  at  as  high  a rate  which  could  be  easily  controlled.  These  last 
two  modifications  were  made  in  an  effort  to  shorten  the  length  of  the  experiments. 

Results  of  the  500  MPa  hydrostatic  compression  experiment  (9902B012)  are 
presented  in  Figure  2-28.  A total  of  7 creep-unload-reload  cycles  were  conducted, 
with  the  creep  pressure  held  for  a duration  of  15  minutes.  The  material  appears  to 
behave  quite  isotropically  up  to  the  fifth  creep-unload-reload  cycle.  At  this  point 
the  response  of  the  radial  gauges  indicates  much  stiffer  response  than  is  indicated 
by  the  axial  gauges.  Overall,  the  quality  of  the  data  appears  excellent  for  both 
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axial  and  radial  measurements.  Figure  2-29  displays  the  first  two  creep-unload- 
reload  cycles  in  detail.  Once  again,  there  is  good  correlation  between  the  two 
signals.  The  absence  of  a discernable  creep  plateau  at  50  MPa  is  indicative  that 
the  material  is  responding  largely  elastically  at  this  pressure  range.  Finally,  Figure 
2-30  displays  a plot  of  volume  strain  rate  as  a function  of  time  for  the  500  MPa 
hydrostatic  compression  experiment.  The  plot  indicates  that  for  the  four  cycles, 
the  volume  strain  rate  has  dropped  to  very  nearly  zero  by  the  end  of  the  creep 
phase.  As  predicted,  there  is  virtually  no  hysteresis  apparent  in  the  unload-reload 
phase  following  these  initial  cycles.  In  the  final  three  cycles,  the  volume  strain  rate 
does  not  reach  zero  prior  to  the  unload-reload  phase  (there  is  no  reload  following 
the  seventh  creep  cycle),  and  there  is  indeed  a small  degree  of  hysteresis  apparent 
following  the  sixth  cycle.  The  overall  duration  of  the  experiment  was  3 hours  and  3 
minutes. 


Strain  (mm/mm) 


Figure  2-28:  Hydrostatic  compression  Karman  experiment  to  500  MPa 


Volume  Strain  Rate 
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Figure  2-29:  Hydrostatic  compression  Karrnan  experiment  to  500  MPa,  first  2 
cycles 


Figure  2-30:  Volume  strain  rate  vs.  time  for  500  MPa  hydrostatic  compression 
experiment 
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2.9  50  MPa  Experiment 

Following  the  500  MPa  hydrostatic  compression  test,  experiments  were  run  at 
50  MPa,  200  MPa,  300  MPa,  100  MPa  and  375  MPa  radial  confinement.  Results  of 
the  50  MPa  experiment  (9902B013A)  are  presented  in  Figures  2-31  through  2-35. 
Figure  2-31  displays  the  results  of  the  monotonic  hydrostatic  compression  phase 
of  the  experiment.  Close  examination  of  the  data  shows  a small  initial  negative 
offset  to  the  axial  data,  indicative  of  specimen  seating.  There  does  not  appear 
to  be  good  correlation  between  the  axial  and  radial  signal  for  this  phase  of  the 
experiment. One  should  consider  however  that  the  magnitudes  of  the  measured 
strain  is  quite  small  with  respect  to  the  calibration  of  the  instrumentation.  Figure 
2-32  displays  the  axial  strain  signal  as  a function  of  principal  stress  difference  for 
the  deviator ic  shear  portion  of  the  experiment.  The  quality  of  the  data  appears 
quite  good  with  no  hysteresis  evident  in  the  unload- reload  cycles.  There  is  a good 
deal  of  “jitter”  in  the  signal  for  the  final  unload-reload  cycle  and  the  subsequent 
load  to  failure.  This  is  perhaps  indicative  of  material  damage.  Figure  2-33  shows 
the  radial  strain  signal  for  the  deviatoric  shear  phase.  Once  again,  there  is  evidence 
of  radial  expansion  in  the  first  creep  cycle.  The  overall  quality  of  the  data  appears 
good  with  no  significant  hysteresis  in  the  unload-reload  cycles.  The  first  two  creep- 
unload-reload  cycles  are  detailed  in  Figure  2-34.  Here  the  radial  expansion  in  the 
first  cycle  can  be  clearly  seen,  as  well  as  oscillatory  behavior  in  the  second  cycle. 
Close  inspection  also  indicates  that  there  is  some  lag  in  the  response  of  the  radial 
signal  at  the  onset  of  the  deviatoric  loading.  The  initial  lack  of  appreciable  radial 
strain  at  the  onset  of  the  deviatoric  phase,  and  following  the  reload  portion  of  the 
first  cycle  is  evidence  of  this  lag  in  response.  Finally,  Figure  2-35  displays  all  the 
signals  together.  The  transition  from  compressibility  to  dilatancy  is  apparent  in  the 
volumetric  strain  signal  only  near  failure.  The  material  exhibits  strain  hardening 
behavior  up  to  failure,  indicative  of  a brittle  material.  The  sample  achieved  a peak 
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principal  stress  difference  of  221  MPa  prior  to  failure.  The  overall  length  of  the 
experiment  was  2 hours  and  51  minutes. 


Strain  (mm/mm) 

Figure  2—31 : Hydrostatic  compression  phase,  50  MPa  Karman  experiment 


Principal  Stress  Difference  (MPa)  5 Principal  Stress  Difference  (MPa) 
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2-32:  Deviatoric  phase  of  50  MPa  Karman  experiment,  axial  strain 


Figure  2-33:  Deviatoric  phase  of  50  MPa  Karman  experiment,  radial  strain 
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Figure  2-34:  Deviatoric  phase  of  50  MPa  Karman  experiment,  radial  strain  first  2 
cycles 


Strain  (mm/mm) 


Figure  2-35:  Deviatoric  phase  of  50  MPa  Karman  experiment,  all  strains 
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2.10  200  MPa  Experiment 

The  next  series  of  figures  (2-36  through  2-40)  displays  the  results  of  the  200 
MPa  radial  confinement  experiment  (9902B014).  Figure  2 36  shows  the  axial 
and  radial  strain  signals  for  the  monotonic  hydrostatic  compression  phase  of  the 
experiment.  Excellent  correlation  is  indicated  between  the  two  signals  indicating 
a high  degree  of  material  isotropy.  The  next  figure  presents  the  axial  strain  signal 
as  a function  of  principal  stress  difference  for  the  deviatoric  shear  phase  of  the 
experiment.  Again  the  quality  of  the  data  appears  very  good  with  no  hysteresis 
apparent  in  the  unload-reload  cycles.  The  radial  response  for  the  deviatoric  phase 
of  the  experiment  is  presented  in  Figure  2-38.  Again,  there  are  indications  that  the 
gauges  are  not  initially  responding  to  the  applied  deviatoric  load.  The  response  lag 
is  far  more  pronounced  in  this  experiment  than  it  was  in  the  50  MPa  experiment. 
The  first  creep  cycle  again  displays  radial  expansion.  There  does  not  appear  to 
be  any  appreciable  level  of  hysteresis  in  the  final  two  creep-unload-reload  cycles. 

An  expanded  view  of  the  first  two  cycles  is  shown  in  Figure  2-39.  Here,  the  lag  in 
response  is  clear  not  only  at  the  onset  of  deviatoric  loading,  but  also  following  the 
first  unload  cycle.  It  is  most  likely  that  this  lag  is  due  to  a mechanical  problem 
in  the  LVDT  cage  apparatus.  The  slope  of  the  reload  side,  once  the  gage  begins 
responding,  is  very  close  to  that  from  the  unload  side.  The  second  cycle  once  again 
displays  oscillatory  behavior  in  the  creep  cycle  and  there  is  again  indication  of  the 
LVDT  not  responding  for  the  reload  portion  of  the  cycle,  though  to  much  less  a 
degree.  The  volumetric,  axial  and  radial  signals  are  presented  together  in  Figure 
2-40.  Once  again,  the  transition  from  compressibility  to  dilatancy  is  only  apparent 
at  the  onset  of  failure.  The  material  does  not  show  signs  of  strain  hardening  in  this 
experiment,  rather,  near  failure,  the  stress  level  is  constant  for  increasing  strain. 
This  is  indicative  of  a transition  from  brittle  to  ductile  behavior.  A peak  principal 
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stress  difference  of  457  MPa  was  achieved  in  the  experiment,  which  was  conducted 
over  a period  of  3 hours  and  44  minutes. 


Strain  (mm/mm) 


Figure  2-36:  Hydrostatic  compression  phase,  200  MPa  Karman  experiment 


Principal  Stress  Difference  (MPa)  ^ Principal  Stress  Difference  (MPa) 
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Axial  Strain  (mm/mm) 

2-37:  Deviatoric  phase  of  200  MPa  Karman  experiment,  axial  strain 


Figure  2-  38:  Deviatoric  phase  of  200  MPa  Karman  experiment,  radial  strain 
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Figure  2-39:  Deviatoric  phase  of  200  MPa  Karman  experiment,  radial  strain  first  2 
cycles 


Figure  2-40:  Deviatoric  phase  of  200  MPa  Karman  experiment,  all  strains 
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2.11  300  MPa  Experiment  (Second) 

The  next  experiment  was  conducted  at  a radial  confinement  of  300  MPa 
(9902B015).  This  second  experiment  under  the  300  MPa  confinement  conditions 
was  intended  to  serve  two  purposes.  First,  it  was  desirable  to  conduct  a second 
experiment  under  nearly  identical  conditions  due  to  the  obvious  problems  with  the 
radial  deforometers  in  the  first  experiment.  Second,  it  would  give  some  qualitative 
sense  regarding  the  reproducibility  of  the  experiments.  Figure  2-41  displays  the 
results  of  the  monotonic  hydrostatic  compression  phase  of  the  experiment.  There 
is  generally  very  good  agreement  between  the  two  signals.  Axial  strain  results  for 
the  deviatoric  shear  phase  of  the  experiment  are  presented  in  Figure  2 -42.  Once 
again,  the  quality  of  the  data  is  excellent  with  no  significant  hysteresis  in  the 
unload-reload  cycle.  Following  the  final  creep  cycle,  the  data  exhibits  a significant 
amount  of  noise  which  is  likely  due  to  aggregates  breaking  in  the  specimen.  Figure 
2-43  displays  the  radial  strain  signal  as  a function  of  principal  stress  difference. 
The  familiar  radial  expansion  and  subsequent  oscillating  behavior  in  the  first  two 
creep  cycles  is  again  evident.  Following  the  final  creep-unload-reload  cycle,  there 
is  a significant  “glitch”  in  the  data  indicated  by  a sudden  radial  expansion.  This 
point  corresponds  to  the  onset  of  noise  in  the  axial  signal  and  is  likely  due  to  a 
localized  aggregate  fracture  or  void  collapse  in  the  vicinity  of  the  radial  gauge 
mounts.  Figure  2-44  presents  a detailed  view  of  the  radial  strain  for  the  first 
two  creep-unload-reload  cycles.  There  is  no  evidence  here  of  the  response  delay 
observed  in  the  previous  two  experiments.  All  the  signals  are  presented  together 
in  Figure  2-45.  Again,  the  transition  from  compressibility  to  dilatancy  is  only 
very  near  material  failure.  The  entire  experiment  was  conducted  over  a period  of 
4 hours  and  15  minutes,  a significant  reduction  from  the  6 hours  and  38  minutes 
in  the  first  300  MPa  experiment.  A peak  principal  stress  difference  of  570  MPa 
was  achieved  in  the  experiment,  somewhat  higher  than  the  504  MPa  measure  in 


the  first  300  MPa  experiment.  This  difference  is  no  doubt  due  partly  to  natural 
experimental  scatter,  but  also  to  the  additional  work  performed  on  the  specimen 
the  hydrostatic  compression  creep  cycles  of  the  first  experiment. 


Strain  (mm/mm) 


Figure  2-41:  Hydrostatic  compression  phase,  300  MPa  Karman  experiment 


Principal  Stress  Difference  (MPa)  % Principal  Stress  Difference  (MPa) 
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2—42:  Deviatoric  phase  of  300  MPa  Karman  experiment,  axial  strain 


Figure  2-43:  Deviatoric  phase  of  300  MPa  Karman  experiment,  radial  strain 
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Figure  2-44:  Deviatoric  phase  of  300  MPa  Karman  experiment,  radial  strain  first  2 
cycles 


Figure  2-45:  Deviatoric  phase  of  300  MPa  Karman  experiment,  all  strains 
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2.12  100  MPa  Experiment 

The  next  experiment  conducted  was  at  100  MPa  radial  confinement  (9902B016). 
This  experiment  differed  from  the  previous  ones  in  that  relaxation  cycles  were  in- 
cluded in  lieu  of  creep  cycles.  In  a relaxation  cycle,  the  applied  stress  is  controlled 
in  such  a way  that  the  strain  remains  constant.  For  this  experiment,  only  the  axial 
stress  was  feedback  controlled,  the  radial  confining  pressure  remained  constant. 

This  experiment  was  included  in  the  matrix  to  provide  an  alternate  stress  path  for 
later  constitutive  model  validation.  Results  of  the  monotonic  hydrostatic  compres- 
sion phase  of  the  experiment  are  depicted  in  Figure  2 46.  Correlation  between  the 
axial  and  radial  strain  signals  are  again  quite  good.  The  axial  strain  as  a function 
of  principal  stress  difference  for  the  deviatoric  shear  phase  of  the  experiment  is 
presented  in  Figure  2-47.  The  quality  of  the  data  appears  quite  good  with  no 
hysteresis  apparent  in  the  unload-reload  cycles.  Radial  strain  results  are  depicted 
in  Figure  2-48.  Since  the  confining  pressure  was  not  adjusted  to  control  the  radial 
strain,  data  was  not  constant  throughout  the  relaxation  phase.  The  response  ap- 
pears to  be  initially  invariant,  followed  by  radial  expansion  up  to  the  unload  cycle 
where  the  slope  increases  significantly.  During  the  reload  cycle,  the  LVDT  response 
lag  seems  to  have  reappeared.  This  response  lag  can  more  clearly  be  seen  in  Figure 
2-49,  which  depicts  the  first  two  relaxation- unload- reload  cycles  in  detail.  Finally, 
all  the  signals  are  presented  together  in  Figure  2-50.  In  this  experiment,  the  mate- 
rial was  not  tested  to  failure  and  as  the  volumetric  strain  curve  indicates,  remained 
purely  compressive  throughout  the  test.  The  overall  duration  of  the  experiment 
was  2 hours  and  45  minutes. 


Principal  Stress  Difference  (MPa)  ^ Mean  Normal  Stress  (MPa) 
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Strain  (mm/mm) 

2-46:  Hydrostatic  compression  phase,  100  MPa  Karrnan  experiment 


Axial  Strain  (mm/mm) 


Figure  2-47:  Deviatoric  phase  of  100  MPa  Karrnan  experiment,  axial  strain 
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Figure  2-48:  Deviatoric  phase  of  100  MPa  Karman  experiment,  radial  strain 


Figure  2-49:  Deviatoric  phase  of  100  MPa  Karman  experiment,  radial  strain  first  2 
cycles 
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Figure  2-50:  Deviatoric  phase  of  100  MPa  Karrnan  experiment,  all  strains 

2.13  375  MPa  Experiment 

The  final  experiment  in  our  matrix  was  conducted  at  375  MPa  radial  confine- 
ment (9902B017).  Results  of  the  monotonic  hydrostatic  compression  phase  are 
given  in  Figure  2-51.  Relatively  good  correlation  is  seen  in  the  axial  and  radial 
signals.  The  axial  strain  as  a function  of  principal  stress  difference  is  presented  in 
Figure  2-52.  The  quality  of  the  axial  data  is  quite  good  with  little  hysteresis  in 
the  unload-reload  cycle.  The  signal  again  becomes  noisy  towards  the  end  of  the 
deviatoric  loading  as  the  specimen  approached  failure.  The  radial  strain  for  the 
deviatoric  shear  phase  is  presented  in  Figure  2-53.  Oscillatory  behavior  is  again 
evident  in  the  single  creep  cycle  conducted.  A detailed  view  of  this  creep- unload- 
reload  cycle  is  presented  in  Figure  2-54.  The  response  lag  for  the  radial  LVDT  does 
not  appear  to  be  an  issue  in  this  experiment.  All  the  signals  are  presented  together 
in  Figure  2-55.  Again,  the  transition  from  compressibility  to  dilatancy  only  occurs 
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near  failure.  A peak  principal  stress  difference  of  617  MPa  was  achieved.  The 
overall  duration  of  the  experiment  was  3 hours  and  3 minutes. 


Figure  2-51:  Hydrostatic  compression  phase,  375  MPa  Karman  experiment 
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Axial  Strain  (mm/mm) 


Figure  2-52:  Deviatoric  phase  of  375  MPa  Karman  experiment,  axial  strain 


Figure  2-53:  Deviatoric  phase  of  375  MPa  Karman  experiment,  radial  strain 
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Figure  2-54:  Deviatoric  phase  of  375  MPa  Karman  experiment,  radial  strain  first  2 
cycles 


Figure  2-55:  Deviatoric  phase  of  375  MPa  Karman  experiment,  all  strains 
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2.14  Elastic  Parameter  Evaluation 

The  elastic  parameters  are  taken  from  the  slopes  of  the  unload-reload  cycles  in 
the  Karman  experiments.  The  Young’s  modulus,  E,  is  determined  from  the  slope 
of  the  axial  strain  vs.  principal  stress  difference  curve  unload-reload  cycles.  By 
examining  the  recorded  value  of  the  principal  stress  difference,  the  commencement 
of  the  reload  cycle  was  unmistakable.  For  this  reason,  E was  determined  from  this 
portion  of  the  cycle.  In  evaluating  the  slope,  a least  squares  fit  was  made  to  the 
data.  Figure  2-56  presents  a typical  result,  in  this  case,  the  data  was  taken  from 
the  first  unload-reload  cycle  of  the  200  MPa  deviatoric  experiment.  Generally,  the 
r2  value  for  the  least  squares  fit  showed  excellent  correlation  with  the  value  of  0.998 
presented  here  being  quite  typical.  The  term  r2  as  used  here  is  a goodness  of  fit 
statistic  such  that: 


r 


2 


SSE 

SSM 


(2.4) 


where  SSE  is  the  sum  of  the  squared  errors  (residuals)  and  SSM  is  the  sum  of 
squares  about  the  mean.  Values  for  E from  the  different  quasi-static  experiments 
are  presented  as  a function  of  mean  normal  stress  in  Figure  2-57.  Visual  inspection 
of  the  figure  indicates  a very  clear  and  consistent  trend  in  the  data,  with  all  the 
experiments  in  general  agreement. 

Values  for  the  bulk  modulus  were  taken  from  the  slope  of  the  unload-reload 
cycles  of  the  volume  strain  versus  mean  normal  stress  curves,  for  both  the  hydro- 
static and  deviatoric  phases  of  the  Karman  experiments.  For  consistency  with  the 
analysis  of  the  Young’s  modulus,  K values  were  also  estimated  from  the  reloading 
portion  of  each  cycle.  This  proved  to  be  somewhat  problematic  due  to  the  response 
lag  problems  discussed  previously  for  the  radial  deformeters  on  the  reload  cycles. 
Inspection  of  the  data  indicated  that  once  the  LVDTs  actually  started  to  respond, 
the  slope  of  the  data  was  nearly  parallel  to  that  of  the  down  load  cycles.  Therefore, 
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Figure  2-56:  Least  squares  analysis  for  determination  of  Young’s  modulus 

when  necessary,  the  least  squares  fit  to  the  data  was  taken  starting  at  the  point 
when  the  gauges  appeared  to  be  responding  normally.  For  the  evaluation  of  K 
from  the  hydrostatic  compression  data,  since  the  axial  signals  were  felt  to  be  more 
reliable  than  the  radial,  only  the  axial  signals  were  used  in  calculating  the  volume 
strain.  This  was  possible  due  to  the  isotropic  nature  of  the  material  (i.e  ev  — 3£q). 
Results  of  the  data  analysis  for  K are  presented  in  Figure  2-58.  With  the  excep- 
tion of  the  values  determined  from  the  first  test  at  300  MPa  confining  pressure,  the 
experiments  appear  to  be  in  very  good  agreement.  It  was  not  expected  that  the 
results  from  this  experiment  would  correlate  well  due  to  the  problems  observed  in 
the  radial  measurements. 

2.15  General  Comments 

Several  conclusions/observations  can  be  made  at  this  point  regarding  the 
WES5000  quasi-static  Karrnan  experiments.  There  was  generally  very  good  corre- 
lation between  the  axial  and  radial  strain  signals  for  the  hydrostatic  compression 
phase  of  the  experiments.  It  is  therefore  quite  reasonable  to  treat  the  material  as 
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Figure  2-57:  Young’s  modulus  (E)  as  a function  of  mean  normal  stress,  experimen- 
tal data 

isotropic  in  nature.  There  were  recurring  problems  with  the  radial  strain  signals, 
most  likely  tied  to  calibration  issues  and  the  mechanical  linkage  used  to  transmit 
radial  deformation  to  the  LVDTs.  Generally,  there  is  a higher  degree  of  confidence 
in  the  axial  signals  recorded  than  the  radial  signals.  There  was  a consistent  transi- 
tion in  slope  observed  in  the  axial  hydrostatic  signals  in  the  vicinity  off'  110  to  125 
MPa  confinement.  According  to  Cargile  [19],  this  can  be  attributed  to  compression 
cracking  and  breakdown  of  the  mortar.  This  phenomena  was  not  observed  in  the 
two  tests  conducted  at  100  MPa  and  50  MPa.  The  radial  oscillatory  behavior 
observed  in  the  first  two  deviatoric  shear  phase  creep  cycles  was  consistent  in  all 
the  experiments  conducted.  As  was  mentioned  previously,  Cristescu  observed 
this  behavior  previously  in  highly  porous  rocks.  Cristescu  and  Cazacu  [20]  have 
attempted  to  model  this  behavior  in  mortar  using  an  elliptic  yield  function.  In  this 
study,  insufficient  data  was  collected  in  the  lower  principal  stress  difference  regime 
to  adequately  map  out  this  behavior.  The  deformation  observed  in  the  creep  cycles, 
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Figure  2-58:  Bulk  modulus  (K)  as  a function  of  mean  normal  stress,  experimental 
data 


as  well  as  the  effect  of  achieving  stabilization  in  resolving  the  hysteresis  issues  in 
the  unload-reloading  cycles  shows  that  there  are  strong  time  dependencies  in  the 
material.  Figures  2-58  and  2-57  show  that  the  elastic  parameters  are  clearly  pres- 
sure dependent.  Finally,  the  peak  value  of  principal  stress  difference  observed  in 
each  experiment  was  plotted  as  a function  of  radial  confining  pressure.  The  results 
are  presented  in  Figure  2-59.  A very  clear,  remarkably  consistent  trend  is  evident 
in  the  data.  Only  the  data  point  for  the  100  MPa  experiment  appears  somewhat 
low,  however  that  experiment  was  not  actually  run  to  failure. 
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Figure  2-59:  Peak  principal  stress  difference  as  a function  of  radial  confinement 


CHAPTER  3 

DYNAMIC  MATERIAL  CHARACTERIZATION 
In  Chapter  2,  techniques  for  characterization  of  geologic  materials  at  elevated 
pressures  and  quasi-static  strain  rates  were  discussed.  Higher  strain  rate  material 
characterization  (i.e.  lO's-1  and  higher)  is  generally  carried  out  via  bar  impact 
experiments,  or  high-velocity  plate  impact  experiments.  Both  experimental  tech- 
niques will  be  reviewed  in  Section  3.1  of  this  chapter,  as  well  as  their  applicability 
to  geologic  materials.  In  Section  3.2,  issues  regarding  dynamic  testing  will  be 
discussed  and  an  experimental  plan  will  be  detailed  for  characterizing  a mortar  and 
the  material  under  study,  WES5000  under  both  unconfined  and  confined  conditions 
using  a Split  Hopkinson  Pressure  Bar.  The  experimental  set-up  and  test  protocol 
will  be  presented  in  Section  3.3.  The  results  of  the  tests  conducted  along  with 
concluding  remarks  are  given  in  Sections  3.4  and  3.5,  respectively. 

3.1  Dynamic  Characterization  Techniques 
3.1.1  The  Hopkinson  Bar 

Higher  strain  rate  testing,  on  the  order  of  101  to  104  s_1,  is  usually  accom- 
plished via  use  of  bar  impact  devices  (see  Figure  2-1).  Bertram  Hopkinson  was 
among  the  first  to  perform  laboratory  investigations  concerning  the  propagation  of 
elastic  waves  in  bars.  In  1914,  he  reported  on  a device  for  measuring  the  pressure 
produced  by  explosive  and  impact  loading  [21],  His  device  consisted  of  a steel 
cylindrical  bar  suspended  horizontally  by  wires  so  that  it  could  swing  in  a pen- 
dulum fashion.  One  end  of  the  bar  was  referred  to  as  the  “firing  end”,  which  was 
subjected  to  explosive  loading  or  ballistic  impact,  causing  a compressive  elastic 
stress  pulse  to  be  propagated  down  the  bar.  To  the  other  end  of  the  bar,  Hopkin- 
son attached  small  disks  constructed  of  the  same  material  as  the  bar.  Hopkinson 


61 


62 


referred  to  these  disks  as  “time-pieces” , and  they  were  held  on  to  the  end  of  the  bar 
via  either  a wrung  fit  or  magnetic  attraction.  As  an  elastic  stress  pulse  propagated 
down  the  bar  from  the  firing  end,  it  would  transmit  into  the  time-piece  and  reflect 
off  of  the  free  surface  as  a tensile  pulse.  The  time-piece  would  then  separate  from 
the  bar  once  the  reflected  tensile  pulse  reached  the  bar  interface.  Momentum  is 
trapped  in  the  time-piece  corresponding  to  a portion  of  the  original  stress  pulse 
twice  the  length  of  the  time-piece.  By  catching  the  time-piece  in  a ballistic  pendu- 
lum the  captured  momentum  could  be  determined,  while  the  amplitude  swing  of 
the  bar  would  indicate  the  momentum  remaining  in  the  bar.  If  the  time-piece  was 
of  sufficient  length  to  trap  the  momentum  in  it’s  entirety,  the  bar  would  remain 
still  upon  separation.  Thus,  by  varying  the  length  of  the  time-piece,  both  the  pulse 
magnitude  and  duration  could  be  determined.  However, the  technique  did  not  allow 
to  determine  the  shape  of  the  stress  pulse.  This  apparatus  came  to  be  known  as 
the  Hopkinson  pressure  bar.  To  overcome  this  limitation,  Davies,  in  1948,  devised 
an  electrical  capacitance  technique  for  measuring  stress  pulses  in  a Hopkinson  bar 
[22].  The  use  of  electrical  resistance  strain  gauges  to  measure  the  elastic  stress 
pulse  in  bars  did  not  become  common  until  the  late  1950’s  and  early  1960s  [23]. 
Still,  Hopkinson  pressure  bar  is  utilized  today  in  the  measurement  of  explosive 
blast  because  the  environment  is  too  harsh  for  a conventional  gauge  to  survive. 

3.1.2  Compressive  Split  Hopkinson  Pressure  Bar  Testing 

The  Split  Hopkinson  Pressure  Bar  (SHPB)  is  a modified  Hopkinson  pressure 
bar  which  was  proposed  by  Kolsky  [24].  For  this  reason,  in  the  literature  the  SHPB 
is  sometimes  also  referred  to  as  a Kolsky  bar.  The  apparatus  setup  is  depicted  in 
Figure  3-1.  The  SHPB  consists  of  a short  impactor  (striker  bar)  and  two  long  high- 
strength  steel  bars  referred  to  as  the  incident  and  transmitter  bars,  respectively. 

The  specimen  is  sandwiched  between  the  incident  and  transmitter  bar.  The  striker 
bar  impacts  the  incident  bar  and  produces  a nearly  square  elastic  compressive 
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stress  pulse  in  the  incident  bar.  The  pulse  propagates  down  the  incident  bar 
and  at  the  specimen  bar  interface,  a portion  of  the  pulse  is  transmitted  to  the 
specimen,  and  a portion  is  reflected  back  due  to  an  impedance  mismatch.  At  the 
other  end  of  the  specimen,  a portion  of  the  wave  is  likewise  transmitted  to  the 
“transmitter”  bar  and  a portion  reflected  back  into  the  specimen.  Electrical  strain 
gauges  mounted  at  the  mid  point  of  the  incident  and  transmitter  bars  record  the 
passage  and  magnitude  of  the  incident,  reflected,  and  transmitted  stress  pulses.  A 
brief  summary  of  the  SHPB  equations  for  obtaining  applied  load,  strain,  and  strain 
rate  is  given  below. 


Striker  Bar 


Incident  Bar 


Transmitter  Bar 


vr 


strain  gages 


Test  Specimen 
(cylinder) 


strain  gages 


Figure  3-1:  Split  Hopkinson  pressure  bar 


Consider  the  elastic  impact  of  two  circular  cross-section  bars  as  depicted 
in  Figure  3-2.  Let  E,  A,  p,  L,  V,  and  C represent  the  bars  respective  Young’s 
modulus,  cross-sectional  area,  density,  length,  velocity,  and  acoustic  wave  speed. 
Bar-1  impacts  bar-2  at  its  left  end,  thus  V\  > V2  . Assuming  that  the  two  bars 
remain  in  contact  for  a finite  period  of  time  r,  then  the  load  at  the  interface  of 
each  bar  is  the  same,  i.e. 


P\  — o\A\  — cr2A2  — P2 
Conservation  of  momentum  yields: 


(3.1) 


— P\A\C\(V\  — v)  — — p2A2C2(v  — V2) 


(3.2) 
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Figure  3-2:  Elastic  impact  of  two  circular  bars 
where  v therefore  is: 


P\A\C\X\  + P2A2C2V2 
P\A\C\  + P2A2C2 


Hence,  the  stresses  in  the  two  bars  may  be  expressed  as: 


_ . n P2A2C2(V\  — V2) 

(7l~~Pl  lCl  p^c,  + P2A2C2 


(3.4) 


<7  2 — —P2A2C2 


PiA1Cl(V,-V2) 

P\A\C\  + P2A2C2 


(3.5) 


If  the  two  bars  are  of  the  same  material  and  cross  sectional  area,  then: 


P\A\C\  — P2A2C2  — pACo 


(3.6) 
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and 


0\  — Oi  — a. 


The  peak  stress,  a,  then  is: 


a = pC0 


(Vi  ~ ^2) 
2 


(3.7) 


(3.8) 


Let  bar-1  represent  the  striker  and  bar-2  the  incident  in  a SHPB  apparatus. 
Since  the  incident  bar  is  always  initially  at  rest,  V2  = 0,  and  the  striker  bar  velocity 
may  simply  be  denoted  by  V.  Then,  Equation  3.8  becomes 

a = pC„|,  (3.9) 

The  compressive  pulse,  of  magnitude  given  by  Equation  3.9,  propagates  back 
down  the  striker  bar  away  from  the  impact  point.  The  leading  edge  of  the  pulse 
reflects  from  the  free  surface  at  the  back  face  of  the  striker  bar,  L\,  as  a tensile 
pulse.  Upon  it’s  arrival  back  at  the  impact  interface,  the  two  bars  separate,  leaving 
a compressive  pulse  traveling  in  the  incident  bar,  the  length  of  the  pulse  being: 


r = 


2 Lx 
Co" 


(3.10) 


Note  that  the  pulse  length  and  magnitude  is  determined  by  the  striker  bar  length 
and  velocity  respectively  (Equations  3.9  and  3.10). 

Load  and  Strain  Analysis 

Figure  3-3  shows  an  expanded  view  of  a material  specimen  sandwiched 
between  the  incident  and  transmitted  bar  for  a SHPB  compression  test.  Here 
the  subscripts  1 and  2 are  used  to  denote  the  incident  and  transmitted  bar  ends, 
respectively.  Let  £iy  er,  and  et  be  the  incident,  reflected,  and  transmitted  strains 
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Figure  3-3:  Split  Hopkinson  pressure  bar  expanded  view 


respectively,  while  U\  and  u2  denote  the  displacements  of  the  ends  of  the  incident 
and  transmitted  bars,  respectively  . Under  the  assumption  that  wave  propagation 
in  the  bars  is  one  dimensional  and  that  the  behavior  is  elastic,  it  follows  that: 


d2u  2 d2u 

OF  = c°d72 


(3.11) 


where  t = time,  x = the  1-D  Lagrangian  coordinate,  and  Co  is  the  elastic  wave 
speed  in  the  bar  defined  as: 

Co  = (3.12) 

P 

The  classical  d’Alembert  solution  to  this  1-D  wave  equation  is: 


u(x,  t)  = f(x  - C0t ) + F(x  + C0t)  (3.13) 

where  / and  F are  arbitrary  functions  depending  on  the  initial  conditions. 

The  function  / represents  a wave  travelling  in  the  positive  x direction(to  the  right 
in  our  frame  of  reference),  while  F corresponds  to  a wave  traveling  in  the  opposite 
direction.  Thus,  / and  F may  be  thought  of  as  the  incident  and  reflected  waves. 
Equation  3.13  may  then  be  written: 


u(x,  t ) = f(x  — C0t ) + F(x  + Cot)  = Ui  + ur. 


(3.14) 
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Since  the  strain  is  : 


£ 


(hi 

dx 


(3.15) 


Substitution  into  Equation  3.14  yields: 


du(x,  t ) 
dx 


o r\ 

= ~Q^j(X  ~ ^O^)  + = ei  + er 


(3.16) 


where  £i  and  er  represent  the  incident  and  reflected  strains  respectively. 

The  deflection  of  the  incident  and  transmitter  bar  interfaces  are  given  by: 


and 


U\ 


Co£\dt 


CoE^dt 


(3.17) 


(3.18) 


Assuming  stresses,  strains,  and  loads  are  positive  in  compression,  the  Equations 
3.17  and  3.18  may  be  written  in  terms  of  the  incident  (e*),  reflected  (er),  and 
transmitted  (et)strains. 


u\  — Co 


+ er)  dt 


(3.19) 


U2  — Cq 


(3.20) 


The  average  strain  in  the  specimen  is: 


U\  — u2 

Es  = 


l 


(3.21) 
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where  ls  is  the  specimen  length.  Substitution  of  Equations  3.19  and  3.20  into 
Equation  3.21  yields: 

£s  = 7^  [ {£i  + £r-et)dt  (3.22) 

‘s  Vo 

Next,  using  Hooke’s  law,  the  load  on  each  end  of  the  specimen  is  obtained: 

P\  — EA  (si  sr)  (3.23) 

and 

P2  = EAet  (3.24) 

where  E an  A are  the  Young’s  modulus  and  cross-sectional  area  of  the  SHPB, 
respectively.  The  average  force  in  the  specimen  then  becomes: 

EA 

Pav  = ~ X-  (£i  + £r  + £t)  (3.25) 

A general  assumption  is  that  the  specimen  is  in  a state  of  uniform  stress  or  in  other 
words,  the  loads  at  the  ends  of  the  specimen  are  equal.  It  follows  that: 


£%  “b  £r  t-t 


(3.26) 


and 


= 7 - / (£t  - £r  - £r  - £t)dt 

ls  Jo 


(3.27) 


In  summary,  for  a direct  compression  test,  the  average  strain,  average  stress, 
and  average  strain  rate  in  the  specimen  are  given  by  the  relations: 

2C0  rt 


Jo 


i sr  dt 


(3.28) 


= E~A£t 


(3.29) 
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. _2C0 

s 


(3.30) 


where  the  subscripts  r and  t refer  to  the  reflected  and  transmitted  signals  re- 
spectively, while  s is  in  reference  to  the  specimen  and  b the  bar.  E is  the  elastic 
modulus  of  the  bar,  Co  is  the  elastic  wave  speed  in  the  bar,  Ls  is  the  length  of 
the  specimen,  and  A refers  to  cross  sectional  area.  In  general,  the  Equations  3.28 
through  3.30  represent  the  basic  features  of  the  SHPB,  the  specimen  stress  is  pro- 
portional to  the  transmitted  strain,  the  strain  rate  is  proportional  to  the  reflected 
strain,  and  the  specimen  strain  is  proportional  to  the  integral  of  the  reflected 
strain.  It  is  important  to  note  that  these  expressions  represent  average  values  and 
are  calculated  from  uniaxial  stress-state  assumptions  in  the  specimen  and  elastic 
behavior  in  the  bars. 

3.1.3  Dynamic  Tensile  Testing 

Tensile  testing  data,  either  under  quasi-static  or  dynamic  conditions  is  usually 
obtained  by  conducting  a splitting  tension  test,  also  referred  to  as  a Brazilian 
test.  In  fact,  the  Brazilian  test  is  the  only  method  recognized  by  the  ASTM  for 
determining  quasi-static  concrete  tensile  strength.  The  required  geometry  and 
testing  procedure  is  detailed  in  ASTM  C496-96.  Figure  3-4  illustrates  the  testing 
procedure.  A cylindrical  concrete  specimen  is  subjected  to  a diametral  compressive 
force  along  its  longitudinal  axis.  The  load  is  applied  at  a rate  within  prescribed 
limits  until  specimen  failure  occurs.  The  state  of  stress  within  the  specimen  is 
characterized  by  high  compressive  forces  at  the  point  of  contact  between  the 
specimen  and  the  platens  of  the  load  frame,  and  a uniform  tensile  zone  towards 
the  center.  According  to  Neville  [15],  this  tension  zone  runs  along  approximately 
75%  of  the  diameter  of  the  specimen.  Although  the  compressive  forces  are  of  far 
greater  magnitude,  due  to  the  fact  that  concrete  is  much  weaker  in  tension  than 
in  compression,  the  specimen  fails  in  tension  before  the  compressive  loads  reach  a 
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critical  value.  Figure  3-5  illustrates  the  stress  distribution  within  the  specimen.  At 
failure,  the  magnitude  of  the  tension  force  is  given  by: 


a 


y 


2 P 
nLD 


(3.31) 


Ross  [25]  was  the  first  to  conduct  a dynamic  Brazilian  test  in  a SHPB.  Making 
the  assumption  that  the  maximum  load  sustained  in  the  specimen  is  equal  to  the 
measured  transmitted  load,  the  dynamic  tensile  strength  is: 
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(3.32) 


where  omax  is  the  maximum  sustained  tensile  load,  Rb  is  the  bar  radius,  L is  the 
specimen  length,  and  D is  the  specimen  diameter.  The  average  strain  rate  is 
determined  from  the  slope  of  the  transmitted  stress-time  signal: 


A a 

avg  = EcAt 


(3.33) 


where:  Act  is  the  change  in  stress  over  time  interval  At  and  Ec  is  the  concrete 
elastic  modulus  obtained  from  static  strength  tests.  Finite  element  analysis  was 
used  to  verify  that  the  stress  distribution  in  the  Brazilian  sample  under  dynamic 
loading  is  nearly  equivalent  to  that  in  the  static  case  [26,  27,  28]. 

Classical  Hopkinson  pressure  bar  could  also  be  used  to  generate  tensile  data 
on  concrete.  The  concrete  sample  is  cemented  to  the  end  of  the  bar  opposite  the 
firing  end.  The  load  is  then  applied  to  the  Hopkinson  bar,  usually  with  a striker 
bar  but  sometimes  through  direct  application  of  an  explosive  charge  to  the  firing 
end.  A compressive  pulse  travels  down  the  bar,  and  at  the  bar-specimen  interface, 
a portion  of  the  pulse  is  reflected  back  into  the  bar  and  a portion  transmitted  into 
the  specimen.  At  the  free  surface  of  the  specimen,  the  compressive  pulse  is  reflected 
back  as  a tensile  one.  If  the  magnitude  of  the  resulting  tensile  pulse  is  greater 
than  the  dynamic  strength  of  the  specimen,  tensile  failure,  referred  to  as  spall,  will 
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Figure  3-4:  Brazilian  test 
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Figure  3-5:  Stresses  and  stress  distribution  for  the  Brazilian  test 
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occur  trapping  the  momentum  from  the  portion  of  the  pulse  exceeding  the  dynamic 
tensile  strength  in  the  spalled  off  section.  Direct  application  of  strain  gauges  to 
the  specimen  allows  for  measurement  of  both  the  initial  compressive  pulse  and  the 
reflected  pulse  corresponding  to  the  dynamic  tensile  strength.  One  advantage  to 
this  technique  is  that  the  specimen  is  under  true  uniaxial  tension.  However,  the 
initial  compressive  load  may  also  impart  some  level  of  damage  into  the  material. 
Recent  applications  of  the  Hopkinson  bar  for  dynamic  tensile  testing  of  brittle 
materials  has  been  reported  by  Diaz- Rubio  et  al.  [29]. 

3.1.4  Confined  Dynamic  Testing 

All  of  the  techniques  discussed  so  far  have  been  used  for  dynamic  testing 
under  unconfined  conditions.  As  has  been  discussed  previously,  the  behavior  of 
most  geologic  materials  is  highly  sensitive  to  hydrostatic  pressure.  It  is  therefore 
desirable  to  gather  dynamic  data  under  confined  conditions.  One  method  of 
achieving  confining  pressure  in  a dynamic  SHPB  test  is  by  attempting  to  limit  the 
lateral  expansion  of  the  specimen.  Lankford  et  al.  [30]  Walker  et  al.  [31],  and  Chen 
and  Ravichandran  [32]  have  reported  on  efforts  to  apply  autofret.taged  (shrink  fit) 
rings  to  ceramic  specimens.  Very  high  confining  pressures  were  thus  achieved  at 
strain  rates  in  excess  of  1000  s-1.  The  disadvantages  of  this  technique  include  the 
difficulty  of  specimen  fabrication  and  the  impossibility  of  maintaining  a constant 
pressure  throughout  the  experiment.  While  confining  pressure  as  a function 
of  time  could  be  inferred  by  instrumenting  the  rings,  it  is  generally  difficult  to 
accurately  establish  what  the  initial  state  of  confinement  was.  An  alternative 
method  proposed  by  Bhushan  and  Jahsman  [33]  makes  use  of  a massive  steel  collar 
to  confine  the  specimen.  Though  specimen  preparation  appears  easier  than  in  the 
autofrettaged  approach,  due  to  tight  tolerances,  the  SHPB  specimens  had  to  be 
kept  in  a freezer  prior  to  insertion  in  the  collar.  Analysis  is  performed  under  the 
assumption  of  uniaxial  strain  conditions.  Once  again,  the  confining  pressure  on 
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the  specimen  is  not  constant  throughout  the  experiment  and  pressure  levels  must 
be  inferred.  Note  that  the  influence  of  frictional  forces  would  also  be  very  difficult 
to  quantify  in  both  techniques.  An  alternative  method  was  proposed  by  Malvern 
[34],  who  coupled  a hydrostatic  confinement  cell  to  a SffPB.  The  confinement  cell, 
shown  in  Figure  3-6  coupled  with  an  axial  preload  device  is  capable  of  placing 
the  specimen  under  true  hydrostatic  confinement  prior  to  application  of  the 
dynamic  axial  load.  Confining  pressures  were  held  fairly  constant  due  to  internal 
accumulator  action  provided  by  air  trapped  in  the  cell.  This  technique  is  quite 
close  to  a dynamic  Karman  test.  The  primary  disadvantage  of  this  technique  is 
that  the  range  of  confining  pressures  achievable  with  this  cell  is  relatively  low  (3000 
psi). 


Figure  3-6:  Hydrostatic  confinement  cell 


3.1.5  Plate  impact  Experiments 

In  a plate  impact  experiment,  a thin  disk  sample  of  concrete  is  fired  from  a 
gas  gun  and  normally  impacted  into  a hard  steel  target  plate.  Elastic  response 
of  the  steel  plate  and  uniaxial  strain  conditions  in  the  concrete  specimen  are 
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assumed.  The  transient  response  of  the  back  face  of  the  plate  is  recorded  using  a 
VISAR  [35]  laser  interferometer.  Knowing  the  density  and  acoustic  velocity  of  the 
target  material,  it  is  possible  to  calculate  the  average  stress  at  the  specimen-target 
interface  as  cr  — | pcVfs , where  cr,  V/s,  p,  and  c are  the  average  longitudinal  stress 
at  the  specimen/ target  interface,  particle  velocity  at  the  rear  surface  of  the  target, 
mass  density  and  longitudinal  elastic  wave  speed  of  the  target  material.  These  flyer 
plate  experiments  have  been  demonstrated  to  yield  constitutive  data  at  strain  rates 
from  approximately  103s_1  to  105s-1  under  confining  pressures  on  the  order  of  1 to 
1.5  GPa  [36].  Kipp  et  al.  [37]  reported  on  a series  of  reverse  ballistic  flying  plate 
experiments  in  which  WES5000  was  utilized.  In  a reverse  ballistic  experiment,  the 
concrete  material  is  the  target  plate.  Unfortunately,  since  specimen  strains  are 
not  directly  measured  in  the  flyer  plate  experiment,  it  is  not  possible  to  directly 
develop  a stress  strain  model  based  on  results  of  such  tests.  Flyer  plate  data  can 
thus  serve  only  for  validation  purposes.  When  using  the  VISAR  technique  on  a 
heterogenous  material  like  concrete,  there  appear  instrumentation  issues  as  well. 
Because  the  VISAR  measures  only  the  free  surface  response  in  a highly  localized 
area,  either  the  average  response  from  a number  of  experiments  must  be  utilized, 
or  several  VISAR  employed  in  a single  experiment  to  ensure  that  the  measurement 
is  representative  of  the  macroscopic  level  material  properties  and  not  of  local 
aggregates  or  other  local  defects. 

3.2  Dynamic  Experimental  Program 

In  order  to  dynamically  characterize  the  WES5000  concrete,  a series  of  Split 
Hopkinson  Pressure  Bar  (SHPB)  experiments  were  conducted  under  both  confined 
and  unconfined  conditions.  For  comparison  purposes,  a second  mortar  material 
was  also  characterized.  The  experimental  matrix  called  for  a total  of  three  different 
dynamic  loading  rates  at  four  different  confining  conditions.  The  experimental 
matrix  is  given  in  Table  3-1. 
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Table  3-1:  Dynamic  Karman  test  matrix 


Confined  dynamic  Karman  tests 

Con  fining  pressure 

0 

Lp 

Mp 

Hp 

Striker  bar  velocity 

C 

M 

c 

M 

C 

M 

C 

M 

Quasi  — static 

3 

3 

3 

3 

3 

3 

3 

3 

Lv 

3 

3 

3 

3 

3 

3 

3 

3 

Mv 

3 

3 

3 

3 

3 

3 

3 

3 

Hv 

3 

3 

3 

3 

3 

3 

3 

3 

Subtotal 

24 

24 

24 

24 

A total  of  96  experiments  are  included  in  this  matrix.  The  quasi-static  loading 
conditions  are  intended  to  provide  a baseline  characterization  for  comparison 
purposes.  The  confining  conditions  Lp,  Mp,  and  Hp  refer  to  low,  medium  and 
high  levels  of  confining  pressure  respectively.  Likewise,  in  the  striker  bar  velocity 
column,  Lv,  Mv,  and  Hv  stands  for  low,  medium,  and  high  velocities.  The  proposed 
matrix  represents  a large  number  of  experiments.  It  was  felt  that  this  was  required 
due  to  the  inherent  variability  of  the  concrete  specimens  as  well  as  to  the  novel 
nature  of  the  experiments.  The  duplicate  conditions  also  allowed  to  check  that  the 
data  generated  are  representative  of  the  material  behavior. 

3.3  Experimental  Setup 

3.3.1  SHPB 

The  issue  of  ensuring  that  tests  were  conducted  with  specimens  that  can 
be  assumed  to  be  homogeneous  (representative  volume  elements)  was  discussed 
earlier  in  this  chapter.  It  was  pointed  out  that  for  WES5000,  a minimum  SHPB 
diameter  of  approximately  40mm  was  required.  The  tests  were  carried  out  at  the 
Material  Properties  Laboratory  at  the  University  of  Florida’s  Graduate  Engineering 
Research  Center  (GERC)  using  a 76.2mm  diameter  SHPB.  The  GERC’s  SHPB 
is  ideally  suited  for  this  effort  and  represents  one  of  only  a few  devices  of  its  size 
available  in  the  world.  The  incident  and  transmitter  bars  are  each  3.10m  (122in) 
in  length.  Two  striker  bars  of  0.762m  (30in)  and  lm  (40in)  are  available.  They  are 
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capable  of  generating  a loading  pulse  of  approximately  300 gs  and  400 gs  duration, 
respectively.  The  incident  and  transmitter  pressure  bars,  as  well  as  the  strikers 
are  fabricated  from  ETD-150  steel,  whose  density,  bar-wave  speed,  and  Poisson’s 
ratio  are  7,650  kg/m3,  5,150  m/s  and  0.32,  respectively.  The  strikers  are  launched 
with  a gas  gun,  fabricated  by  Terra- Tek  Systems,  Inc.  The  firing  chamber  of  the 
gun  has  a maximum  chamber  pressure  of  4.1  MPa  ( 2000  psi).  The  propellant 
gas  is  furnished  by  a nitrogen  bottle.  An  optical  device  for  measuring  striker  bar 
velocity  is  positioned  at  the  mouth  of  the  gas  gun  barrel.  A full  strain-gauged 
bridge  is  permanently  mounted  at  the  midpoint  of  each  pressure  bar.  Strain 
gauge  conditioning  is  provided  via  Measurements  Group  model  2311  strain  gauge 
conditioning  amplifiers  with  a bandwidth  of  1.25  MHz.  Data  signals  are  recorded 
by  a Nicolet  4094  digital  oscilloscope  capable  of  recording  at  0.5  gs  intervals.  The 
GERC’s  SHPB  is  shown  in  Figure  3-7. 


Figure  3-7:  University  of  Florida  GERC  SHPB 
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Radial  confinement  may  be  applied  to  the  material  sample  using  a confining 
cell  originally  developed  by  Malvern  [34] . The  confinement  cell  is  constructed  of 
17-4  PH  Stainless  Steel  (Condition  A)with  a maximum  operating  pressure  of  20.7 
MPa  (3000  psi).  The  design  includes  a factor  of  safety  of  10.  The  cell  contains  4 
access  ports,  including  one  which  accommodates  a pressure  gauge.  A cut  away 
drawing  of  the  pressure  cell  is  shown  in  Figure  3-8 


Figure  3-8:  Cut  away  view  of  GERC  confinement  cell  (all  dimensions  in  inches) 

In  order  to  measure  the  pressure  in  the  confining  cell,  a Kulite  model  HKM- 
375-5000-A  transducer,  with  a natural  frequency  of  1.12  MHz,  was  selected  for  the 
study.  To  achieve  a triaxial  pre-test  load  on  the  test  specimens,  a hand  operated 
hydraulic  load  cell  was  utilized.  Figure  3-  9 shows  the  load  cell  in  position  at 
the  far  end  of  the  transmitter  bar.  In  order  to  restrain  the  incident  bar,  a “bar 
restraint”  was  designed  and  positioned  near  the  mouth  of  the  gas  gun.  The  device 
consisted  of  a steel  plate  with  a counter-sunk  hole  cut  into  it  (see  Figure  3-10).  A 
small  bar  or  “stop”  was  machined  from  ETD-150  steel  to  fit  into  the  counter-sunk 
hole.  The  smaller  inner  diameter  of  the  stop  matched  the  outer  diameter  of  both 
the  striker  and  incident  bars.  Figures  3 11  and  3 12  show  the  stop  and  the  end 
restraint  with  the  stop  in  place,  respectively. 
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Figure  3-9:  Hydraulic  load  cell  for  SHPB 


Figure  3-10:  End  restraint  for  SHPB 
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Figure  3-11:  Bar  stop  for  SHPB 
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Figure  3-12:  End  restraint  with  stop  in  place  for  SHPB 
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3.3.2  Quasi-Static  Setup 

For  the  baseline  quasi-static  characterization,  an  MTS  815  rock  mechanics 
system  load  frame  was  utilized.  Since  the  load  frame  did  not  possess  a confinement 
cell,  a means  was  required  to  adapt  the  SHPB  confinement  cell  to  a conventional 
hydraulic  load  frame.  Two  attachments  were  designed  for  the  MTS  machine, 
essentially  consisting  of  high-strength  tool  steel  bars  with  special  attachments  for 
connecting  to  the  upper  and  lower  platens  of  the  load  frame.  The  adaptors  are 
shown  in  position  on  the  MTS  machine  with  the  SHPB  cell  attached  in  Figure 
3-13. 

3.3.3  Specimen  Design/Preparation 

According  to  Gray  [38],  radial  and  longitudinal  inertia  errors,  as  well  as 
friction  effects  can  be  reduced  by  choosing  a sample  size  that  minimizes  the  areal 
mismatch  between  the  sample  and  the  bar  diameter,  and  maintains  a sample  L/D 
ratio  of  0.5  to  1.0.  It  is  also  desirable  from  an  homogeneity  standpoint  to  have  an 
equal  number  of  aggregates  across  both  the  specimen  length  and  diameter.  For 
these  reasons,  a nominal  L/D  ratio  of  1-to-l  was  selected.  TheWES5000  samples 
were  cored  from  the  same  material  batch  utilized  in  the  quasi-static  experiments 
of  Chapter  2.  The  cores  taken  were  nominally  75  mm  in  diameter  and  were  cut 
into  lengths  of  approximately  76  mm.  The  end  faces  on  each  sample  were  then 
ground  to  assure  parallelism.  The  diameter  and  length  of  each  specimen  was 
measured  with  a set  of  digital  calipers  and  weighed  on  a digital  balance.  The 
samples  utilized  were  screened  according  to  end  face  parallelism,  roundness, 
and  obvious  large  surface  flaws.  Due  to  the  porous  nature  of  WES5000,  surface 
preparation  was  necessary  for  both  strain  gage  application  and  to  ensure  that  the 
latex  membranes  utilized  to  isolate  the  specimens  from  the  fluid  in  the  confined 
experiments  were  not  forced  into  pores  or  voids  causing  punctures.  Preparing  the 
samples  was  a multi-step  process.  The  specimen  was  first  probed  with  a pick  to 
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Figure  3-13:  Load  frame  modified  for  confined  test 
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identify  and  open  larger  surface  voids.  These  voids  were  then  filled  with  a gypsum 
paste.  The  use  of  a gypsum  past  was  recommended  by  Waterways  Experimental 
Station  and  it  was  believed  it  would  not  influence  the  strength  properties  of  the 
material.  The  paste  was  then  sanded  down  even  with  the  surface  of  the  specimen 
and  excess  dust  removed  with  a damp  cloth.  Prior  to  strain  gage  instrumentation, 
an  initial  application  of  epoxy  is  required.  This  initial  application  is  sanded  down 
smooth  and  the  specimen  again  cleaned  and  neutralized.  A second  application 
of  epoxy  is  then  applied  for  actual  gage  application.  However,  trial  showed  that 
the  epoxy  tended  to  be  absorbed  into  the  gypsum  and  upon  drying  resulted  in 
depressions  in  the  specimen  unsuitable  for  gage  application.  In  order  to  avoid  this, 
an  intermediate  layer  of  cement  epoxy  was  applied.  This  was  then  sanded  down 
smooth,  cleaned  as  called  for,  and  the  gage  epoxy  applied.  It  was  felt  that  this 
would  not  affect  the  overall  material  behavior.  Following  strain  gauge  application, 
clamps  were  placed  on  the  specimens  to  apply  a uniform  pressure  on  the  gauge. 
The  specimens  were  then  allowed  to  cure  for  several  days  in  a warm  area  prior 
to  conducting  the  tests.  For  comparison  purposes,  mortar  was  also  tested.  This 
material  was  formulated  and  prepared  for  the  study  by  the  Air  Force  Research 
Laboratory  Materials  Directorate’s  Pavement  Division  at  Tyndall  Air  Force  Base. 
Samples  of  this  material  were  cast  into  molds  on  8 May  1998.  Mixture  properties 
for  the  mortar  are  provided  in  Table  3-2. 

Table  3-2:  Ingredients  and  mixture  proportions  for  mortar 


Item 

Mixture  proportions, 
saturated  surface-dry 

Type  I Portland  cement 

450# 

Sand,  sieve  #4 

1620s 

WRDA-19  0.5% 

3.7  g 

“F”  fly  ash 

297 g 

Water,  w/c  = 0.55 

411  g 
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Measurements  Group  precision  strain  gauges  were  utilized  in  the  experiments. 
Gauge  type  EA-06-10CBE-120,  with  a gauge  length  of  one-inch  was  selected 
for  the  WES5000  to  average  out  the  local  aggregate  effects,  while  type  EA-06- 
0500BH-120,  with  a gauge  length  of  0.5  inch  was  selected  for  the  mortar  because 
the  mix  doesn’t  contain  any  aggregates.  Two  axial  and  two  radial  gauges  were 
applied  to  each  specimen.  The  application  of  both  axial  and  radial  gauges  allowed 
measurement  of  volume  changes  under  dynamic  conditions.  Figures  314  and  3 15 
show  instrumented  WES5000  and  a mortar  samples. 


Figure  3-14:  Strain  gauges  applied  to  WES5000  SHPB  specimen 
3.3.4  Experimental  Protocol 

In  this  section,  the  protocol  observed  for  the  confined  SHPB  experiments 
is  presented.  First,  special  care  is  required  to  ensure  that  the  striker  bar  is  fully 
retracted  into  the  barrel  of  the  gas  gun.  Next,  the  bar  stop  is  inserted  into  the 
restraint  and  the  incident  bar  moved  towards  the  stop  until  the  two  are  in  contact. 
The  UF/GERC  confining  cell,  along  with  the  locking  flange  o-ring  and  collar  are 
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Figure  3-15:  Mortar  SHPB  specimen  with  strain  gauges  applied 

then  positioned  onto  the  transmitter  bar.  Proper  orientation  is  required  in  that  the 
incident  side  flange  is  of  slightly  greater  inner  diameter  to  allow  feed-through  of  the 
instrumentation  wires,  and  the  respective  flanges  are  matched  to  specific  sides  of 
the  cell.  Next,  thin  ETD-150  steel  disks  of  the  same  outer  diameter  as  the  pressure 
bars  are  placed  on  both  the  incident  and  transmitter  faces  and  held  in  place  with  a 
piece  of  electrical  tape.  A thin  layer  of  grease  is  applied  to  the  interface  to  assure 
good  contact  is  achieved.  The  disks  serve  to  protect  the  incident  and  transmitter 
bar  faces  from  any  damage  due  to  contact  with  the  hard  aggregate.  The  disks  were 
polished  after  each  test  to  remove  any  burrs  or  scratches.  A latex  membrane  is 
then  placed  over  the  transmitter  bar  and  secured  at  one  end  with  electrical  tape. 
The  latex  membrane  is  to  be  utilized  as  a barrier  between  the  specimen  and  the  the 
confining  fluid.  Figure  3-16  shows  this  first  stage  of  experimental  setup. 

The  specimen  is  next  placed  between  the  incident  and  transmitter  bars.  To 
minimize  frictional  effects,  a thin  layer  of  grease  is  applied  to  the  face  of  the  disks 
covering  the  bar  faces.  The  specimen  is  held  in  place  with  one  wrap  of  electrical 
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Figure  3-16:  Setup:  cell  in  place  on  SHPB 


tape  at  each  end  (Figure  3-17).  The  strain  gauges  are  then  connected  to  electrical 
leads  made  of  34-gage  solid  core  copper  wire.  A very  fine  wire  was  necessary  so 
that  it  would  pass  between  the  incident  side  flange  and  the  pressure  bar.  Special 
care  was  required  so  as  not  to  tangle  or  kink  the  very  hue  wires.  The  wires  are  led 
down  the  bar  and  held  in  place  with  a strip  of  electrical  tape.  A strip  of  tape  was 
also  placed  over  each  gauge  (Figure  3-18).  Finally,  the  latex  membrane  is  pulled 
over  the  specimen  and  secured  in  place  with  a piece  of  electrical  tape  (Figure  3 19). 

The  confinement  cell  is  then  carefully  slid  into  position  over  the  specimen. 

Due  to  close  tolerances,  a water  based  lubricant  is  first  applied  to  the  membrane 
(Figure  3-20).  The  o-rings  and  flange  are  then  placed  in  position.  Locking  collars 
are  threaded  onto  the  cell  and  tightened  with  a set  of  specially  designed  spanner 
wrenches.  The  pressure  applied  to  the  flange,  and  subsequently  the  o-rings,  by 
the  locking  collar  causes  the  o-rings  to  seal  around  the  SHPB.  The  pressure  gauge 
may  now  be  threaded  into  the  instrumentation  port  and  securely  tightened. 
Pressurization  is  applied  to  the  confining  cell  via  a 20.7  MPa  (3000  psi)  capacity 
hydraulic  accumulator.  The  accumulator  is  next  filled  with  water  and  connected 
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Figure  3-  17:  Setup:  specimen  in  place  and  secured 


to  the  confinement  cell  through  an  opening  in  the  SHPB  support  structure  (Figure 
3 21).  The  pressurization  fluid  selected  was  water  because  the  latex  membrane 
would  have  degraded  through  contact  with  normal  petroleum  based  hydraulic  fluid. 
Water  is  forced  into  the  confinement  cell  when  the  accumulator  is  pressurized.  The 
air  trapped  in  the  cell  was  expected  to  provide  internal  accumulator  action  for 
damping  of  the  pressure  during  the  actual  test  [34].  The  external  accumulator  was 
not  expected  to  provide  any  pressure  damping  since  its  response  time  was  believed 
far  too  slow  considering  that  the  duration  of  the  experiment  was  on  the  order  of 
a few  hundred  micro-seconds.  It  did  prove  however  to  be  an  excellent  means  of 
maintaining  a nearly  constant  pressure  prior  to  launching  the  striker  bar  as  well  as 
compensating  for  minor  leaks  prior  to  testing.  The  accumulator  is  charged  using 
a high  pressure  hose  connected  to  a high  pressure  nitrogen  bottle.  Pressurization 
of  the  system  is  a multi-step  process.  A small  axial  load  is  applied  using  the 
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Figure  3 18:  Setup:  electrical  leads  connected 


hydraulic  hand  pump,  then  the  accumulator  is  pressurized  via  the  nitrogen  tank 
to  the  same  level.  The  axial  load  is  monitored  with  a pressure  gauge  (Kulite 
HKS- 11-375-2000  SG)  attached  to  the  pump  and  connected  to  a digital  readout. 
The  cell  pressure  was  monitored  by  a pressure  gauge,  whose  output  was  read  on 
one  of  the  oscilloscopes.  It  was  necessary  to  apply  the  axial  load  first  to  assure 
that  the  radial  pressure  did  not  force  the  bars  apart,  causing  the  specimen  to  lose 
contact  with  the  SHPB.  By  continuing  the  step-wise  process  of  increasing  the  axial 
load,  then  matching  the  radial  pressure,  the  specimen  could  be  safely  brought  up 
to  the  desired  hydrostatic  load.  Once  the  specimen  is  under  the  desired  level  of 
hydrostatic  confinement,  the  strain  gauges  are  balanced  and  the  instrumentation 
set  to  record.  The  gas  gun  could  then  be  charged  to  the  desired  pressure  and  fired. 
Data  was  recorded  on  the  Nicolet  4094  B oscilloscopes. 
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Figure  3-19:  Setup:  latex  membrane  in  place  and  secured 


Setup  for  the  quasi-static  experiments  utilizing  the  MTS  machine  was  almost 
identical.  The  only  exception  is  that  the  step  wise  hydrostatic  loading  was  done  by 
programming  the  MTS  machine  . Once  the  desired  level  of  hydrostatic  confinement 
was  achieved,  it  was  kept  constant  and  the  axial  load  was  applied  via  the  platens  of 
MTS  machine.  Data  was  recorded  with  a computer  based  data  acquisition  system. 

3.4  Experimental  Results 

3.4.1  Quasi-Static  Baseline 

Characterization  experiments  were  conducted  on  mortar  and  WES5000 
samples  under  quasi-static  loadings  conditions  (strain  rate  ~ 10-6)  for  four 
different  confining  pressures.  The  specimen  sizes  were  of  the  same  dimensions  as 
those  used  for  the  SHPB  tests.  Each  specimen  was  strain  gauged  both  axially 
and  radially  (two  each)  so  that  the  volumetric  response  could  be  measured.  Each 
experiment  was  given  a unique  designator,  either  QSMCXY,  QSMUXX,  QSCCXY, 
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Figure  3-20:  Setup:  confinement  cell  in  position  over  specimen 


Figure  3-21:  Setup:  accumulator  and  end  caps  connected 


91 


or  QSCUAA,  where  QSMC,  QSMU,  QSCC  and  QSCU  refer  to  quasi-static  mortar 
confined,  quasi-static  mortar  unconfined,  quasi-static  concrete  confined  and  quasi- 
static concrete  unconfined  respectively.  The  XX  is  a unique  numerical  designator. 
The  following  two  sections  summarize  the  quasi-static  results  for  mortar  and 
WES5000,  respectively. 

Mortar 

A total  of  13  quasi-static  experiments  were  conducted  on  mortar  samples 
under  both  confined  and  unconfined  conditions.  The  confining  pressures  corre- 
sponding to  Lp,  Mp,  and  Hp  in  Table  3-1  are  1.72  MPa  (250  psi),  3.45  MPa  (500 
psi)  and  6.89  MPa  (1000  psi)  respectively.  Nominal  loading  rate  was  133  New- 
tons/s (30  lbf/s).  The  experiment  identifier,  confining  condition  under  which  they 

were  run,  and  the  peak  principal  stress  difference  are  summarized  in  Table  3-3. 

Table  3-3:  Quasi-static  mortar  results 


Experiment  # 

Confining  pressure 
(MPa) 

Qmax  {MPa) 

Comments 

QSMU01 

- 

49.5 

QSMU02 

- 

50.0 

QSMU03 

- 

56.0 

unload-creep-reload  cycles 

QSMU04 

- 

54.5 

unload-creep-reload  cycles 

QSMC01 

3.45 

70.2 

QSMC02 

3.45 

63.1 

QSMC03 

3.45 

63.8 

QSMC04 

3.45 

69.0 

QSMC05 

1.72 

63.6 

QSMC06 

1.72 

59.0 

QSMC07 

6.89 

78.1 

QSMC08 

6.89 

64.0 

membrane  failure 

QSMC09 

6.89 

81.6 

Figures  3-22  through  3-25  show  the  principal  stress  difference  versus  strain 
curves  for  the  various  loading  conditions.  Only  one  representative  curve  is  shown 
for  each  condition.  The  plots  show  initial  compressibility,  followed  by  a gradual 
transition  to  dilatancy.  The  dilatancy  threshold  depends  on  the  confinement,  the 
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higher  the  confinement,  the  closer  is  the  dilatancy  threshold  to  the  peak  stress 
reached  in  the  test.  Representative  principal  stress  difference  vs.  volume  strain 
curves  are  presented  for  the  four  levels  of  confinement  in  Figure  3-26.  Here  the 
effect  of  confinement  on  strength  and  dilatancy  can  be  clearly  seen.  The  right 
axis  of  Figure  3-26  is  expressed  in  terms  of  Confinement  Increase  Factor,  or  CIF. 
The  CIF  is  defined  as  the  principal  stress  difference  normalized  by  the  unconfined 
compressive  strength.  Even  at  the  relatively  low  confinement  level  of  6.89  MPa, 
a strength  enhancement  of  almost  1.5  times  is  observed.  The  average  values  of 
the  peak  stresses  (strength)  are  shown  in  Table  3-4  and  a plot  of  peak  stresses 
versus  confinement  pressure  is  shown  in  Figure  3-27.  Note  that  over  the  confining 
pressure  range  investigated,  the  relationship  between  strength  and  confinement 
pressure  is  nearly  linear. 


Figure  3-22:  Quasi-static  mortar  unconfined  compression  experiment  - QSMU01 


Principal  Stress  Difference  (MPa)  co  Principal  Stress  Difference  (MPa) 
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23:  Quasi-static  mortar  confined  (QSMC05)  - 1.72  MPa  confinement 


Figure  3-24:  Quasi-static  mortar  confined  (QSMC04)  - 3.45  MPa  confinement 
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Table  3-4:  Average  peak  principal  stress  difference  at  failure  for  mortar. 


Confining  pressure  (MPa) 

0 

1.72 

3.45 

6.89 

Avg.  peak  principal  stress  diff.  (MPa) 

52.5 

61.3 

66.5 

79.85 

100 


-0.004  -0.002  0.000  0.002  0.004  0.006  0.008  0.010 


Strain  (mm/mm) 


Figure  3-25:  Quasi-static  mortar  confined  (QSMC07)  - 6.89  MPa  confinement 


Figure  3-26:  Confinement  effect  on  quasi-static  mortar  behavior 


CIF  (MPa/MPa) 
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Figure  3-27:  Peak  principal  stress  difference  at  failure  as  a function  of  confinement 
for  mortar 
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WES50Q0 

A total  of  14  quasi-static  experiments  were  conducted  on  WES5000  samples 
under  both  confined  and  unconfined  conditions.  The  confining  pressures  were 
the  same  as  those  utilized  in  the  quasi-static  mortar  experiments  . The  nominal 
loading  rate  was  133  Newtons/s  (30  lbj/s).  The  experiment  identifiers,  confining 
conditions  under  which  they  were  run,  and  the  peak  principal  stress  differences  are 

summarized  in  Table  3-5. 

Table  3-5:  Quasi-static  WES5000  results 


Experiment  # 

Confining  pressure  (MPa) 

Qmax  ( MPd ) 

Comments 

QSCU01 

- 

62.9 

QSCU02 

- 

64.9 

QSCC01 

3.45 

107.7 

QSCC02 

3.45 

100.5 

QSCC03 

3.45 

100.9 

QSCC04 

3.45 

112.4 

QSCC05 

6.89 

68.5 

membrane  failure 

QSCC06 

6.89 

120.7 

QSCC07 

6.89 

41.0 

membrane  failure 

QSCC08 

6.89 

75.3 

membrane  failure 

QSCC09 

1.72 

52.5 

membrane  failure 

QSCC10 

1.72 

48.6 

membrane  failure 

QSCC11 

1.72 

86.7 

QSCC12 

1.72 

91.7 

Performing  these  experiments  with  WES5000  proved  to  be  very  problematic. 
The  one  successful  experiment  at  6.89  MPa  (1000  psi)  required  68  minutes  to 
complete.  In  Chapter  2,  it  was  briefly  discussed  that  because  of  the  heterogeneous 
nature  of  WES5000,  under  high  confinement,  strain  gage  measurements  are 
nearly  impossible.  The  same  problems  were  encountered  during  testing  under 
very  low  confining  pressures.  As  the  experiment  progressed,  and  the  axial  load 
on  the  specimen  increased,  the  likelihood  of  a membrane  failure  and  gage  failure 
dramatically  increased.  Due  to  the  tight  tolerances  of  the  GERC  confinement 
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cell,  the  use  of  multiple  membranes  and  layers  of  putty  surrounding  the  specimen 
(the  technique  utilized  by  WES  and  discussed  in  Chapter  2)  was  not  an  option. 
However,  the  correlation  between  the  successful  experiments  was  quite  good. 

Figures  3-28  through  3-31  present  the  principal  stress  difference  versus  strain 
curves  for  the  various  loading  conditions.  Only  one  representative  plot  is  shown 
for  each  condition.  The  data  generated  at  the  GERC  is  consistent  with  the  data 
generated  at  WES  for  unconfined  conditions  and  confining  pressures  below  50 
MPa  (see  Chapter  2).  The  curves  show  initial  compaction,  followed  by  a gradual 
transition  to  dilatancy  with  the  transition  to  dilatancy  moving  closer  to  the  failure 
strength  as  confinement  is  increased.  Representative  principal  stress  difference  vs. 
volume  strain  curves  are  presented  for  the  four  levels  of  confinement  in  Figure  3 32. 
Here  the  effect  of  confinement  on  strength  and  dilatancy  can  be  clearly  seen.  For 
confinement  levels  of  6.89  MPa,  C1F  values  of  almost  2 were  observed.  Average 
values  of  peak  principal  stress  difference  were  taken  for  each  loading  condition. 

The  average  values  are  shown  in  Table  3-6  and  are  plotted  in  Figure  3-33.  Over 
the  confining  pressure  range  investigated,  the  relationship  between  strength  and 
confinement  pressure  is  nearly  linear. 

Table  3-6:  Average  peak  principal  stress  difference  at  failure  for  WES5000 


Confining  pressure  (MPa) 

0 

1.72 

3.45 

6.89 

Avg.  peak  principal  stress  diff.  (MPa) 

63.9 

89.2 

105.4 

120.7 
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Figure  3-28:  Quasi-static  concrete  unconfined  (QSCU01) 


Figure  3-29:  Quasi-static  concrete  confined  (QSCCll) 


F72  MPa  confinement 
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Figure  3-30:  Quasi-static  concrete  confined  (QSCC04)  - 3.45  MPa  confinement 


Figure  3-31:  Quasi-static  concrete  confined  (QSCC06)  - 6.89  MPa  confinement 
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Figure  3-32:  Confinement  effect  on  quasi-static  WES5000  behavior 


130 


2 4 6 
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Figure  3-33:  Peak  principal  stress  difference  at  failure  as  a function  of  confinement 
for  WES5000 


CIF  (MPa/MPa) 
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3.4.2  SHPB  Results 

In  order  to  characterize  the  dynamic  behavior  of  WES5000,  and  the  mortar 
material,  a series  of  unconfined  SHPB  and  confined  experiments  at  1.72  MPa 
(250  psi),  3.45  MPa  (500  psi),  and  6.89  MPa  (1000  psi)were  conducted  following 
the  protocol  outlined  in  Section  (3.3.4)  of  this  chapter.  The  range  of  striker  bar 
velocities  was  selected  based  on  the  capacity  of  the  gas  gun  firing  chamber.  The 
highest  velocity  was  taken  to  be  the  speed  obtained  by  the  0.762m  bar  when 
the  gun  was  pressurized  to  3.45  MPa,  the  lowest  was  selected  to  correspond  to  a 
chamber  pressure  of  1.39  MPa  (200  psi),  while  a mid-range  velocity  corresponded 
to  2.41  MPa  (350  psi).  For  convenience,  the  low,  medium,  and  high  confining 
conditions  are  referred  to  as  Lp,  Mp,  and  Hp , respectively  while  the  low,  medium, 
and  high  velocity  conditions  are  denoted  by  Lv,  Mv,  and  Hv. 

Dynamic  Characterization  of  Mortar 

A total  of  36  SHPB  experiments  were  performed  on  the  mortar  material 
under  confined  and  unconfined  conditions.  The  test  conditions  are  given  in  Table 
3-7.  DCMC  and  DCMU  stands  for  Dynamic  Compression  Mortar  Confined  and 
Unconfined  respectively.  Each  experiment  was  also  given  a two  digit  numerical 
identifier.  Experiments  DCMC01  through  DCMC12  were  performed  for  the 
purposes  of  developing  the  experimental  protocol  and  the  results  are  not  presented 
here. 

A one-wave  analysis  was  utilized,  thus  the  peak  principal  stress  difference 
values  (referred  to  as  qmax ) were  derived  solely  from  the  transmitted  pulse.  Since 
strains  were  measured  directly  on  the  specimens,  there  was  no  need  to  time  shift 
and  combine  the  incident  and  reflected  pulses  to  determine  the  axial  strains.  Due 
to  the  large  number  of  tests  there  was  a great  deal  of  motivation  to  develop  a 
technique  for  partially  automating  the  process.  Most  existing  techniques  rely 
on  the  precise  timing  of  the  pulses.  Kaiser  [39]  developed  a technique  based  on 
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Table  3-7:  Mortar  SHPB  test  conditions 


Experiment  # 

Design 

condition 

Striker  (m) 

Pg(MPa) 

Pc(MPa) 

Date 

DCMC13 

LvMp 

0.762 

1.38 

3.45 

8-Mar- 1999 

DCMC14 

LvMp 

0.762 

1.38 

3.45 

8-Mar- 1999 

DCMC15 

LVMP 

0.762 

1.38 

3.45 

17-Mar- 1999 

DCMC16 

HyMp 

0.762 

3.45 

3.45 

17-Mar- 1999 

DCMC17 

HVMP 

0.762 

3.45 

3.45 

24-Mar- 1999 

DCMC18 

HVMP 

0.762 

3.45 

3.45 

29-Mar- 1999 

DCMC19 

MVMP 

0.762 

2.41 

3.45 

2- Apr- 1999 

DCMC20 

MVMP 

0.762 

2.41 

3.45 

2-Apr-1999 

DCMC21 

MVMP 

0.762 

2.41 

3.45 

4- Apr- 1999 

DCMC22 

MVHV 

0.762 

2.41 

6.14 

4- Apr- 1999 

DCMC23 

Ly  H p 

0.762 

1.38 

6.89 

29-Apr-1999 

DCMC24 

LyHp 

0.762 

1.38 

6.89 

30- Apr- 1999 

DCMC25 

HVHP 

0.762 

3.45 

6.89 

3-May- 1999 

DCMC26 

HVHP 

0.762 

3.45 

6.89 

3-May- 1999 

DCMC27 

HyLp 

0.762 

3.45 

1.72 

3-May- 1999 

DCMC28 

HyLp 

0.762 

3.45 

1.72 

1-May- 1999 

DCMC29 

LyHp 

1.0 

2.07 

6.89 

13-May- 1999 

DCMC30 

LyHp 

1.0 

2.07 

6.89 

14-May- 1999 

DCMC31 

MvLp 

1.0 

4.14 

1.72 

14-May- 1999 

DCMC32 

MvLp 

1.0 

4.27 

1.72 

17-May- 1999 

DCMC33 

MVHP 

1.0 

4.27 

6.89 

20-May- 1999 

DCMC34 

MyHp 

1.0 

4.27 

6.89 

24-May- 1999 

DCMC35 

MyHp 

1.0 

4.27 

6.89 

25-May- 1999 

DCMC36 

LyLp 

1.0 

2.07 

1.72 

25-May-1999 

DCMC37 

Ly  Lp 

1.0 

2.07 

1.72 

26-May- 1999 

DCMC38 

Ly  Lp 

1.0 

2.07 

1.72 

27-May- 1999 

DCMC39 

MVLP 

1.0 

4.27 

1.72 

27-May- 1999 

DCMU01 

Ly 

0.762 

1.03 

- 

23-Sep-1998 

DCMU02 

Ly 

0.762 

1.03 

- 

23-Sep-1998 

DCMU03 

Ly 

0.762 

1.38 

- 

22-Oct-1998 

DCMU04 

Ly 

0.762 

1.38 

- 

26-Oct-1998 

DCMU05 

Hv 

0.762 

3.45 

- 

19-Nov-1998 

DCMU06 

Mv 

0.762 

2.41 

- 

7-Feb-2000 

DCMU07 

Mv 

0.762 

2.41 

- 

7-Feb-2000 

DCMU08 

Mv 

0.762 

2.41 

- 

7-Feb-2000 

DCMU09 

Hv 

0.762 

3.45 

- 

7-Feb-2000 
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taking  numerical  derivatives  of  the  signal  and  applying  a rule  regarding  a minimum 
number  of  same  sign  derivatives  to  identify  the  start  of  the  pulse.  The  technique 
utilized  in  this  research  was  a statistical  one.  The  leading  pre-trigger  signal 
background  noise  was  evaluated  and  the  standard  deviation  and  mean  determined. 
Using  the  mean,  the  signal  could  be  automatically  base-line  shifted,  accounting  for 
any  instrumentation  offset.  Based  on  the  background  noise  analysis,  the  incident 
and  transmitted  signals  were  compared  to  the  standard  deviation.  A consistent 
positive  deviation  from  the  standard  deviation  was  determined  to  indicate  the  start 
of  the  pulses.  Since  the  mid-point  of  the  SHPB  specimen  was  half  way  between  the 
two  SHPB  strain  gauges,  the  specimen  strain  gauges  were  set  to  zero  corresponding 
to  a point  half  the  difference  between  the  start  of  the  incident  and  transmitted 
pulse.  The  results  of  the  SHPB  mortar  experiments  are  summarized  in  Table  3-  8. 
The  column  heading  labeled  “Comments”  always  contains  a letter  which  is  either 
A,  B,  C,  of  D.  These  letter  designators  are  an  attempt  to  provide  a qualitative 
measure  of  the  damage  observed  in  the  specimen  post  test.  A designation  of 
“A”  means  that  there  was  either  no,  or  only  superficial  visible  damage  to  the 
specimen,  “B”  corresponded  to  cracking  in  the  specimen.  Crack  orientation  was 
generally  at  a 45-degree  orientation  to  the  vertical  axis,  indicative  of  shear  failure. 

A damage  state  of  “C”  was  meant  to  indicate  that  the  specimen  was  no  longer 
capable  of  sustaining  any  appreciable  load.  This  generally  meant  that  the  specimen 
was  in  several  pieces  post  test.  The  designation  of  “D”  was  indicative  of  total 
material  failure,  and  was  only  observed  in  the  unconfined  experiments.  Figures 
3-34  through  3-36  are  representative  of  the  visible  damage  observed  post  test 
for  mortar.  There  were  no  instances  of  the  severe  level-C  failure  in  the  mortar 
experiments.  In  fact,  most  of  the  specimens  showed  very  little  damage  post  test 
other  than  some  chipping  around  the  top  edges,  most  likely  due  to  localized  tensile 


effects. 
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Table  3-8:  Mortar  SHPB  results 


Experiment  # 

V,  {m/s) 

e(s  !) 

Qmax  [MPa) 

Comments 

DCMC13 

- 

- 

- 

A (no  data) 

DCMC14 

- 

74.9 

40.0 

A 

DCMC15 

- 

67.1 

62.34 

A 

DCMC16 

- 

120 

94.1 

B 

DCMC17 

- 

117 

117.97 

A 

DCMC18 

- 

130 

107.5 

A 

DCMC19 

11.14 

110 

96.0 

B 

DCMC20 

11.23 

90 

106.7 

A 

DCMC21 

11.23 

105 

101.2 

B (no  press.) 

DCMC22 

7.79 

62 

88.75 

A 

DCMC23 

7.6 

54 

90.3 

A 

DCMC24 

7.76 

60 

93.4 

A 

DCMC25 

13.78 

122 

136.9 

A 

DCMC26 

13.71 

- 

- 

B (no  press.) 

DCMC27 

13.61 

120 

104 

B 

DCMC28 

13.93 

131 

103.6 

B 

DCMC29 

7.51 

- 

- 

A 

DCMC30 

- 

41 

107.6 

A 

DCMC31 

10.8 

87 

107.6 

B 

DCMC32 

11.01 

51 

95.9 

B 

DCMC33 

11.31 

48.7 

115.6 

A 

DCMC34 

11.19 

38.3 

117 

A 

DCMC35 

11.04 

54 

113.9 

A 

DCMC36 

7.23 

46 

89.3 

A 

DCMC37 

7.27 

46.2 

76.7 

B 

DCMC38 

7.29 

40 

84.0 

B 

DCMC39 

11.16 

88 

100.2 

B 

DCMU01 

- 

71.2 

70.18 

D 

DCMU02 

- 

65.2 

65 

D 

DCMU03 

- 

75 

75.5 

D 

DCMU04 

- 

60 

76.9 

D 

DCMU05 

- 

156 

98.6 

D 

DCMU06 

11.36 

108 

117.1 

D 

DCMU07 

11.32 

- 

- 

D 

DCMU08 

11.39 

Ill 

113.6 

D 

DCMU09 

14.01 

125 

126.6 

D 
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Figure  3-34:  Damage  state  A - DCMC14 


Figure  3-35:  Damage  state  B - DCMC27 
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In  order  to  evaluate  the  influence  of  loading  rate  and  confinement  on  the 
materials  volumetric  behavior  and  strength,  plots  were  made  for  each  experiment 
of  axial,  radial,  and  volume  strain  as  a function  of  principal  stress  difference,  or  q. 
Representative  examples  for  each  loading  case  are  presented  below.  Figure  3-37 
shows  the  results  of  the  unconfined  test  DCMU04,  conducted  under  low  striker 
bar  velocity  conditions.  The  strain  rate  determined  from  the  reflected  pulse  was 
60/s.  The  response  of  the  specimen  is  qualitatively  similar  to  the  unconfined 
response,  in  that  it  exhibits  both  compressive  and  dilatant,  behavior.  Significant 
strength  enhancement  over  the  unconfined  compressive  strength  is  also  observed. 
Only  the  data  recorded  at  25/us  intervals  are  shown  in  the  plots.  According  to 
Davies  and  Hunter  [22]  n reverberations  are  required  for  the  specimen  to  “ring 
up”  and  achieve  a uniform  stress  state.  Since  the  elastic  wave  speed  of  mortar 
is  approximately  4.8  km/s,  about  100  /us  would  be  required  for  the  specimen 
to  achieve  a uniform  stress  state.  Thus,  for  the  lowest  loading  rate  condition, 
stress  uniformity  can  be  assumed  only  at  the  onset  of  failure.  Therefore,  from  the 
unconfined  experiments,  only  information  concerning  the  peak  load  as  a function 
of  strain  rate  may  be  deemed  reliable.  Volumetric  strain  curves  as  a function  of 
principal  stress  difference  for  unconfined  conditions  and  at  three  different  strain 
rates  (from  quasi-static  to  160/s)  are  presented  in  Figure  3-38.  The  concept  of  a 
Dynamic  Increase  Factor  (DIF)  is  introduced  here  to  help  quantify  the  influence  of 
loading  rate  on  strength.  The  DIF  is  defined  as  the  dynamic  strength  normalized 
by  the  unconfined  compressive  strength  [40].  The  DIF  is  given  on  the  right  hand 
axis  of  Figure  3-38.  For  the  examples  shown,  over  the  range  of  strain  rates,  the 
material  continues  to  exhibit  both  compressible  and  dilatant  behavior.  Strength 
enhancement  is  obvious  with  a DIF  value  of  approximately  1.8  apparent  for  the 
160/s  strain  rate  curve.  The  transition  from  compressibility  to  dilatancy  also 
appears  to  be  rate  sensitive. 
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Figure  3-36:  Damage  state  D - DCMU07 


Figure  3-37:  Direct  compression  mortar  unconfined  (DCMU04)  - Lv  condition 
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Next  presented  are  results  of  the  experiments  at  1.72  MPa  confinement 
corresponding  to  low,  medium,  and  high  striker  velocity.  In  all  three  cases,  after 
the  initial  stress  rise,  the  majority  of  the  material  deformation  took  place  under 
uniform  stress  conditions.  Under  each  condition,  mortar  exhibits  both  compressive 
and  dilatant  behavior.  Damage  was  not  observed  for  the  LVLP  case,  though  it  was 
present  in  the  other  two  cases.  Note  that  rate  effects  are  less  pronounced  than 
under  unconfined  conditions. 

Figures  3-42  to  3-44  are  representative  of  the  results  observed  at  moderate 
confining  conditions.  Once  again,  the  mortar  exhibits  compressibility  and  dilatancy. 
Rate  effects  again  do  not  appear  as  pronounced  as  in  the  unconfined  case.  The 
major  portion  of  each  experiment  also  took  place  under  uniform  stress  conditions. 
No  damage  was  observed  in  any  of  the  three  cases  presented. 

The  behavior  in  the  experiments  under  6.89  MPa  confinement  is  purely 
compressive,  hence  markedly  different  than  in  those  under  lower  confinement.  No 
cracking  was  observed.  The  uniform  stress  assumptions  again  appears  to  have 
applied  in  all  tests. 

In  order  to  quantify  the  influence  of  loading  rate  on  strength  at  the  various 
levels  of  confinement,  the  dynamic  strengths  given  in  Table  3-8  were  normalized 
by  the  quasi-static  strengths  corresponding  to  the  same  level  of  confinement,  thus 
yielding  the  DIF  values  as  a function  of  confinement.  In  Figure  3-63  the  DIF  is 
plotted  as  a function  of  the  log10(e).  It  appears  that  under  unconfined  conditions, 
mortar  follows  the  same  rate  sensitive  trends  as  those  reported  in  the  literature  for 
cementitious  materials  [40,  41,  42].  However,  under  confined  conditions  the  DIF 
values  for  mortar  are  significantly  lower. 

Several  conclusions  may  be  drawn  at  this  point.  Though  mortar  is  highly  rate 
sensitive  at  unconfined  conditions,  the  rate  sensitivity  appears  to  diminish  with 
confinement.  This  same  observation  was  also  reported  for  concrete  by  Malvern  and 
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Volume  Strain  (mm/mm) 


Figure  3-38:  Unconfined  mortar  response  at  three  different  strain  rates 


Figure  3-39:  Direct  compression  mortar  confined  (DCMC36)  - Lv,  Lp  condition 
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40:  Direct  compression  mortar  confined  (DCMC32)  - Mv,  Lp  condition 


Figure  3-41:  Direct  compression  mortar  confined  (DCMC28)  - Hv,  Lp  condition 


Ill 


Figure  3-42:  Direct  compression  mortar  confined  (DCMC15)  - Lv,  Mp  condition 


Strain  (mm/mm) 

Figure  3-43:  Direct  compression  mortar  confined  (DCMC20)  - Mv,  Mp  condition 
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Figure  3-44:  Direct  compression  mortar  confined  (DCMC18)  - Hv,  Mp  condition 

Jenkins  [34]  and  Malvern  et  al.  [43].  The  materials  susceptibility  to  damage  under 
shear  loads  appears  to  be  more  a function  of  confinement  than  loading  rate.  Also, 
though  the  uniform  stress  requirements  were  not  met  under  uniaxial  conditions, 
minimal  levels  of  confinement  yielded  valid  constitutive  property  data  over  the 
majority  of  the  duration  of  the  experiments. 
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45:  Direct  compression  mortar  confined  (DCMC29)  - L„,  Hp  condition 


Figure  3-46:  Direct  compression  mortar  confined  (DCMC35)  - Mv,  Hp  condition 
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Figure  3-47:  Direct  compression  mortar  confined  (DCMC26)  - Hv,  Hp  condition 


lo9io  e 


Figure  3-48:  Effect  of  confinement  on  mortar  DIF 
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WES500Q 

The  following  section  presents  the  dynamic  characterization  effort  undertaken 
for  WES5000.  The  same  analysis  techniques  utilized  for  the  mortar  samples  were 
employed.  A total  of  35  SHPB  experiments  were  conducted  under  both  confined 
and  unconfined  conditions.  Table  3-9  summarizes  the  conditions  under  which  the 
experiments  were  conducted.  The  naming  convention  utilized  is  the  same  as  for 
mortar,  obviously  DCCC  and  DCCU  stand  for  Dynamic  Compression  Concrete 
Confined  and  Unconfined,  respectively.  Results  as  well  as  comments  concerning  the 
experiments  are  summarized  in  Table  3 10.  The  same  damage  designations  letters 
used  in  discussing  the  results  on  mortar  are  utilized  for  WES5000. 

Figures  3-49  through  3-52  show  the  observed  damage.  Unlike  mortar,  the 
more  severe  type  “C”  response  was  observed  in  a number  of  cases  under  confined 
conditions.  In  order  to  evaluate  the  material  response  to  the  different  loading 
conditions,  the  principal  stress  difference  versus  axial,  radial,  and  volumetric 
strains  curves  were  plotted  . Representative  results  corresponding  to  each  loading 
condition  are  shown  in  Figures  3-53  through  3-62. 

Figure  3-53  shows  the  results  of  an  unconfined  SHPB  experiment  (DCCU03) 
conducted  under  the  low  striker  bar  velocity  conditions  (Lv).  The  lines  represent 
all  the  recorded  data,  while  the  symbols  (circles  and  triangles)  represent  data 
points  recorded  at  every  25  /j,s.  The  volumetric  stress  strain  curve  shows  that 
the  material  exhibits  both  compressive  and  dilatant  behavior.  Note  also  that 
WES5000  exhibits  strong  rate  dependency,  the  unconfined  strength  at  a strain 
rate  of  55.6  s-1  is  nearly  the  double  of  the  quasi-static  strength.  A damage  level  of 
“D”  was  observed  for  unconfined  conditions  regardless  of  the  striker  bar  velocity. 
Under  the  Davies-Hunter  criteria,  approximately  100  /is  would  again  be  required 
for  achievement  of  stress  uniformity.  Careful  examination  of  the  data  (Figure 
3-53)  shows  that  the  specimen  fails  in  a little  over  75  [.is,  thus  the  specimen  never 
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Table  3—9:  Test  conditions  for  WES5000  SHPB  experiments 


Experiment  # 

Design 

condition 

Striker(m) 

Pq(MPa) 

Pc(MPa) 

Date 

DCCC01 

LvMp 

1.0 

2.07 

3.45 

17-May-1999 

DCCC02 

LyMp 

1.0 

2.07 

3.45 

19-May- 1999 

DCCC03 

LvMp 

1.0 

2.07 

3.45 

19-May- 1999 

DCCC04 

MVMP 

1.0 

4.27 

3.45 

19- May- 1999 

DCCC05 

LVMP 

1.0 

2.07 

3.45 

25-May- 1999 

DCCC06 

MVMP 

1.0 

4.27 

3.45 

25-May- 1999 

DCCC07 

LyHp 

1.0 

2.07 

6.89 

28- May- 1999 

DCCC08 

LyHp 

1.0 

2.07 

6.89 

28-May-1999 

DCCC09 

H y Lp 

0.762 

3.45 

1.72 

24- Nov- 1999 

DCCC10 

Hy  Lp 

0.762 

3.45 

1.72 

30-Nov-1999 

DCCC11 

LyLp 

0.762 

1.38 

1.72 

6- Jan-2000 

DCCC12 

LyLp 

0.762 

1.38 

1.72 

13- Jan-2000 

DCCC13 

MyLp 

0.762 

2.41 

1.72 

17- Jan-2000 

DCCC14 

My  L p 

0.762 

2.41 

1.72 

18- Jan- 2000 

DCCC15 

HVHP 

0.762 

3.45 

3.45 

25-Jan-2000 

DCCC16 

HyMp 

0.762 

3.45 

3.45 

25- Jan-2000 

DCCC17 

MyHp 

0.762 

3.45 

6.89 

26- Jan- 2000 

DCCC18 

LyHp 

0.762 

1.38 

6.89 

26- Jan-2000 

DCCC19 

MyHp 

0.762 

2.41 

6.89 

27-Jan-2000 

DCCC20 

HyHp 

0.762 

3.45 

6.89 

31-Jan-2000 

DCCC21 

MyHp 

0.762 

2.41 

6.89 

31-Jan-2000 

DCCC22 

LyHp 

0.762 

1.38 

6.69 

14-Feb-2000 

DCCU01 

Ly 

0.762 

1.38 

- 

22-Dec- 1998 

DCCU02 

Ly 

0.762 

1.38 

- 

22-Dec- 1998 

DCCU03 

Ly 

0.762 

1.38 

- 

22-Dec- 1998 

DCCU04 

Ly 

0.762 

1.38 

- 

22- Jan- 1999 

DCCU05 

My 

0.762 

2.41 

- 

22- Jan- 1999 

DCCU06 

My 

0.762 

2.06 

- 

3-Feb-2000 

DCCU07 

My 

0.762 

2.06 

- 

3- Feb- 2000 

DCCU08 

My 

0.762 

2.06 

- 

3-Feb-2000 

DCCU09 

My 

0.762 

2.41 

- 

7-Feb-2000 

DCCU10 

My 

0.762 

2.41 

- 

7-Feb-2000 

DCCU11 

Hy 

0.762 

3.45 

- 

7-Feb-2000 

DCCU12 

Hy 

0.762 

3.45 

- 

14-Feb-2000 

DCCU13 

My 

0.762 

3.45 

- 

14-Feb-2000 
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Table  3-10:  Results  for  WES5000  SHPB  experiments 


Experiment  # 

V,  (m/s) 

£(S  X) 

Qmax  (MPa) 

Comments 

DCCC01 

— 

- 

100.18 

A 

DCCC02 

7.44 

48 

121.7 

A 

DCCC03 

7.52 

47 

121.1 

B 

DCCC04 

11.34 

54 

136.4 

B 

DCCC05 

7.27 

46.5 

128.6 

B 

DCCC06 

11.06 

77.8 

164 

C 

DCCC07 

7.26 

49.8 

85.1 

A 

DCCC08 

7.37 

50.0 

91.1 

A 

DCCC09 

- 

- 

- 

B (no  press,  signal) 

DCCC10 

- 

139.7 

166.2 

B 

DCCC11 

8.18 

57 

121.5 

B 

DCCC12 

8.18 

59 

123.57 

B 

DCCC13 

- 

104 

139 

C 

DCCC14 

11.51 

113.8 

138.8 

C 

DCCC15 

- 

113.9 

123.6 

B 

DCCC16 

- 

113.8 

123.6 

B 

DCCC17 

- 

0.762 

138.8 

B 

DCCC18 

8.19 

42 

138.1 

A 

DCCC19 

- 

80.8 

160.1 

B 

DCCC20 

14.27 

103.5 

176.8 

C 

DCCC21 

- 

77.4 

168.6 

B 

DCCC22 

8.01 

47.6 

123.6 

A 

DCCU01 

- 

63.3 

67.8 

D 

DCCU02 

- 

65.8 

89.2 

D 

DCCU03 

- 

55.6 

114.56 

D 

DCCU04 

- 

60.7 

89.64 

D 

DCCU05 

- 

109.3 

120.786 

D 

DCCU06 

10.37 

96.8 

121.978 

D 

DCCU07 

10.37 

101.3 

119.467 

D 

DCCU08 

10.38 

94.49 

146.455 

D 

DCCU09 

11.42 

105.37 

128.587 

D 

DCCU10 

11.39 

88.5 

153.008 

D 

DCCU11 

14.06 

129.15 

155.45 

D 

DCCU12 

13.98 

126.4 

163.769 

D 

DCCU13 

13.93 

133.39 

160.623 

D 
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Figure  3-49:  Damage  state  A - WES5000 


Figure  3-50:  Damage  state  B - WES5000 
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Figure  3-51:  Damage  state  C - WES5000 


Figure  3-52:  Damage  state  D - WES5000 
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achieves  a uniform  stress  state.  Thus,  only  the  peak  stress  as  a function  of  time 
may  be  considered  as  being  representative  of  the  material  behavior.  The  same 
conclusion  can  be  drawn  concerning  the  results  obtained  for  the  higher  impact 
velocities.  Also,  as  for  mortar  at  higher  velocities,  transition  to  dilatant  behavior  is 
not  observed. 

The  next  three  figures  (3-54  to  3-56)  present  representative  results  under  low 
confinement  corresponding  to  low,  medium,  and  high  loading  rates,  respectively. 

In  all  three  cases,  the  material  first  exhibits  compressive  behavior,  followed  by 
a transition  to  dilatancy.  Significant  dynamic  strength  enhancement  is  again 
observed,  though  not  on  the  same  order  as  in  the  unconfined  case.  Except  for  a 
small  portion  of  the  experiment  conducted  under  Lv  conditions,  the  specimens 
fail  prior  to  achieving  stress  uniformity  and  can  not  be  considered  valid  SHPB 
experiments  except  for  discerning  the  influence  of  loading  rate  on  failure  strength. 
For  all  three  conditions  presented,  severe  diagonal  cracking  was  observed. 

The  following  three  cases  concern  moderate  confinement  at  low,  medium,  and 
high  loading  rates.  Stress  uniformity  was  achieved  prior  to  failure.  Regardless  of 
the  strain  rate  , compressibility  followed  by  dilatant  behavior  is  observed.  Strain 
rate  influence  is  less  pronounced  than  under  unconfined  conditions.  No  visible 
damage  was  observed  for  the  lowest  strain  rate  conditions,  however  significant 
diagonal  cracking  was  observed  under  higher  strain  rates. 

Finally,  Figures  3-60  to  3-62  show  the  test  results  under  6.89  MPa  conditions 
corresponding  to  three  different  loading  rates.  Under  moderate  and  high  loading 
rate  conditions  the  material  exhibits  both  compressive  and  dilatant  behavior,  but 
at  the  lowest  loading  rate  tested,  only  compressive  behavior  was  observed.  No 
damage  was  observed  at  the  lowest  loading  rate,  while  significant  cracking  occurred 
under  higher  loading  rate  conditions.  In  order  to  quantify  the  influence  of  loading 
rate  on  strength,  the  dynamic  strengths  given  in  Table  3-10  were  normalized  by  the 


q (MPa)  „ q (MPa) 
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-53:  Direct  compression  concrete  unconfined  (DCCU03)  - Lv  condition 


Figure  3-54:  Direct  compression  concrete  confined  (DCCC11)  - Lv,  Lp  condition 
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Strain  (mm/mm) 


55:  Direct  compression  concrete  confined  (DCCC13)  - Mv,  Lp  condition 
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Figure  3-56:  Direct  compression  concrete  confined  (DCCC10)  - Hv , Lv  condition 
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57:  Direct  compression  concrete  confined  (DCCC02)  - Lv , Mp  condition 


Figure  3-58:  Direct  compression  concrete  confined  (DCCC04)  - Mv,  Mp  condition 
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Figure  3-59:  Direct  compression  concrete  confined  (DCCC15)  - Hv,  Mp  condition 

quasi-static  strengths  corresponding  to  the  same  level  of  confinement,  thus  yielding 
the  DIF  values  as  a function  of  confinement.  In  Figure  3-63  the  DIF  is  plotted  as 
a function  of  the  log10(e).  It  can  be  concluded  that  under  unconfined  conditions, 
WES5000  follow  the  same  trends  as  mortar  and  those  reported  in  the  literature 
[40,  41,  42]  for  other  cementitious  materials.  However,  under  confined  conditions 
the  DIF  values  for  WES5000  are  significantly  lower.  This  is  again  in  agreement 
with  the  findings  of  Malvern  and  Jenkins  [34]  and  Malvern  et  al.  [43].  Note 
also  that  WES 5000  generally  exhibited  far  more  cracking  under  similar  loading 
conditions  than  mortar.  Thus,  while  WES5000  is  a significantly  stronger  material, 
it  also  exhibits  far  less  ductility.  Finally,  it  was  again  observed  that  minimal  levels 
of  confinement  can  be  utilized  to  generate  valid  constitutive  property  data  in  a 
SHPB  even  for  brittle  materials  such  as  concrete. 
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60:  Direct  compression  concrete  confined  (DCCC18)  - Lv , Hp  condition 


Figure  3-61:  Direct  compression  concrete  confined  (DCCC21)  - Mv,  Hp  condition 
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Figure  3-62:  Direct  compression  concrete  confined  (DCCC17)  - Hv,  Hp  condition 
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Figure  3-63:  Effect  of  confinement  on  WES5000  DIF 


CHAPTER  4 

CONSTITUTIVE  MODELING 
4.1  Introduction 

A constitutive  model  consists  of  a set  of  stress-strain  (time)  relationships.  The 
choice  of  a material  model  formulation  is  generally  guided  by  three  factors,  these 
being: 

• The  scale  of  the  problem 

• The  type  of  material  and  behavior  of  interest  for  a given  problem 

• Available  computational  resources. 

In  this  chapter,  after  a review  of  concrete  continuum  based  constitutive 
models,  an  elastic/ viscoplastic  model  for  WES5000  concrete  is  presented.  The 
model  developed  is  based  on  data  obtained  in  the  series  of  tests  reported  in 
Chapters  2 and  3.  The  procedure  for  determination  of  the  constitutive  functions 
and  the  parameters  involved  in  their  expressions  is  outlined.  Finally,  comparisons 
are  made  between  the  theoretical  predictions  and  experimental  results. 

4.2  Concrete  Constitutive  Modeling  Review 
4.2.1  Introductory  Note 

The  development  of  material  models  is  inseparably  linked  to  the  development 
of  advanced  numerical  analysis  techniques  and  computational  capabilities.  The  ac- 
curacy and  utility  of  numerical  simulations  are  limited  by  the  material  descriptions 
(Zukas  et  al.  [1])  whereas  the  degree  of  complexity  of  constitutive  models  is  limited 
by  the  computation  time  and  robustness  of  numerical  algorithms  used  in  the  sim- 
ulations. Nowadays,  general  purpose  finite  element  and  finite  difference  methods 
provide  means  for  simulating  penetration  events  and  allow  a wide  range  of  thought 
experiments  that  would  not  be  otherwise  possible.  One  of  the  earliest  works  in  the 
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area  of  Finite  Element  Methodology  was  that  of  mathematician  Richard  Courant, 
who  in  1943  used  piecewise  linear  functions  on  triangles  [44]to  obtain  numerical 
solution  of  partial  differential  equations.  The  triangular  linear  element  is  referred  to 
as  the  Courant  element.  The  work  of  Turner,  Clough,  Martin  and  Topp  in  1956  is 
generally  regarded  as  the  first  paper  on  Finite  Element  Methodology  in  engineering 
[45].  In  this  landmark  paper,  the  authors  proposed  to  use  the  displacement  method 
to  derive  the  stiffness  of  plane  elements.  In  1960,  Clough  [46]  coined  the  phrase  Fi- 
nite Element  Method  (FEM)  in  a paper  presented  at  the  2nd  ASCE  Conference  on 
Electronic  Computation.  The  first  application  of  FEM  to  the  analysis  of  reinforced 
concrete  beams  was  by  Ngo  and  Scordelis  [47]  in  1967.  Due  to  the  lack  of  generally 
accepted  constitutive  equations  for  concrete  [48],  in  most  of  the  early  FEM  anal- 
yses concrete  was  considered  to  obey  a linear  elastic  stress-strain  law.  Also,  FEM 
methods  were  mostly  applied  to  simulate  quasi-static  structural  problems  while 
finite  difference  methods  were  applied  to  hypervelocity  impact  problems.  However, 
early  finite  difference  formulations  were  used  in  conjunction  with  material  models 
that  did  not  account  for  strength  effects.  All  materials  were  treated  as  fluids,  thus 
the  term  “hydrodynamic  code,”  or  “hydrocode,”  came  into  use.  To  better  describe 
lower  velocity  regime  problems,  later  versions  of  hydrocodes  included  strength  ef- 
fects. More  details  concerning  hydrocode’s  development  can  be  found  in  the  review 
paper  by  Anderson  [49] 

4.2.2  Micromechanical  Modeling 

There  is  a tremendous  amount  of  information  available  in  the  open  literature 
concerning  modeling  the  behavior  of  concrete.  In  micromechanics  based  effective 
medium  models,  the  overall  continuum  properties  of  concrete  are  estimated  using 
different  averaging  (homogenization)  procedures  from  the  properties  and  structure 
of  each  of  its  microconstituents.  Specifically,  concrete  is  considered  to  be  a com- 
posite consisting  of  randomly  distributed  inclusions  (aggregate  particles)embedded 
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in  a continuous,  porous  matrix  (hardened  cement  paste) . It  is  to  be  noted  that 
the  most  widely  used  effective  medium  homogenization  methods  such  as  the  self- 
consistent  method,  the  differential  scheme,  the  generalized  self-consistent  method 
and  Mori-Tanaka  methods  [50]  do  not  apply  to  arbitrarily  located  and  oriented  mi- 
crocracks. Embedded  in  the  mentioned  averaging  techniques  is  the  restriction  that 
there  is  a dilute  concentration  of  defects  and  that  those  defects  are  non-interacting 
[51].  Improved  formulations  that  account  for  interacting  microcracks  have  been 
developed  by  Ju  and  co-workers  (e.g.  Ju  and  Tseng  [52])  for  2-D  conditions.  The 
great  majority  of  these  micromechanics  based  models  are  applicable  only  to  short- 
term quasi-static  stress  conditions.  Very  few  models  account  for  high-strain  rate, 
high-pressure  behavior.  Examples  include  Xia  [53]  and  Thoma  et  al.  [54],  An  al- 
ternative method  to  obtain  the  macro-level  properties  in  terms  of  the  properties  of 
the  constituents  is  offered  by  the  mixture  theory  [55].  Such  an  approach  was  used 
by  Ortiz  [56]  for  concrete.  In  all  effective  medium  models,  the  material  is  treated 
as  quasi-homogeneous  (statistical  homogeneity  assumptions)at  the  macro-scale.  At 
difference  with  effective  medium  models,  in  Discrete  Element  models  (DEM)  the 
material  is  treated  as  heterogeneous  at  the  macro-scale.  Specifically,  the  material 
is  discretized  into  individual  elements  or  blocks,  which  are  allowed  to  interact  as 
stress  is  applied  [57].  Void  growth  and  cracking  are  represented  by  separation  of 
adjoining  elements  due  to  shear  induced  sliding.  Hentz  et  al.  [58]  has  recently 
reported  on  a concrete  DEM  model  in  which  elastic  forces  with  a local  rupture 
criteria  are  applied  between  interacting  elements.  Macro-scale  level  experiments 
and  numerical  simulations  are  used  to  calibrate  the  model  parameters.  Applica- 
tions of  DEM  to  modeling  the  dynamic  behavior  of  concrete  have  been  reported  by 
Camborde  et  al.  [59]  and  Donze  et  al.  [60],  who  used  the  technique  in  the  analysis 
of  flyer  plate  and  SHPB  experiments,  respectively.  Very  good  agreement  with 
experiments  was  obtained.  However,  DEM  approach  is  generally  not  applicable 
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to  large  scale  problems.  The  microstructural  variables  that  must  be  recorded  and 
used  in  the  computation  of  the  deformations  may  require  an  unduly  large  amount 
of  storage  capacity  and  manipulation  of  data.  Computation  time  required  to  solve 
even  simple  problems  can  be  excessive.  Material  boundary  interface  treatments, 
numerical  accuracy  and  mesh  dependencies  may  also  dominate  the  problem. 

4.2.3  Macromechanical  Models 

A promising  approach  for  practical  large-scale  problems  is  to  treat  concrete  as 
a homogeneous  continuum.  Two  main  methods  have  been  adopted:  (1)  empirical 
methods  and  (2)  constitutive  modeling.  Empirical  models  result  from  curve-fits  of 
experiments  by  means  of  simple  functions.  The  characteristics  of  these  functions 
are  chosen  just  to  meet  the  shape  of  the  experimentally  observed  stress-strain 
curves  and  do  not  have  any  physical  basis.  Generally,  empirical  laws  are  one- 
dimensional and  apply  only  to  problems  having  boundary  conditions  and  loading 
paths  similar  to  those  used  in  the  tests  based  on  which  the  laws  were  derived.  It  is 
possible  to  take  into  account  the  macroscopic  experimental  results  by  constructing 
a continuum  model  with  a set  of  thermodynamic  internal  state  variables  that 
characterize  the  microstructure  of  the  material.  In  the  case  of  cementitious 
materials,  damage  related  parameters  are  used  extensively  as  internal  variables 
accounting  for  the  effect  of  microcracking  on  the  macroscopic  behavior.  Most  of 
the  early  work  and  most  of  the  models  that  have  been  actually  applied  to  concrete 
are  based  on  a single  scalar  representation  of  the  damage  (e.g.  Kachanov  [61], 
Lemaitre  and  Chaboche[62],  Lemaitre[63],  etc).  However,  test  results  under  both 
quasi-static  and  dynamic  conditions  show  that  damage  itself  and  its  effect  on  the 
material  effective  compliance  are  anisotropic.  To  account  for  the  directionality  of 
damage,  different  higher  order  tensorial  representations  were  proposed.  The  list 
runs  from  vectors  (e.g.  Davison  and  Stevens  [64],  Krajcinovic  and  Fonseka  [65], 
Talreja  [66]),  second  order  tensors  (e.g.  Vakulenko  and  Kachanov  [67],  Dragon 
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and  Mroz  [68],  Kachanov  [69,  70],  Cordebois  and  Sidoroff  [71],  Murakami  [72], 
Murakami  and  Kamiya  [73],  etc.)  to  fourth-order  tensors  (e.g.  Chaboche  [74], 

Simo  and  Ju  [75],  Chow  and  Wang  [76],  Lubarda  and  Krajcinovic  [77],  Halm 
and  Dragon  [78],  etc.).  However,  the  existing  continuum  damage  mechanics 
(CDM)  models  are  rate-independent  and  unable  to  describe  simultaneously: 

(1)  the  directionality  of  damage  and  the  corresponding  induced  anisotropy  in 
behavior  and  (2)  damage  activation-deactivation  process.  Very  few  CDM  models 
accounting  for  the  dynamic  behavior  of  cementitious  materials  have  been  proposed. 
At  Sandia,  Taylor  et  al.  [79],  have  proposed  a CDM  model  which  uses  a scalar 
damage  parameter  and  combines  the  theory  of  fracture  mechanics  governing  the 
behavior  of  individual  cracks  with  the  fragment  size  equation,  which  is  based  on 
kinetic  energy  considerations  (Grady  [80]).  However,  the  model  allows  for  damage 
growth  only  under  hydrostatic  tension;  under  compression  loading  the  model  is 
rate  independent.  Considerable  efforts  have  also  been  undertaken  to  adapt  models 
developed  for  metallic  materials  to  describe  the  non-linear  irreversible  behavior  of 
concrete.  The  most  widely  used  formulations  are  Mohr-Coulomb  or  Drueker-Prager 
[48]  which  account  for  the  experimentally  observed  influence  of  hydrostatic  pressure 
on  the  behavior.  It  is  thus  assumed  that  concrete  is  perfectly  plastic,  i.e.  the  yield 
locus  is  fixed  in  the  stress  space.  Hence,  the  yield  surface  coincides  with  the  failure 
surface.  According  to  Mohr-Coulomb  theory  failure  occurs  under  compression  when 

r = an  tan  0O  + C0  (4.1) 

where:  r is  the  shear  stress  on  the  failure  plane  , on  is  the  normal  stress  on  that 
plane,  and  4>0  and  Co  are  material  constants.  The  angle  of  internal  friction,  0O , and 
the  cohesion  Co  can  be  determined  from  triaxial  compression  data.  Under  tensile 
conditions  it  is  generally  assumed  that  the  material  fails  if  a predetermined  level  is 
achieved.  This  so  called  ’’tension  cut-off”  is  estimated  from  unconfined  tests.  The 
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Drucker-Prager  criterion  is  given  by: 

Fs  = \/~Ji  ~ P tan  (3  - d (4.2) 

where  J 2 is  the  second  invariant  of  the  stress  deviator,  p is  the  mean  stress,  (3 
is  a measure  of  the  internal  angle  of  friction,  and  d is  the  cohesion.  Note  that 
Mohr-Coulomb  and  Drucker-Prager  surfaces  are  open  in  the  direction  of  increasing 
hydrostatic  compression.  If  the  associativity  assumption  is  adopted,  i.e.  the  plastic 
potential  coincides  with  the  yield  function,  these  criteria  predict  levels  of  dilatancy 
far  in  excess  of  what  it  is  observed  experimentally.  The  reason  for  this  can  be  seen 
graphically  in  Figure  4-2.  Indeed,  the  incremental  volumetric  strain  rate  vector  is 
always  in  the  direction  of  volumetric  expansion.  The  Osborn  model  and  Osborn 
and  Matuska  model  for  concrete  which  is  implemented  in  HULL  [81],  a widely  used 
and  accepted  continuum  mechanics  code,  employs  a Mohr-Coulomb  yield  criterion. 
The  material  parameters  are  obtained  based  on  triaxial  data  generated  by  Chinn 
and  Zimmerman  [6]  in  the  1960s.  It  is  stated  that  the  Osborn  model  has  been  used 
“successfully  in  scores  of  calculations  of  penetrations  into  concrete  at  velocities 
from  a few  hundred  to  over  1000  feet/second.”  However,  experimental  evidence 
clearly  suggests  that  for  cementitious  materials  a clear  distinction  exists  between 
the  stress  conditions  at  which  inelastic  behavior  initiates  and  those  corresponding 
to  failure.  Among  the  models  that  account  for  pre-failure  irreversible  response 
double-surface  models  have  become  quite  popular.  These  models  are  characterized 
by  a hardening  yield  locus,  referred  to  as  cap,  and  by  a fixed  failure  surface  which 
is  a limit  of  the  admissible  stress  states.  Typically,  cap  models  are  completed 
by  providing  a non-associative  flow  rule  and  elastic  moduli  which  are  stress-level 
dependent.  Figure  4-3  illustrates  the  modified  Drucker/Prager  cap  constitutive 
model.  The  shear  failure  surface  is  Drucker-Prager  surface  given  by  Equation  4.2. 
The  cap  has  an  elliptical  shape  in  the  plane  ( , I\ ) and  hardens  or  softens  as 


a function  of  the  volumetric  plastic  strain:  compaction  produces  hardening,  while 
dilation  produces  shear  induced  softening.  It  is  defined  by: 


Fc  = < (p  ~ Pa)2  + 


RVT2 


-i  2' 


1/2 


(1  + a — cr/cos  P) 


R(d  + pa  tan/?)  = 0 


(4.3) 


where  R,  a are  constants  and  pa  is  a function  of  the  plastic  volumetric  strain. 
Plastic  flow  is  associated  on  the  cap,  and  nonassociated  on  the  failure  surface. 

Among  the  authors  that  applied  this  cap  type  model  to  describe  concrete 
behavior  under  quasi-static  conditions  are  Faruque  and  Chang  [82],  Simo  et  al. 

[83],  etc.  Very  recently,  Warren  et  al.  [84]  applied  the  same  type  cap  model  to 
describe  the  response  of  concrete  targets  to  high-velocity  impact.  Schwer  and 
Murray  [85]  proposed  a cap  model  that  depends  on  all  three  stress  invariants.  The 
numeric  singularities  at  the  cap  and  failure  surface  intersections  are  handled  by 
smoothing  the  intersections.  This  model  has  been  implemented  into  the  FEM 
codes  DYNA2D  and  DYNA3D.  While  cap  type  models  have  become  widely 
accepted  because  they  pair  a relative  simplicity,  and  therefore  a limited  number 
of  parameters,  to  a description  of  the  main  features  of  the  behavior  of  concrete, 
their  applicability  to  high-velocity  problems  is  questionable.  This  is  due  to  the  fact 
that  such  model  are  rate-insensitive.  A review  of  rate-type  constitutive  equations  is 
presented  in  the  next  section. 

Rate  Dependent  Models 

To  characterize  the  temporal  growth  of  permanent  deformation  (flow  of 
matter  by  creep)  and  strain  rate  effects  on  the  behavior  of  certain  materials 
viscoplastic  models  have  been  proposed.  The  basis  of  what  is  today  called  the 
mathematical  theory  of  viscoplasticity  was  laid  by  Bingham  [86].  According  to 
Bingham,  viscoplasticity  is  made  up  of  two  fundamental  properties  which  are  ’’yield 
value”  and  ” mobility” . The  yield  value  is  dependent  on  the  shear  stress  required  to 
start  the  deformation  while  the  mobility  is  proportional  to  the  rate  of  deformation 
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Figure  4-1:  Mohr-Coulomb  yield  criterion 


after  the  yield  value  has  been  exceeded.  Hohenemser  and  Prager  [87]  proposed  the 


first  elastic/ viscoplastic  model.  The  basic  hypotheses  of  this  model  are  that  the 
total  strain  rate  is  the  sum  of  two  components:  the  elastic  component  ef3  and  the 
viscoplastic  component  such  that 


where  k is  the  yield  value  in  shear  and  a'l}  is  the  stress  deviator. 

Concerning  the  development  of  the  mathematical  modeling  of  strain  rate 
effects  on  the  irreversible  strain  properties  of  materials,  an  important  date  is  1948, 
with  the  introduction  of  Sokolovskii’s  [88]  rate- type  constitutive  equation  of  the 
general  form: 


(4.4) 


and  the  viscoplastic  strain  rate  being: 


(4.5) 


& = 0(e,  a )e  + ip(e,  a) 


(4.6) 
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Figure  4-2:  Dilatancy  associated  with  yield  surface  open  in  positive  direction 


with  0(e,  a)  = E,  and 


ip(e,a) 


0 if  |(j|  < <Jy 

(—sign  a)F(\a\  — aY)  if  \a\  > aY 


(4.7) 


where  ay  > 0 is  the  yield  limit  and  F is  a smooth  function  with  F(r ) > 0, 

F'(r ) > 0 for  r > 0 and  F( 0)  = 0.  For  any  state  (e.a)  with  \a\  < aY,  the 

response  is  elastic.  From  any  state  (e,a)  with  |cr|  > aY,  the  material  relaxes  to  the 

state  (e,aY)  if  a > aY  or  to  the  state  (e,  -aY)  if  a < -aY.  Malvern  [89,  90]  also 

proposed  a rate-type  law  of  the  type  given  by  Equation  4.6  with  </>(e,  a)  = E and 


il)(e,a)  = 


0 if  0 < a < f(e),  e > 0 

-kF(a  - f(e))  if  a > /(e),  e > 0 


(4.8) 
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Figure  4-3:  Modified  Drucker/Prager  cap  model 


where  k > 0 is  a viscosity  constant.  Malvern  considered  F(r)  = r.  Other  particular 
forms  of  the  function  F can  be  found  in  the  monograph  by  Cristescu  and  Suliciu 
[88].  The  curve  a = /(e)  is  the  quasi-static  loading  curve. 

Another  important  date  is  1951,  with  Malvern’s  law  of  logarithmic  dependence 
of  strength  on  strain  rate  [91]: 


where  a — /(e)  is  the  static  stress-strain  relation,  e is  the  strain  rate  (irreversible 
strain  rate),  whereas  a and  b are  material  constants. 

The  Johnson-Cook  model,  which  is  one  of  the  most  widely  used  models  for 
describing  high  strain  behavior  in  metals,  is  essentially  a variation  on  Equation  4.9. 


a — 


/(e)  + aln  (1  + be) 


(4.9) 


(4.10) 


In  this  equation,  £q  is  a reference  strain  rate  equal  to  1.0s  x,  T*  is  the  homologous 
temperature,  and  A,  B,  n,  m and  C are  material  constants.  The  Johnson-Cook 
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model  was  modified  for  use  on  brittle  and  geologic  materials  by  Holmquist.  Johnson 
and  Cook  [92].  This  model  states  that: 

<7*=  [A{1- D)  + BP*n]  [1  + CTne*]  (4.11) 

where:  a*  is  the  normalized  equivalent  stress  defined  as  cr*  = a//',  D is  a damage 
parameter,  P is  the  hydrostatic  pressure,  P*  = P/ f'c  is  the  normalized  pressure, 
and  A,  B,  N and  C are  material  constants.  A depiction  of  the  model  is  presented 
in  Figure  4-4.  The  tensile  portion  of  the  model  is  defined  by  T*  = T / /',  where 
T*  is  the  normalized  maximum  tensile  hydrostatic  pressure  and  T is  the  maximum 
hydrostatic  pressure  the  material  can  withstand.  The  material  constant  SMAX 
is  the  normalized  maximum  strength  that  can  be  developed.  SMAX  essentially 
defines  a von-Mises  yield  surface.  The  HJC  model  and  similar  formulations  have 


Figure  4-4:  Holmquist- Johnson-Cook  Model 

been  implemented  into  a number  of  continuum  mechanics  codes  including  EPIC 
and  CTH  and  utilized  successfully  for  a broad  range  of  penetration  calculations 
[93,  94,  95,  96,  97].  As  can  be  seen  in  Figure  4-4  the  model  does  not  account  for 
increased  material  strength  at  elevated  pressures. 
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The  early  1960’s  mark  the  development  of  the  first  3-D  elastic/ viscoplastie 
equations  by  Perzyna  [98]  and  Cristescu  [99].  Following  Perzyna  [98],  the  structure 
and  basic  assumptions  involved  in  the  formulation  of  an  elastic/ viscoplastic  type 
constitutive  equation  is  presented.  The  formulation  is  similar  to  that  of  plasticity. 
Decoupling  between  elastic  and  viscoplastic  effects  is  assumed  such  that: 

eij=£ijE  + £^  (4.12) 

where  efj  is  the  elastic  component  and  is  the  permanent  (irreversible)  part  of 
the  strain.  The  existence  of  a yield  function  which  depends  on  the  stress  state  and 
on  the  accumulated  viscoplastic  deformation  is  postulated: 

= — p^-1  (4.13) 

where  / and  k are  scalar  valued  functions.  The  yield  value,  k , depends  on  the 
history  of  strain.  The  elasticity  domain  is  defined  as:  F < 0 or  / < k.  The 
problem  of  activation  of  viscoplastic  flow  is  much  simpler  than  in  plasticity  since 
any  state  outside  of  the  elastic  domain  is  viscoplastic  and  the  viscoplastic  strain 
increment  depends  only  on  the  current  values  of  the  stress  and  hardening  and  not 
on  the  stress  increment.  No  consistency  condition  needs  to  be  reinforced.  In  the 
case  of  associated  viscoplasticity,  the  flow  rule  is  thus  given  by: 

e”p  = M(*(F)>|£  (4.14) 

where  /i  is  the  fluidity  or  viscosity  parameter,  <E>  is  a monotonic  function  and  the 
symbols,  <>,  represent  the  so  called  Macaulay  bracket,  used  to  denote  the  positive 
part  of  a function  (i.e.  (A)  = (1/2)  ( A + |j4|)). 

In  the  case  of  non- associated  viscoplasticity  it  is  assumed  that  the  potential 
from  which  the  viscoplastic  strain  rate  is  derived  does  not  coincide  with  the  yield 
function.  Generally,  the  modeling  strategy  used  is  to  consider  that  the  expressions 
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of  the  yield  and  viscoplastic  potentials  are  the  same  as  the  quasi-static  (inviscid) 
expressions.  For  example,  in  Simo  et  al.  [83]  associated  flow  viscoplastic  model  for 
concrete,  the  expressions  of  the  constitutive  functions  coincide  with  the  expressions 
of  the  inviscid  cap-model.  This  model  is  reportedly  easy  to  implement  and  provides 
more  accurate  results  than  the  inviscid  form.  The  associated  flow  assumption 
is  however  responsible  for  non-representative  compressive  and  dilatant  behavior 
predictions. 

In  contrast  to  this  approach,  Cristescu  [7,  8]  proposed  to  determine  for  a 
given  material  the  specific  expressions  of  the  constitutive  functions  based  on 
data  corresponding  to  several  strain  rates.  The  general  form  of  the  viscoplastic 
constitutive  formulation  proposed  by  Cristescu  is: 


_ , f 1 

ij  2 G \3K 


2^7 ) ptiij  + kT 


A w{t) 

\ H(amn) 


dF 

doij 


(4.15) 


where: 


W = Irreversible  Stress  Work 
F = Viscoplastic  Potential 
H = Yield  Function 
kT  = Viscosity  Coefficient 
p = mean  normal  stress 


I<  = Bulk  Modulus 
G = Shear  Modulus 
£ij  = Strain  Tensor 
Oij  = Stress  Tensor 
5ij  = Kronecker  Delta 


{A)  = i(A+\A\). 

No  a priori  assumptions  concerning  the  specific  mathematical  expressions  of 
the  yield  function  and  viscoplastic  potential  are  made.  The  formulation  does  not 
assume  a decoupling  between  shear  and  volumetric  effects  and  does  not  require 
a separate  equation  of  state  in  order  to  handle  the  pressure  volume  relationship. 
Since  geomaterials  and  cementitious  materials  may  display  both  compressible  and 
dilatant  behavior,  the  concept  of  compressibility/dilatancy  boundary  is  integrated 
in  the  Cristescu  model.  This  allows  one  to  describe  both  compressibility  and 
dilatancy  using  a unique  continuous  yield  function  and  viscoplastic  potential, 
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respectively.  This  is  a significant  advantage  over  rate  sensitive  cap  models  (e.g. 
Simo  et  al.)  since  it  captures  in  a continuous  way  the  transition  between  the 
two  regimes  of  behavior.  The  model  has  been  implemented  for  a wide  range 
of  materials  [7,  100,  101,  12].  The  formulation  also  does  not  require  a separate 
equation  of  state  in  order  to  handle  the  pressure  volume  relationship.  It  has 
however,  not  been  utilized  for  standard  high-strength  concrete  or  applied  to 
highly  dynamic  loading  environments.  Experimental  data  on  some  geomaterials 
and  granular  materials  do  not  support  the  hypothesis  concerning  the  existence 
of  a viscoplastic  potential.  An  alternative  flow  rule  formulation  was  proposed 
by  Cazacu  et  al.  [9]  using  representation  theorems  for  isotropic  scalar  valued 
tensor  functions.  The  remainder  of  this  chapter  deals  with  the  development  of 
a an  elastic/viscoplastic  model  for  WES5000  concrete.  The  elastic- viscoplastic 
formulation  of  Cristescu  is  chosen  due  to  its  apparent  capability  to  capture 
the  behavior  of  interest  and  the  fact  that  there  are  no  a priori  limitations  or 
restrictions  regarding  the  yield  functions.  To  better  describe  the  peculiarities  of  the 
viscoplastic  flow  behavior  of  WES5000,  the  approach  proposed  by  Cazacu  et  al.  [9] 
was  adopted. 

4.3  Elastic-Viscoplastic  Model  Development 
The  general  form  of  the  elastic-viscoplastic  constitutive  equation  proposed  in 
the  previous  section  is  repeated  below  4.15: 


■ _ 

13  ~ 2G 


3K~2G)p5ij  + kT 


W{t)  \ dF 

F {(7rnn')  ! dOi j 


(4.15) 


If  material  displacements  and  rotations  are  assumed  to  be  small,  an  additive 
decomposition  of  the  rate  of  deformation  into  an  elastic  instantaneous  part  and  an 
irreversible  part  can  be  made. 


(4.16) 
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The  reversible  deformation  rate  is  given  by  the  generalized  Hooke’s  law: 


(4.17) 


in  which  both  K and  G are  variable.  The  irreversible  deformation  rate  is: 


Recalling  the  definition  of  the  Macaulay  bracket  from  the  preceding  section,  if  the 
stress  state  is  such  that  H(am^n)  < W(t),  then  the  material  response  is  elastic.  If 
H (o"m,n)  > IT(t),  further  irreversible  deformation  may  occur.  In  this  model,  no 
a priori  assumptions  concerning  the  mathematical  expressions  of  the  viscoplastic 
potential,  F,  and  yield  function,  H are  made.  These  expressions  are  determined 
from  data  collected  on  a given  geologic/cementitious  material  in  Karman  type 
experiments  or  true  triaxial  tests. 

The  first  step  in  the  development  of  an  elastic-viscoplastic  model  for  a par- 
ticular material  is  the  evaluation  of  the  elastic  parameters  (i.e.  K and  G ).  The 
accurate  determination  of  these  moduli  is  critical.  The  importance  for  evaluation 
of  the  reversible  deformation  rate  (Equation  4.17)  is  obvious.  In  addition,  since 
E1  = i — ee,  & correct  estimation  is  crucial  in  the  determination  of  every  aspect 
of  the  model.  These  parameters  can  be  determined  based  on  triaxial  or  Karman 
type  experiments.  The  elastic  parameters  are  evaluated  from  the  slopes  of  unload- 
reload  cycles  of  such  an  experiment.  The  bulk  modulus,  K , is  the  slope  of  the 
unload-reload  portion  of  the  volume  strain  vs.  mean  normal  stress  curve  of  either 
the  hydrostatic  or  deviatoric  part  of  the  experiment.  The  Young’s  modulus,  E , is 
taken  from  the  slope  of  the  axial  stress  vs.  strain  unload- reload  cycle.  The  shear 


(4.18) 


where  the  irreversible  stress  work,  W(t),  is  calculated  as: 


(4.19) 
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modulus,  G , may  be  calculated  from  the  relationship: 


G = 


3KE 
9 K - E 


(4.20) 


Thus,  all  the  elastic  parameters  of  the  model  can  be  derived  from  a series  of 
carefully  controlled  hydrostatic  and  deviatoric  creep  tests. 

4.4  Elastic  Parameter  Evaluation 

As  was  mentioned  in  Section  4.2  (Elastic-Viscoplastic  Modeling),  the  Young’s 
modulus,  E,  is  determined  from  the  slope  of  the  axial  strain  vs.  principal  stress 
difference  curve  unload-reload  cycles.  Values  for  E from  the  different  quasi-static 
experiments  were  previously  presented  in  Section  2.15  as  a function  of  mean 
normal  stress  in  Figure  2-57.  Visual  inspection  of  the  figure  indicates  a very  clear 
and  consistent  trend  in  the  data,  with  all  the  experiments  in  general  agreement. 
According  to  elastic  wave  theory,  the  maximum  value  of  E is  given  by: 


E = pv% 


(4.21) 


where  p is  the  material  density  and  vb  is  the  axial  bar  velocity,  or  P-wave  velocity, 
using  the  notation  of  Table  2-3.  Using  values  of  density  and  wave  speed  tabulated 
for  WES5000  in  Table  2-3,  the  maximum  value  for  E should  be  on  the  order  of 
61000  MPa.  The  data  in  Figure  2-57  does  indeed  appear  to  be  asymptotically 
approaching  such  a value.  A decaying  exponential  function  of  the  form: 

E(p)  = a — be~',  with  a,b,c>  0 (4.22) 


was  found  to  approximate  the  data  quite  well.  Here,  p is  the  mean  normal  stress, 
or  pressure,  which  is  defined  as: 


(4.23) 


while  a represents  the  asymptotic  limit  towards  which  E tends  as  p 
P = 0, 


oc.  For 
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£'(0)  = a — b, 


(4.24) 


so,  b — E( oo)  — E( 0).  Taking  the  derivative  of  Equation  4.22  with  respect  to  p 
yields: 


Thus,  c is  directly  related  to  the  rate  of  increase  of  E with  respect  to  p.  Values  for 
the  parameters  were  found  to  be:  a = 61000  MPa,  b — 18920  MPa,  and  c = 329.065 
MPa.  Comparison  between  theoretical  and  experimental  variation  of  E with  the 
mean  normal  stress,  p,  is  presented  in  Figure  4-5,  where  the  solid  line  represents 
the  theoretical  model  results.  The  correlation  between  theoretical  and  experimental 
values  for  E is  excellent. 

Values  for  the  bulk  modulus  were  taken  from  the  slope  of  the  unload-reload 
cycles  of  the  volume  strain  versus  mean  normal  stress  curves,  for  both  the  hydro- 
static and  deviatoric  phases  of  the  Karman  experiments  as  discussed  in  Section 
2.14.  Results  of  the  data  analysis  for  K were  presented  in  Figure  2-58.  Though  the 
experimentally  determined  values  for  K are  in  good  agreement  (with  the  exception 
of  the  first  experiment  conducted  at  300  MPa  confining  pressure),  due  to  the  un- 
certainty in  estimating  K,  it  was  decided  to  take  a mean  value  of  v = 0.175  for  the 
Poisson  ratio  and  calculate  K from  the  elasticity  relationship: 


(4.25) 


For  p — 0: 


E'(  0)  = -. 

c 


(4.26) 


E 


3(1  - 2u) 


K = 


(4.27) 
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Figure  4-5:  Young’s  modulus  ( E ) as  a function  of  mean  normal  stress 


Several  other  factors  entered  into  the  making  of  this  decision.  For  the  sake  of 
simplicity  in  this  first  implementation  of  the  model,  it  was  seen  as  advisable  to 
limit  the  number  of  coefficients  involved  in  the  model.  Also,  due  to  the  overall 
higher  level  of  confidence  in  the  axial  measurements,  if  possible  it  would  be  best 
to  have  E as  the  sole  pressure  dependent  elastic  function.  The  value  selected  for 
the  mean  Poisson’s  ratio  is  well  within  the  range  of  values  typically  observed  for 
concrete  (see  Neville  [15]).  Figure  4-6  presents  the  experimentally  determined 
values  for  K , along  with  the  function  defined  by  Equation  4.27.  The  assumption 
of  a constant  Poisson’s  ratio  was  quite  reasonable  and  resulted  in  a good  match 
between  the  theoretical  and  experimental  values. 

The  elastic-viscoplastic  model  as  defined  by  Equation  4.15  does  not  explicitly 
use  E in  the  formulation,  rather  the  shear  modulus  ( G ) is  used.  The  shear  modulus 
may  be  evaluated  from  the  elasticity  relationship: 
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Figure  4-6:  Bulk  modulus  (K)  as  a function  of  mean  normal  stress 


G 


2(1  + u') 


(4.28) 


4.5  Yield  Function  Evaluation 
Recall  that  the  irreversible  deformation  rate  is  given  by: 

W(t)  \ dF 


£h  kT^  H{amn)  / do 


(4.29) 


i] 


In  the  special  case  where  the  irreversible  deformation  rate  approaches  zero, 
Equation  4.29  becomes: 


0 — kf 


Reducing  further  yields: 


W{t)  \ dF 

H{amn)  / do yj 


(4.30) 


W{t)  < H(omn), 


(4.31) 
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The  corresponding  stress  state,  amn,  belongs  to  the  stabilization  boundary. 

Such  conditions,  where  the  irreversible  deformation  rate  approaches  zero,  occur 
in  creep  experiments.  Thus,  this  boundary  can  be  experimentally  mapped  out 
through  a series  of  carefully  controlled  creep  experiments.  In  a creep  experiment, 
the  stress  state  is  kept  constant  for  a certain  interval  while  the  strains  are  con- 
tinuously measured.  When  the  deformation  over  time  practically  stops,  it  can 
be  assumed  that  a stabilized  state  has  been  reached.  Figure  4 7 depicts  how  the 
stabilization  boundary  can  be  mapped  out  experimentally  in  a (a,  e)  plane.  The 
stair-step  line  represents  the  loading  and  creep  steps  of  such  an  experiment.  In  any 
given  creep  stage,  the  stress  is  held  constant  until  the  strain  rate  approaches  zero. 
The  stabilization  boundary  is  defined  by  the  stress-strain  points  taken  at  the  end  of 
each  creep  cycle. 


£ 


Figure  4-7:  Stabilization  boundary 
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The  material  being  isotropic,  the  yield  function,  H,  depends  on  through  it’s 
invariants: 


b 

•40 

III 

(4.32) 

J2  = 7,{tra'i3)2, 

(4.33) 

Ji  = 

(4.34) 

where  I\  is  the  first  stress  invariant  while  J2  and  J3  are 
invariants  of  the  stress  deviator.  By  definition, 

1 the  second  and  third 

1,  ^ lr 

P = ^{traij)  = -I 

(4.35) 

Let  q be  defined  as: 

q = v/3  J2  = : a'ip 

(4.36) 

For  simplicity,  let  us  assume  that  the  influence  of  «J3  on  H may  be  neglected,  then 
H is  a function  of  p and  q only. 

For  the  special  case  of  a Karman  experiment,  q reduces  to  the  principal  stress 
difference,  q = cri  — <73  and  p = 4 (<7!  + 2<72)-  Let  us  assume  that  the  yield  function 
can  be  expressed  as  the  sum  of  two  terms: 

H(p,  q)  = HH(p)  + HD(p,  q), 

(4.37) 

with 

hd{p,  0)  = 0. 

(4.38) 

Hence,  Hh  and  HD  are  the  hydrostatic  and  deviatoric  contributions  respectively. 
Since  under  hydrostatic  conditions,  the  yield  function  reduces  to  HH(p),  it  may  be 
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determined  directly  from  data  obtained  in  the  hydrostatic  compression  phase  of  the 
Karman  tests.  Under  hydrostatic  loading,  (TtJ  = pS^ , the  irreversible  stress  work 
reduces  to: 


W[t)  = WH 


p(t)i‘v(t)dt 


(4.39) 


Here,  e[  is  the  irreversible  volume  strain  rate  and  tn  is  the  time  corresponding  to 
the  end  of  the  hydrostatic  loading  phase.  Values  for  W(t ) were  calculated  from 
the  experimental  data,  and  the  points  corresponding  to  the  end  of  the  creep  cycles 
plotted  as  a function  of  mean  normal  stress.  Results  from  the  150,  300,  450  and 
500  MPa  hydrostatic  tests  are  presented  in  Figure  4-8. 


Figure  4-8:  Irreversible  hydrostatic  work 


The  values  for  Wh  are  initially  nearly  zero,  which  is  expected  since  the 
material  is  initially  responding  primarily  elastically.  At  approximately  150  MPa 
mean  normal  stress,  there  is  an  upturn  in  the  data  after  which  the  Wh  increases 
almost  linearly  with  increasing  pressure.  Eventually,  as  the  pressure  is  continually 
increased,  all  of  the  microcracks  and  pores  should  close  and  the  material  will 
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achieve  a fully  compressed  state  and  respond  to  increased  hydrostatic  loading  as 
a perfectly  elastic  body.  At  this  point,  the  curve  should  flatten  out  since  with 
increasing  pressure,  further  irreversible  volumetric  work  will  not  be  possible.  This 
plateau  obviously  lies  beyond  the  hydrostatic  pressure  limit  of  the  equipment 
available  for  this  modeling  effort.  In  order  to  achieve  higher  pressure  volumetric 
compression  results,  it  may  be  necessary  to  utilize  a uniaxial  strain  compression  cell 
as  described  in  Cristescu  [8].  A cumulative  symmetric  double  sigmoidal  transition 
function  was  selected  as  providing  the  best  match  to  the  experimental  data,  as  well 
as  replicating  the  expected  plateau  for  the  fully  compressed  material  behavior.  The 
expression  is  shown  below  as  Equation  4.40 


2(4  In  ( exp 


2 p + ch 

2(4 


-2(4  In  ( exp  ( — ) + exp 


P 


+ exp 
2bft  + Ch 

2(4 


h 

dh 


+ Ch 


+ Ch 


(4.40) 


where  the  coefficient  “e4’  assures  zero  intercept,  uah"  defines  the  transition  height, 
“ bh ” the  transition  center  while  the  transition  width  is  defined  by  the  expression: 


width  = 41n(1-exp(il))“<i'‘ln(expfe)’exp(^:)) 

~dh  ln  (‘  - exp  (tf))  + <Un  (exp  (^)  “ exp  (5f))('4'41) 

Figure  4-9  displays  the  experimental  data  along  with  the  theoretical  expression 
for  Hh(p)  evaluated  using  the  parameter  values:  ah  — 12.602  MPa,  bh  = 364.272 
MPa,  Ch  = —364.935  MPa,  dh  — 30.696  MPa,  and  = 0.0028347  MPa.  These 
coefficients  provided  an  excellent  fit  to  the  data  with  an  r2  of  0.9998.  The  plateau 
level  was  selected  rather  arbitrarily  by  choosing  what  seemed  a reasonable  value 
and  including  a heavily  weighted  point  of  this  value  at  1 GPa  mean  normal  stress. 

In  order  to  fine  tune  the  model,  it  would  be  highly  desirable  to  acquire  data  at 
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confining  pressures  up  to  1 GPa,  and  in  close  proximity  to  the  two  inflection  points 
in  the  curve. 
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Figure  4-9:  Hydrostatic  yield  function 


In  a similar  manner,  the  deviatoric  yield  function  {Hp)  is  determined  from  the 
deviatoric  phase  of  the  Karman  experiments.  In  the  deviatoric  phase,  the  lateral 
confining  pressure,  cr3,  is  held  constant,  while  the  axial  stress,  cp,  is  increased  up  to 
the  point  of  material  failure.  The  state  of  stress  under  these  conditions  is  given  in 
Equation  4.42. 


a,{t) 


oi 


, t 


(4.42) 


V 


(73 


Under  the  experimental  conditions,  the  irreversible  stress  work  may  be  written: 


a3(t)i1v(t)dt+  q(t)e'\(t)dt,  t > tH 


’tH 


' tH 


W(t)  = WD(t)  = 


(4.43) 
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Wo  was  calculated  for  the  deviatoric  portions  of  each  experiment.  The  values 
for  Wo  corresponding  to  the  end  of  the  creep  cycles  were  then  plotted  as  a function 
of  q (see  Figure  4-10).  The  result  was  a family  of  curves  dependent  on  both  q and 
the  radial  confining  pressure  <r3.  At  any  given  confining  pressure  the  data  trend  is 
the  same:  i.e.  initially  linear,  eventually  curving  upward  and  asymptoting  towards 
infinity  as  the  specimen  underwent  failure.  In  order  to  capture  this  behavior,  a 
function  was  selected  which  was  a summation  of  linear,  quadratic  and  exponential 
terms.  Thus,  for  any  given  cr3,  the  Wo  vs.  q curve  was  approximated  with  the 
function: 


Hd(o 3,  q)  = d0q  + dxq2  + d2  exp  (d3q2)  - d2,  (4.44) 


where  d*  are  the  constants  for  a given  test.  The  expression  is  shown  to  be  a 
function  of  both  <r3  and  q , though  only  q is  explicitly  represented.  The  other  term 
enters  in  that  the  4 parameters  from  each  of  the  curve  fits,  must  now  be  expressed 
as  a function  of  ct3. 


Figure  4-10:  Deviatoric  yield  function 
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Next,  each  of  the  parameters  (d0,d\,d2,  and  d3)  was  plotted  as  a function  of 
radial  confining  pressure.  The  following  expressions  were  found  to  approximate  the 
variation  of  these  parameters  as  a function  of  confining  pressure. 

d0(a3)  = d0a  • (1  - exp  (d0ba3)) 

(4.45) 

d\(o3)  = dla  • (1  - exp  (dha3)) 

(4.46) 

d2(a3)  = d2a+d2b-{  1 - exp  ( d2ca3 )) 

(4.47) 

H ( ~ \ ^3 a 

d3(cr3)  - 

1 + d3a  ■ d3b  ■ a3 

(4.48) 

Values  for  all  the  parameters  given  in  Equations  4.45  through  4.48  are  pro- 
vided in  Table  4-1. 

Table  4-1:  Parameter  values  for  Ho 


^0a 

dob 

dla 

d\b 

dla 

d2h 

d2c 

d'ia 

dib 

0.0141 

-0.007512 

4.08087-10-5 

-0.00489 

1.5935-10-4 

0.05640 

-0.008213 

0.00135 

175.26 

A comparison  between  the  analytical  model  for  the  deviatoric  yield  function 
and  the  experimental  results  corresponding  to  200  MPa  confining  pressure  are 
presented  in  Figure  4-11.  The  200  MPa  experiment  was  not  utilized  in  the 
determination  of  the  the  function,  so  it  should  provide  an  excellent  test  case.  As 
can  be  seen,  correlation  between  experiment  and  theory  is  excellent. 

In  order  that  the  expression  for  HD  be  as  general  as  possible,  it  is  desirable 
that  it  be  expressed  purely  in  terms  of  the  stress  invariants.  Under  the  special 
conditions  of  the  Karman  test,  a3  = p — |.  Thus  in  general  terms,  Equation  4.44 
may  be  written: 


HD(p,  q)  = d0q  + dxq2  + d2  exp  (d3q2)  - d2. 


(4.49) 
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4.6  Strain  Rate  Orientation  Tensor 

All  elements  of  the  model  have  now  been  formulated  with  the  exception  of 
the  potential  function,  Determination  of  the  viscoplastic  potential  function 
can  be  quite  challenging  and  in  some  cases  not  possible.  As  an  alternative,  the 
concept  of  an  irreversible  strain  rate  orientation  tensor  was  introduced  by  Cazacu 
[102],  Let  us  define  N(oij)  a tensorial  function  of  the  stress  tensor  which  governs 
the  orientation  of  the  irreversible  strain  rate,  eL.  For  an  isotropic  material,  N must 
satisfy  the  invariance  requirement 

N(QalJQT)  = QN{alJ)QT  (4.50) 

where  Q represents  any  orthogonal  transformation.  From  classical  results  re- 
garding the  representation  of  isotropic  tensor  functions,  it  follows  that  N can  be 
represented  as, 


N{&ij)  — N\I  + N2ol3  + 


(4.51) 
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where  Nt  are  scalar-valued  functions  of  the  stress  invariants.  In  order  to  reduce  the 
complexity  of  the  analysis,  the  simplifying  assumption  that  the  contributions  of  N3 
and  the  third  stress  invariant  are  negligible  is  made.  Thus,  the  expression  for  N 
reduces  to: 

gl 

N((Jij)  = Ah(p,<?)/  + N2{p,q )-^,  (4.52) 

where  / is  the  identity  matrix  (for  more  details,  see  [102]).  The  irreversible  portion 
of  the  constitutive  equation  can  now  be  written: 

i'ij  = kT(l-  ^ | Ah(p,g)/  + N2(p,q)~^j  (4.53) 

For  the  special  case  of  the  triaxial  Karman  experiment,  e2  = £3.  Writing  out 
the  individual  terms  of  ek  under  these  conditions  yields: 

e[  = kT(l-  |^i (P. 9)  + ^N2(p,q)\  (4.54) 

^ = kT  ('  - ^))  • <4-55> 

Combining  Equations  4.54  and  4.55  in  order  to  isolate  the  Ah  and  N2  terms 
yields: 

el  = 3kT(l-  {Nx {p,q)}  (4.56) 

(4.57) 

Rearranging  Equations  4.56  and  4.57,  expressions  for  kTN\  and  krN2  can  be 
written  strictly  in  terms  that  are  either  measurable,  or  previously  determined  in  a 
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Karman  experiment. 


kTN\(p,  q)  = 


3(1- 


W{t) 


(4.58) 


kTN2{p,q) 


W(t)  \ 
y)/ 


(4.59) 


The  form  of  Equation  4.58  indicates  that  N\  controls  the  irreversible  volumet- 
ric response,  while  1 V2  (Equation  4.59)  will  control  the  irreversible  shear  response. 
The  ktN\  term  may  be  thought  of  as  consisting  of  a hydrostatic  and  deviatoric  part 
as  indicated  in  Equation  4.60. 


krN\(p,  q ) = 0(p)  + q)  (4.60) 

Equation  4.58  is  utilized,  along  with  data  from  the  hydrostatic  compression 
phase  of  the  experiments  to  determine  0(p).  Figure  4-12  presents  results  calcu- 
lated from  the  experiments  conducted  at  50,  200,  375,  and  500  MPa  hydrostatic 
confinement.  There  is  a high  level  of  correlation  between  the  three  monotonically 
conducted  experiments.  The  500  MPa  experiment  was  conducted  with  a number 
of  creep  cycles.  Due  to  the  nature  of  the  expression,  within  the  creep  cycle,  a nu- 
merical singularity  would  occur  as  the  material  approached  creep  stabilization  (i.e. 
W(t)  — H(oij)).  Thus  the  data  presented  for  this  experiment  has  had  the  creep 
and  unload  cycles  removed.  For  each  cycle,  following  an  initial  “ramp  up” , the  ma- 
terial appears  to  follow  the  same  surface  as  the  monotonic  experiments.  The  data 
from  the  375  MPa  experiment  increases  monotonically,  while  for  the  final  displayed 
cycle  of  the  500  MPa  experiment,  the  data  appears  to  level  out.  As  was  discussed 
in  the  derivation  of  the  hydrostatic  yield  function,  Hu,  as  the  mean  normal  stress 
is  increased,  eventually  the  material  will  reach  a fully  compacted  state  where  no 
further  plastic  deformation  is  possible.  This  state  is  indicated  by  the  plateau  on  the 
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Hh  curve  (Figure  4-9).  In  this  fully  compressed  state  e[  = 0,  thus,  the  expression 
for  0(p)  should  peak  and  then  begin  to  fall  returning  to  zero.  Assuming  the  plateau 
for  Hh  was  estimated  correctly,  e[  should  achieve  a value  of  zero  at  approximately 
600  MPa  and  remain  zero  for  p > 600 MPa.  Figure  4-13  displays  the  values  of 
s(,(p)  for  the  375  MPa  and  500  MPa  hydrostatic  compression  experiments.  Correla- 
tion between  these  two  experiments  is  remarkably  good.  The  value  for  e[,  reaches  a 
peak  value  at  approximately  290  MPa,  and  begins  to  plateau  for  a period,  and  then 
fall.  The  agreement  in  e[  between  the  two  experiments  gives  one  confidence  that 
the  leveling  of  the  0(p)  results  is  real.  It  also  lends  confidence  that  the  estimation 
of  the  plateau  in  Hh  is  reasonable,  since  the  rate  of  volumetric  compression  has 
begun  to  slow  approximately  half  way  into  the  linear  portion  of  the  function.  In 
order  to  mathematically  model  such  behavior,  it  would  likely  be  necessary  to  resort 
to  a two  part  expression,  one  for  increasing  0(p)  and  one  for  decreasing  values.  Due 
to  the  lack  of  data,  the  shape  of  the  descending  function  would  be  pure  conjecture. 
Therefore,  at  this  point,  the  decision  was  made  to  model  the  behavior  with  a single 
smooth  continuous  function  which  asymptotes  to  a constant  value  corresponding  to 
the  plateau  indicated  in  the  experimental  data.  A sigmoid  function  of  the  form: 


<Kv)  = 4>o  + 


0 

i 

1 + exp 

p-<j>2 

03 

(4.61) 


was  found  to  approximate  the  data  well,  where:  0O  = -7.163e-5,  0]  = 0.008137,  02 
= 295.866,  and  03  = 62.632.  Only  the  data  from  the  375  MPa  experiment  and  the 
450  MPa  cycle  of  the  500  MPa  hydrostatic  experiment  were  utilized  to  determine 
0o  to  03.  A comparison  of  the  model  results  to  the  experimental  data  is  given  in 
Figure  4 14. 

To  determine  xp(p,q),  Equation  4.58  along  with  the  data  from  the  deviatoric 
shear  portion  of  the  experiments  were  used.  Results  for  the  200  MPa  experiment 
are  presented  in  Figure  4-15.  Once  again,  only  the  results  from  the  monotonic 
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Figure  4-12:  Strain  rate  orientation  tensor  - kNl  hydro 

loading  cycles  are  included.  The  first  loading  cycle  shows  a clear  monotonically 
increasing  trend  in  the  data,  however,  on  subsequent  reload  cycles,  there  does  not 
appear  to  be  sufficient  spacing  between  the  creep  cycles  for  the  specimen  to  ever 
reach  the  “surface”  upon  reloading.  If  it  does  reach  the  surface,  it  is  only  very 
shortly  prior  to  initiation  of  the  subsequent  creep  phase.  On  the  final  cycle,  the 
material  ramps  up  only  slightly,  then  transitions  to  a negative  slope,  eventually 
becoming  negative.  Close  inspection  of  the  value  for  q at  the  point  kN\  crosses  the 
q axis  indicated  it  corresponds  to  the  transition  from  compressibility  to  dilatancy 
indicated  in  the  volume  strain  curve  (Figure  2-40).  The  data  as  presented  does  not 
appear  to  be  adequate  for  selecting  a suitable  function  representative  of  ip(p.q). 
What  can  be  deduced  from  the  trends  in  the  data,  and  the  characteristics  of 
Equation  4.58  is  that  the  function  should  initially  increase  with  positive  slope, 
reach  a transition  peak  where  the  slope  becomes  negative,  and  then  achieve 
negative  value  at  a point  corresponding  to  the  compressibility-dilatancy  boundary. 
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Figure  4-13:  Irreversible  volumetric  strain  rate 

Finally,  kTN2{q,p)  was  determined  by  utilizing  Equation  4.59  along  with 
data  from  the  deviatoric  shear  experiments.  Results  for  the  200  MPa  analysis  are 
presented  in  Figure  4-16.  Once  again,  after  the  first  loading  cycle,  the  data  on 
subsequent  reload  cycles  does  not  clearly  indicate  a kj'N2  surface.  It  does  appear 
however  that  a monotonically  increasing  function  would  adequately  represent  the 
material  behavior. 

Ideally,  in  order  to  adequately  determine  ip(p,q),  several  monotonic  deviatoric 
experiments  should  have  been  included  in  the  experimental  matrix.  Under  the 
assumption  that  adequate  information  is  available  within  the  experimental  data,  a 
process  was  devised  and  utilized  in  an  attempt  to  extract  the  information.  The  375 
MPa  test  was  initially  chosen  for  analysis  since  it  only  contained  one  creep-unload- 
reload  cycle.  The  process  began  by  hrst  extracting  all  of  the  data  rows  containing 
the  creep  cycles  from  the  data  set.  Ideally,  a monotonic  experiment  conducted 
under  load  control  would  have  a very  nearly  linear  principal  stress  difference  vs. 
time  loading  profile.  Simply  removing  the  creep- unload- reload  cycles,  yielded 
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Figure  4-14:  Comparison  of  data  to  model  for  kNl  hydro 

two  linear  sections  offset  by  time.  A correction  was  made  to  the  time  column 
by  subtracting  out  the  time  period  over  which  the  creep-unload-reload  cycle  was 
conducted  from  the  second  linear  portion.  This  resulted  in  a very  nearly  linear 
data  set  (see  Figure  4-17).  This  resultant  data  set  was  subsequently  curve  fit  as  a 
function  of  time.  The  fine  line  in  Figure  4-17  represents  the  curve  fit  data. 

Following  this  step,  both  the  axial  and  radial  strains  were  plotted  as  a function 
of  time.  Further  data  rows  were  subjectively  removed  in  order  to  provide  a smooth 
data  set  for  curve  fitting.  The  axial  and  radial  strains  were  plotted  as  a function  of 
time  and  a complete  set  of  monotonic  data  was  subsequently  constructed.  Figure 
4-18  contains  the  results  for  the  axial  strain.  Both  the  original  data  set,  as  well  as 
the  retained  adjusted  data  are  presented.  The  fine  line  represents  the  constructed 
monotonic  axial  strain  as  a function  of  constructed  principal  stress  difference. 

Using  this  adjusted  data  set,  values  for  ip(p , q)  and  kr N2  were  calculated.  The 
results  are  presented  in  Figures  4 19  and  4-  20.  The  new  data  set  for  <i)  cap- 
tures the  expected  features  including  the  expected  peak  and  slope  transition,  with 
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Figure  4- 15:  Strain  rate  orientation  tensor  - kNl  deviatoric,  200  MPa  experiment 

a sign  change  of  xp(p,  q)  from  positive  to  negative  corresponding  the  compressibility- 
dilatancy  transition  point.  The  kTN2  results  exhibit  the  monotonically  increasing 
function  predicted. 

Finally,  Figures  4-21  and  4-22  contain  the  results  for  xp(p,  q)  and  kTN2 
from  the  original  and  adjusted  data  sets.  The  results  in  both  cases  indicate  that 
the  adjusted  data  set  accurately  replicates  the  trends  indicated  in  the  original 
experimental  data  and  is  of  suitable  quality  for  function  evaluation.  This  same 
data  analysis  technique  was  subsequently  employed  on  all  of  the  deviatoric  shear 
experiments. 

The  results  of  the  deviatoric  xp  evaluations  from  all  adjusted  data  sets  are 
shown  in  Figure  4-23.  Positive  ip  corresponds  to  compressible  response  of  the 
material,  while  negative  xp  indicates  dilatant  behavior.  The  results  indicate  a family 
of  curves  with  clear  pressure  dependency.  Due  to  scaling,  it  is  very  difficult  to 
see  the  the  unconfined  response  in  Figure  4-23.  Figures  4-  24  and  4-25  present 
the  unconfined  results  in  detail.  From  these  two  figures,  it  can  be  concluded  that 
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Figure  4 16:  Strain  rate  orientation  tensor  - kN2,  200  MPa  experiment 

the  unconfined  compression  test  is  largely  a dilatant  event  with  only  very  little 
early  compressibility.  The  shape  of  the  compressive  portion  of  the  curve  from  the 
50  MPa  experiment  does  not  correlate  well  with  the  other  experiments.  There  is 
insufficient  data  in  this  pressure  regime  to  indicate  whether  this  is  truly  indicative 
of  the  material  behavior,  or  whether  it  is  an  experimental  anomaly.  Therefore,  the 
decision  was  made  not  to  utilize  the  compressive  data  from  this  experiment  in  the 
subsequent  parameter  fits.  In  all  the  experiments,  there  were  three  obvious  phases 
consistently  observed,  a pre-peak  compression  phase,  post-peak  compression  and  a 
dilatant  phase.  It  would  be  very  difficult  to  derive  a function  capable  of  replicating 
this  complex  behavior,  much  less  hope  to  determine  the  variation  of  the  coefficients 
of  such  a function  with  mean  pressure.  Consequently,  the  decision  was  made  to 
model  each  phase  separately.  There  are  three  points  which  might  be  thought  to 
bound  the  compressive  domain,  the  origin,  the  location  of  the  maximum  value  of 
ip,  and  the  compressibility-dilatancy  boundary.  Any  function  chosen  to  represent 
ip  must  pass  through  these  points.  Let  j3  be  defined  as  the  peak  value  of  ip  for  any 
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Figure  4-17:  Time  adjusted  loading  profile 


given  confining  pressure.  Figure  4-26  presents  the  values  of  0 as  a function  of  the 
confining  pressure,  denoted  by  x-  A rational  function  of  the  form: 


approximates  well  the  variation  of  (3  with  x for  (30  = 0.000002073,  0\  = 0.0000335, 
02  = 0.003967,  and  /?3  = 0.00000171.  The  solid  line  in  Figure  4-26  represents 
the  theoretical  predictions  obtained  using  Equation  4.62,  whereas  data  from  the 
different  experiments  are  represented  by  symbols. 

Let  7 be  defined  as  the  value  of  q at  which  the  peak  value  of  ip  occurs  for  any 
given  x-  Figure  4-27  displays  the  values  of  7 as  a function  of  x for  the  different 
Karman  experiments  conducted.  A simple  expression  of  the  form: 


was  found  to  approximate  the  data  quite  well.  The  solid  line  in  Figure  4-27 
indicates  the  theoretical  results  given  by  Equation  4.63.  The  numerical  values  of 


(4.62) 


7 (p,  q)  = 7oX71  + 72 


(4.63) 
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Figure  4-18:  Adjusted  axial  strain  results  - 375  MPa  experiment 

the  parameters  for  Equation  4.63  were  found  to  be:  70  = 7.98,  71  = 0.699,  and  72 
= 1.794. 

The  data  defining  the  compressibility-dilatancy  (C/D)  boundary  is  presented 
in  Figure  4-28.  The  data  points  were  taken  as  the  value  of  q at  which  ip  transi- 
tioned from  positive  to  negative  for  a given  value  of  An  exponential  expression 
of  the  form  given  in  Equation  4.64  was  found  to  approximate  the  data  well. 

a{p,q)  = cco  + Qi  (1  - exp[— a2x])  (4-64) 

The  solid  line  in  Figure  4-28  represents  the  theoretical  values  of  the  C/D 
boundary.  The  numerical  values  of  the  coefficients  involved  in  Equation  4.64  were 
found  to  be:  q0  = 3.596,  oq  = 590.055,  and  a2  = 0.0070659. 

A Gaussian  expression  of  the  form: 

1 / q - ip2  \ 

2\  ip3  J 


ipCi  (p,  q)  = Vb  + tA  exp 


(4.65) 
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Figure  4-19:  Strain  rate  orientation  tensor  - kNl,  375  MPa  experiment 


was  chosen  to  approximate  the  pre-peak  response  of  ip  and  was  found  to  match  the 
experimental  data  well.  The  variable  xpci  is  defined  as  the  pre-peak  behavior  of  ip. 
The  Gaussian  is  a symmetric-peak  function  where,  ip0  is  a parameter  controlling 
the  functions  y intercept,  ip\  is  related  to  the  amplitude  of  the  function,  ip2  defines 
the  location  of  the  functions  peak  value,  while  the  overall  area  is  defined  by 
\/27 }ip\ipz.  The  rational  expression: 


, / N (V'l.+V'uX)  , , 
Vd  (P,  q)  = ~~r.  — + V>i„ 


(4.66) 


(!  + iPicX) 

was  chosen  to  approximate  ip\ . The  numerical  values  of  the  parameters  were 
determined  to  be:  ip\a  = 0.00003392,  ip\h  = 0.00003598,  ip\c  = 0.003899,  and  \p\d  = 
0.00000438.  Note  that 


tpo {p,q)  = P(p,q)  - ,Pi(p,q) 


(4.67) 
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Figure  4-20:  Strain  rate  orientation  tensor  - kN2,  375  MPa  experiment 

while 


fo(p,Q)  = 7(P,$)  (4-68) 

. Thus,  both  functions  ip o and  ip2  are  fully  determined  without  the  introduction  of 
any  additional  parameters.  The  final  parameter,  ip^,  was  modeled  using  the  simple 
expression: 

^(p,  q)  = IbaX***  + ^3r,  (4.69) 

where  ip3a  = 1.95,  ip$b  = 0.79,  and  ip3c  = 13. 

The  post-peak  response  of  ip  is  given  by  Equation  4.70,  the  identical  form  of 
expression  utilized  for  the  pre-peak  response  (Equation  4.65).  This  had  the  effect  of 
minimizing  the  number  of  new  parameters  required. 
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Figure  4-21:  Comparison  of  analytical  to  experimental  results  - kNl,  375  MPa 
experiment 


^C2  (p,  q)  = "04  + V>5  exP 


1 / q-  Ip  6 y 

2 Wt  ; 


The  variable  ip c2  is  defined  as  the  post-peak  behavior  of  i/j,  where 


(4.70) 


^5  (p,g)  = g) 


(4.71) 


and 


^e(p.0)  = 7 (p,?)- 


(4.72) 


Since  'ip  must  cross  the  q axis  at  the  point  corresponding  to  the  compressibil- 
ity/dilatancy  boundary,  using  Equations  4.70  to  4.72  it  follows  that  ip4,  is  expressed 


as: 


Mp,q)  = P{p,q)  - 


P(p-q) 


1 — exp 


_1  / a(p,. 

2 V 


q)-7  (p.g)  \ ' 

) 


(4.73) 
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Figure  4-22:  Comparison  of  analytical  to  experimental  results  - kN2,  375  MPa 
experiment 

A constant  value  of  ip7  — 5000  MPa  was  determined  to  adequately  represent  all  the 
data.  Thus,  the  entire  post-peak  response  was  modeled  with  the  introduction  of 
only  one  additional  parameter. 

Now  only  the  dilatant  portion  of  ip  remains  to  be  modeled.  Looking  again  at 
Figure  4 23,  there  are  only  a limited  number  of  experiments  with  valid  dilatant  ip 
data.  There  is  an  obvious  difference  in  behavior  however  at  low  and  high  pressures. 
Recall  for  a moment  the  volume  strain  curves  from  the  deviatoric  experiments 
of  Chapter  2.  Under  uniaxial  conditions  and  low  confining  pressures,  WES5000 
exhibited  a gradual  transition  from  compressibility  to  dilatancy,  with  a significant 
portion  of  the  experiment  conducted  under  dilatant  conditions.  At  higher  pressures 
(300  MPa  and  greater),  the  onset  of  dilatancy  only  occurred  very  near  specimen 
failure.  An  analogous  response  is  again  apparent  for  the  dilatant  ip  behavior.  In 
the  unconfined  experiment,  the  dilatant  values  of  ip  display  a gradual,  monotonic 
negative  trend  with  asymptotic  behavior  near  failure  (Figure  4 24).  This  initial 
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Figure  4-23:  Strain  rate  orientation  tensor  - ip 


monotonic  behavior  corresponds  to  the  gradual  dilatant  deformation  of  the 
material.  The  higher  pressure  tests  at  375  MPa  and  450  MPa  radial  confinement 
displayed  almost  immediate  asymptotic  behavior,  corresponding  to  the  onset  of 
dilatancy  only  very  near  failure.  A function  of  the  form: 


*l>D{p,q) 


ip8{q-a(p,  q)) 

1 - f 

V9 


(4.74) 


was  chosen  to  model  the  behavior.  The  variable  ipp  is  defined  as  the  dilatant 
behavior  of  ip.  Parameter  ip8  acts  to  control  the  initial  slope  of  the  response,  while 
ipg  controls  the  asymptotic  behavior.  The  simple  linear  expression: 


MP’Q)  = 


(4.75) 


was  found  to  be  adequate  for  modeling  ip$,  where  ip8a  = — 1.1996e  ' and  ip8h  = 
— 3.04e-7.  Parameter  ip9  was  modeled  using  a rational  function  of  the  form: 


Vbfaq) 


V>9a  + fl\X 
1 + V^9CX 


(4.76) 
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Figure  4-24:  Strain  rate  orientation  tensor  - ip,  unconfined  experiment 

Close  correlation  to  the  experimental  data  was  achieved  using  the  values  ipg a = 

60,  ip9b  = 4.3166578,  and  ip9c  = 0.004298336.  All  three  phases  of  ip  have  now  been 
modeled.  Comparisons  of  the  theoretical  results  for  ip  to  the  experimental  data 
are  presented  in  Figures  4-29  through  4-31.  Figure  4-29  displays  the  compressive 
domain  of  ip  for  the  unconfined  experiment.  The  theoretical  and  experimental 
results  are  seem  to  match  almost  exactly.  Figure  4-30  shows  the  complete  range 
of  unconfined  ip  results,  both  experimental  and  theoretical.  The  theoretical  results 
again  compare  quite  well,  capturing  both  the  quasi-linear  and  asymptotic  behavior. 
Figure  4-31  presents  the  theoretical  and  experimental  results  for  the  375  MPa 
radial  confinement  experiment.  The  overall  shape  of  the  375  MPa  ip  results,  in 
particular  the  inflection  occurring  at  x ~ 200  MPa,  is  not  considered  to  be  truly 
indicative  of  material  behavior.  This  is  likely  an  artifact  of  the  creep  step  present 
in  the  actual  experiment.  Still,  the  model  adequately  represents  the  material 
behavior,  in  particular  the  peak,  slope  transition  and  subsequent  asymptotic 
dilatant  response. 
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Figure  4-25:  Strain  rate  orientation  tensor  - unconfined  experiment  low  PSD 
range 

The  complete  function  for  all  three  phases  is  presented  below  in  Equation  4.77. 

xp{p,q)  = xpCl,  q<i(p,q) 

= TpCi,  7 (p,q)  < q < Oi{p,q)  (4.77) 

= ipD,  a(p,  q)  < q 

The  kN2  term  is  the  only  constitutive  function  yet  to  be  determined.  Figure 
4-32  displays  the  global  kN2  at  various  confining  pressures,  while  Figure  4-33 
displays  only  the  data  for  q < 300  MPa.  What  is  seen  are  a family  of  monotonically 
increasing  curves  with  asymptotic  behavior  near  failure.  Closer  inspection  of  Figure 
4 33  shows  that  for  lower  levels  of  confinement  (0  to  150  MPa)  the  slope  of  kN2 
vs.  PSD  increases  along  with  increasing  confinement.  Beyond  confining  levels  of 
150  MPa  however,  the  slopes  decrease  with  increasing  confinement.  The  transition 
appears  to  take  place  for  y % 200  MPa.  This  corresponds  closely  with  the  inflection 
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Figure  4-26:  Peak  0 value  as  a function  of  x 


point  in  the  hydrostatic  yield  function  (Figure  4 9) , thus  what  is  seen  is  likely  the 
result  of  the  breaking  of  the  cementation  bonds  in  the  concrete. 

The  low  q range  (pre-asymptote,  quasi-linear)  response  was  found  to  be 
adequately  approximated  with  the  function: 


kN'jiP:  q)  =£o 91'5- 


(4.78) 


The  values  of  £0  were  determined  for  the  different  confining  conditions.  The 
results  are  displayed  as  a function  of  confinement  (x)  in  Figure  4-34.  The  data 
trend  indicates  a positively  sloped  function  which  peaks  in  value  at  approximately 
X = 200  MPa,  then  begins  to  decrease  in  value,  and  finally  plateaus  to  a constant 
value.  The  plateau  would  correspond  to  the  fully  compacted  state  of  the  material. 
Once  again,  the  behavior  observed  if  very  complex. 
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Figure  4-27:  Principal  stress  difference  at  peak  ip  as  a function  of  y 


The  pre-peak  behavior  could  be  modeled  using  a symmetric  peak  function  such 
as  a Lorentzian: 

{S0(P,  «)  = &.  + 7 ^ <4-79> 

('  + &)) 

where  £0a  controls  the  y intercept  of  the  function,  £0(>  determines  the  amplitude, 

£0c  the  center,  while  the  area  is  determined  by  7r£o6£od-  The  superscript  (1)  applied 
to  the  £0  variable  is  used  here  to  designate  the  Lorentzian  approximation  to  the 
data.  The  values  for  the  parameters  were  determined  to  be:  £0q  = — 2.1396e-5  , 
fo6  - 2.4650e~5,  £0c  = 200,  and  £0d  = 641.5082. 

Theoretical  results  are  displayed  in  Figure  4-35  along  with  the  experimental 
values.  The  theoretical  model  (£q!))  matches  the  pre-peak  data  quiet  well,  and 
doesn’t  begin  to  deviate  significantly  until  x = 375  MPa.  The  £o  values  have 
been  extrapolated  out  to  1 GPa  by  assuming  a constant  value  is  achieved  for  £0 
corresponding  to  the  plateau  in  the  hydrostatic  yield  function.  Let  £q  ' be  defined 
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Figure  4-28:  Compressibility  - dilatancy  boundary 


f<2)(p,  (4.80) 


The  resultant  £q2)  data  is  plotted  in  Figure  4 36  along  with  the  original  £o  values. 
The  post-peak  behavior  of  ^,2)  displays  monotonic  behavior  which  is  modeled  quite 
well  using  a Pulse  Cumulative  function  of  the  form: 


to\p,q)  = £oe  +6/ 


^1  — exp 


X _ £oc  ~ £oh  In  (l  ~ - £o9 

£oh 


(4.81) 


where  controls  the  y intercept  of  the  function,  £o/  defines  the  transition 
height,  £0s  defines  the  transition  center,  while  the  transition  width  is  given  by 
1.3169578969-£oh  • The  function  is  subject  to  the  constraint: 


fo2)  =0  for  \<  £o9  + £o„  In 


+ £o  c- 


(4.82) 
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Figure  4-29:  Theoretical  ^ compared  to  unconfined  experimental  data,  compressive 
domain 

The  values  of  the  parameters  in  Equation  4.81  were  determined  to  be:  £0e  = 

3.6924e-9,  £0/  = — 2.1373e~5,  £o9  = 612.7689,  and  £oh  = 441.1152.  Substituting 

(2) 

these  parameters  into  Equation  4.82  yields  the  constraint  that  Q = 0 for 
X<  271.10. 

In  order  to  implement  this  constraint,  the  concept  of  a unit  step  function  shall 
be  introduced.  Let  U[g]  define  a function  of  the  argument  g such  that: 


U\g]  = 1,  g > 0 

U\g]  = 0.  g < 0 (4.83) 


Rearranging  the  terms  in  Equation  4.80  and  incorporating  the  step  function  C7, 
yields  the  expression: 

V2 


£o  (p,q)  = Zo}  ~ U 


X ~ €og  + £oh  ln  1 — + £o„ 


«J21.  (4.84) 
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Figure  4-30:  Theoretical  1/1  compared  to  unconfined  experimental  data 

A comparison  of  the  theoretical  results  yielded  by  Equation  4.84  to  the 
experimental  data  is  shown  in  Figure  4 37.  In  general,  agreement  between  the  two 
is  excellent. 

Now  the  asymptotic  behavior  of  kN2  must  be  incorporated.  The  expression 
presented  in  Equation  4.85  was  selected  to  model  the  experimental  data. 


kN2(p , q ) = So <715  + ( 1 - exp 


9 ~ ft 
ft 


U (9  — Si) 


-ft 


(4.85) 


The  So91,5  term  is  included  to  model  the  initial  behavior  previously  discussed. 
The  concept  of  the  unit  step  function  was  again  utilized  to  “turn  on”  the  asymp- 
totic behavior,  therefore  Si  determines  the  point  at  which  the  behavior  is  initiated, 
while  parameter  £2  determines  the  asymptote.  The  exponential  term  is  introduced 
to  smooth  out  the  onset  of  the  step  function.  Parameter  S3  therefore  controls  the 
transition  to  asymptotic  behavior.  The  obvious  constraint  is  Si  < S2-  A rational 
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Figure  4-31:  Theoretical  ip  compared  to  375  MPa  experimental  data 


expression  was  again  utilized  to  model  the  £2  response. 

6 (P.«)  = %X^T  (486) 

(1  + bcX) 

The  parameter  values  which  yielded  the  best  results  are:  6n  = 61.6240,  G()  = 
2.8868,  and  £2,  = 0.0017699. 

An  Extreme  Value  Cumulative  expression  was  chosen  to  model  the  £1  behav- 
ior. Given  the  limited  data  set  available,  it  seemed  to  yield  the  best  combination  of 
accuracy,  while  achieving  the  constraint  that  G <6- 


6 (p,  q)  = 6a  exP 


— exp 


X -6c  In  (ln(2))  — Gib 

6c 


(4.87) 


The  parameters  values  for  Equation  4.87  were  evaluated  as:  6a  = 414.1011. 

6fc  = 67.8311,  and  6C  — 61.0303.  Figure  4-38  presents  the  theoretical  results 
of  G and  G,  35  well  as  the  respective  experimental  results.  The  6 data  indicates 
that  the  asymptote  initiation  point  tends  to  level  off  for  higher  values  of  x- 
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Figure  4-32:  Strain  rate  orientation  tensor  - kN2 


This  is  consistent  with  the  results  found  for  the  C/D  boundary  behavior  (see 
Figure  4-28).  There  is  no  experimental  data  at  the  lower  confining  pressures  to 
determine  whether  the  low-y  response  is  representative,  however  the  results  appear 
reasonable. 

Finally,  the  findings  for  £3  again  indicated  the  possibility  of  a ‘‘peak”  type 
response  similar  to  £0-  However,  there  was  insufficient  valid  data  to  fully  map 
out  this  behavior.  Further  study  showed  that  the  accuracy  of  results  yielded  by 
Equation  4.85  appeared  to  be  rather  insensitive  to  the  value  of  <^2  provided  it  was 
of  the  proper  order  of  magnitude.  An  average  value  found  from  the  y = 0,  150  and 
375  MPa  results  was  consequently  used  (£3  = 1881464). 

Figures  4-39,  4-40,  and  4-41  present  the  kN2  experimental  data  along  with 
the  respective  theoretical  results.  In  each  case,  correlation  between  theory  and 
experiment  is  excellent. 
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Figure  4-33:  Strain  rate  orientation  tensor  - kN2,  low  PSD  range 

4.7  Comparison  of  Theoretical  to  Experimental  Results 
Analytic  simulations  were  run  for  three  different  cases,  uniaxial  compression, 
monotonic  hydrostatic  compression  up  to  375  MPa  mean  normal  stress,  and 
a monotonic  deviatoric  case  with  375  MPa  confining  pressure.  For  each  case, 
comparisons  are  made  to  experimental  data.  Figure  4-42  presents  both  theoretical 
and  experimental  results  for  the  uniaxial  compression  case.  The  model  matches 
the  data  quite  well  and  captures  both  the  compressible  and  dilatant  nature  of  the 
material  as  well  as  accurately  predicting  material  failure. 

Results  for  the  375  MPa  hydrostatic  compression  case  are  presented  in  Figure 
4-43.  Correlation  between  experimental  and  analytic  results  is  excellent,  with  the 
analytic  results  literally  falling  between  the  experimental  axial  and  radial  results. 

Finally  in  the  third  case,  the  deviatoric  analytic  and  experimental  results  are 
presented  for  375  MPa  confinement.  Again,  there  is  excellent  agreement  between 
theoretical  and  experimental  axial,  radial  and  volumetric  strains.  Compressibility, 
dilatancy  and  failure  are  all  accurately  simulated. 


179 


5C  (MPa) 


Figure  4-34:  Parameter  £0  as  a function  of  \ 
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Figure  4-35:  Parameters  £o  and  1 as  a function  of  \ 
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Figure  4-37:  Theoretical  £o  compared  to  experimental  data 
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Figure  4-38:  Parameters  £1  and  comparison  of  experimental  and  theoretical 
results 
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Figure  4-39:  Experimental/ theoretical  comparison  of  kN 2,  \ = 0 MPa 
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Figure  4-40:  Experimental/theoretical  comparison  of  k'N2,  x = 150  MPa 


Figure  4-41:  Experimental/theoretical  comparison  of  kN2,  x = 375  MPa 
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Figure  4-42:  Experimental/theoretical  comparison  of  uniaxial  compression  test 
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Figure  4-43:  Experimental/theoretical  comparison  of  375  MPa  hydrostatic  com- 
pression test 
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Figure  4-44:  Experimental/theoretical  comparison  of  375  MPa  deviatoric  compres- 
sion test 


CHAPTER  5 
CONCLUSIONS 

A database  of  mechanical  response  of  mortar  and  concrete  WES5000  was 
generated.  The  mortar  studied  contains  masonry  sand  but  no  large  aggregates. 
The  mortar  experimental  effort  consisted  of  a series  of  quasi-static  and  dynamic 
tests  at  several  strain  rates,  under  both  unconfined  and  confined  conditions. 
Confinement  levels  utilized  were  1.72  MPa  ,3.45  MPa  , and  6.89  MPa.  The  quasi- 
static experiments  were  conducted  at  a strain  rate  of  approximately  10-6/s. 

Same  size  mortar  specimens  were  tested  under  both  unconfined  and  confined 
conditions  in  the  University  of  Florida  7.62  cm  diameter  split  Hopkinson  pressure 
bar  (UFSHPB)  at  three  different  loading  rates  in  the  range  60/s  to  about  160/s. 
The  range  of  strain  rates  was  bound  by  the  capacity  of  the  UFSHPB  gas  gun 
firing  chamber.  The  highest  strain  rates  corresponded  to  a striker  bar  velocity 
of  approximately  14  m/s  obtained  by  use  of  the  0.762m  striker  bar  and  the 
chamber  gun  pressurized  to  3.45  MPa.  The  lowest  dynamic  strain  rates  achieved 
corresponded  to  a chamber  pressure  of  1.39  MPa  (striker  velocity  7.75  m/s). 
Experiments  at  a mid-range  velocity  of  11  m/s,  corresponded  to  2.41  MPa  firing 
pressure  were  also  conducted.  The  data  obtained  in  such  tests  are  of  fundamental 
interest  since  they  allow  the  characterization  of  the  combined  effects  of  confining 
pressure  and  strain  rate  on  the  deformation  and  strength  of  the  material.  The 
confined  dynamic  tests  showed  similar  stress-strain  response  as  the  quasi-static 
tests.  Dilatancy  has  been  observed  at  high  levels  of  the  principal  stress  difference 
for  both  dynamic  and  static  conditions.  The  unconfined  dynamic  compressive 
strengths  are  approximately  double  that  of  the  quasi-static  compressive  strengths. 
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Most  of  the  specimens  showed  very  little  damage  post  test  other  than  some 
chipping  around  the  top  edges,  most  likely  due  to  localized  tensile  effects. 

The  concrete  selected  for  this  research,  WES5000,  was  first  formulated  at 
the  U.S.  Army  Corps  of  Engineer  Research  and  Development  Center  Waterways 
Experiment  Station  (WES,  hence  WES500)  as  a standard  test  concrete  for  phe- 
nomenology studies.  The  mechanical  response  of  WES5000  under  a variety  of 
loading  conditions  was  investigated.  Of  particular  interest  was  to  establish  the 
influence  of  pressure  on  strength,  elastic  properties,  compressibility  and  dila- 
tant  behavior.  A test  matrix  and  test  protocol  for  studying  quasi-static  triaxial 
compression  behavior  of  WES5000  was  proposed.  A step-wise  loading-creep-unload- 
reload  procedure,  which  allowed  accurate  separation  between  elastic  and  viscous 
effects  on  the  material  behavior  was  followed. 

Hydrostatic  tests  up  to  a pressure  of  0.5  GPa  allowed  for  the  accurate  de- 
termination of  the  dependence  of  the  bulk  modulus  on  pressure  and  the  correct 
estimation  the  material’s  compaction  properties  when  subjected  to  pressures  in 
the  range  encountered  in  dynamic  events.  Under  confined  quasi-static  conditions, 
the  material  exhibited  hardening  behavior  up  to  failure.  Both  compressibility  and 
dilatancy  regimes  of  the  volumetric  behavior  were  observed,  the  dilatancy  threshold 
being  highly  dependent  on  the  level  of  confinement.  For  confining  pressures  above 
50  MPa,  transition  from  compressibility  to  dilatancy  is  apparent  only  close  to 
failure.  To  quantify  the  influence  of  loading  rate  on  the  strength  of  WES5000,  the 
dynamic  strengths  were  normalized  by  the  quasi-static  strengths  corresponding  to 
the  same  level  of  confinement.  The  dynamic  increase  factor  (DIF)was  plotted  as  a 
function  of  the  log10(e)  (see  Figure  3-63).  It  can  be  concluded  that  under  uncon- 
fined conditions,  WES5000  follow  the  same  trends  as  mortar  and  those  reported  in 
the  literature  [40,  41,  42]  for  other  cementitious  materials.  However,  under  confined 
conditions  the  DIF  values  for  WES5000  are  significantly  lower.  The  WES5000 
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generally  exhibited  far  more  cracking  under  similar  loading  conditions  than  was 
observed  in  mortar.  A new  elastic/ viscoplastic  model  that  captures  compressibil- 
ity and  dilatancy,  as  well  as  strain  rate  effects  has  been  developed.  As  a general 
framework,  the  elastic-viscoplastic  formulation  of  Cristescu  was  chosen  due  to  its 
apparent  capability  to  capture  the  behavior  of  interest  and  the  fact  that  there  are 
no  a priori  limitations  or  restrictions  regarding  the  specific  expressions  of  the  yield 
function  and  viscoplastic  potential.  However,  it  was  found  that  in  order  to  capture 
the  peculiarities  of  the  volumetric  response,  namely  the  reduced  dilatancy  at  very 
high  pressure,  the  hypothesis  of  existence  of  a viscoplastic  potential  has  to  be 
abandoned.  A new  flow  rule  was  proposed  for  WES5000.  Comparison  between  the 
model  predictions  and  data  showed  that  the  proposed  model  describes  with  very 
good  accuracy  the  high-pressure  behavior  of  concrete.  An  obvious  recommendation 
for  future  work  is  the  application  of  the  model  to  describe  high-velocity  penetration 


into  concrete. 
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